Pytanie:
Optymalizacja geometrii: zależność od geometrii początkowej
Kinformationist
2016-08-08 19:52:04 UTC
view on stackexchange narkive permalink

Można by pomyśleć, że początkowa konfiguracja nie powinna znacząco wpływać na ostateczną zoptymalizowaną geometrię. Zaczynając od dwóch ugrupowań, ostateczna „zoptymalizowana” geometria nie jest wyjątkowa i zależy w dużym stopniu od tego, w jaki sposób zostaną one umieszczone na początku. Czy jest sposób, aby tego uniknąć? W przeciwnym razie nigdy nie można było być pewnym, że skończyło się na globalnym minimum.

Dwa odpowiedzi:
Wildcat
2016-08-08 20:01:28 UTC
view on stackexchange narkive permalink

Zaczynając od dwóch reszt, ostateczna „zoptymalizowana” geometria nie jest unikalna i zależy w dużym stopniu od tego, w jaki sposób zostaną umieszczone na początku. Czy jest sposób, aby tego uniknąć?

Tak i nie. Nie, ponieważ zgodnie z przewidywaniami OP, wychodząc od geometrii wejściowej, algorytm optymalizacji rzeczywiście znajdzie najbliższe (z algorhmowego punktu widzenia) minimum lokalne i niekoniecznie globalne. I tak, ponieważ przeprowadzając obszerną analizę konformacyjną, można w zasadzie znaleźć struktury, które odpowiadają wszystkim możliwym minimom lokalnym, a następnie zidentyfikować minima globalne.

Operacyjne słowa na tym ostatnim zdaniu są *** w zasadzie ***
@F'x, smutne, ale prawdziwe. W praktyce, nawet w przypadku stosunkowo małych cząsteczek, obszerna analiza konformacyjna wymaga zbyt wiele czasu (CPU), tak więc często trzeba uciekać się do pomocy „intuicji chemicznej”, aby obniżyć koszty lub nawet całkowicie uniknąć analizy konformacyjnej.
@Wildcat - Nie jestem pewien, czy powiedziałbym, że „obszerna analiza konformacyjna” wymaga zbyt wiele czasu. Wyczerpujące (systematyczne) poszukiwania są wykładnicze pod względem liczby obracalnych wiązań, ale istnieje wiele strategii wyszukiwania konformerów, które są pragmatyczne.
@GeoffHutchison,, ale takie pragmatyczne, niesystematyczne strategie (wykorzystujące np. Losowe pobieranie próbek lub „intuicję chemiczną”) nie mogły zapewnić 100% gwarancji znalezienia globalnego minimum.
Nie, ale „Nie zamierzam robić * żadnego * próbkowania konformacyjnego, ponieważ nie mam 100% gwarancji znalezienia globalnych minimów” wydaje się nieco ekstremalne. Wypróbowanie * jakiejś * strategii jest lepsze niż nic i może przynajmniej wskazać, czy interesująca nas właściwość w dużym stopniu zależy od geometrii konformera.
Twój pierwszy komentarz sugeruje, że możesz uniknąć analizy konformacyjnej. Mówię, że osobiście nie polecam unikania samplowania konformerów i istnieją pragmatyczne rozwiązania.
Wypróbowanie jakiejś strategii jest zdecydowanie lepsze. I myślę, że w odpowiedzi można również wspomnieć, że niektóre właściwości nie są tak wrażliwe na małe zmiany w geometrii, więc znalezienie globalnego minimum może nie być tak ważne w pierwszej kolejności.
I nie, nie miałem na myśli tego, że należy unikać analizy konformacyjnej. Wręcz przeciwnie, zawsze to robię osobiście (i robię to tak systematycznie, jak to tylko możliwe). W moim pierwszym komentarzu opisałem scenariusz (zastąpienie wyszukiwania konformacyjnego „intuicją chemiczną”), którego nie lubię i nie popieram.
Można by pomyśleć, że po latach rozwoju pakiety obliczeniowe, takie jak Gaussian, musiały uwzględniać sposoby automatycznego unikania lokalnych minimów podczas optymalizacji. Konieczność ręcznej analizy konformacyjnej lub intuicji chemicznej lub jakiejkolwiek interwencji użytkownika można zminimalizować dzięki silnym algorytmom optymalizacji. Moje pytanie dotyczyło konkretnych słów kluczowych lub technik w pakiecie Gaussa, które obejmują silniejszą optymalizację.
@ChemInformationist - Gaussian na pewno nie. O ile mi wiadomo, * nie * standardowy pakiet obliczeniowy zawiera wbudowane wszelkiego rodzaju globalne strategie optymalizacji. W wielu przypadkach ludzie * chcą * lokalnych minimów, ponieważ mają intuicję chemiczną lub cząsteczka jest sztywna (np. Benzen, woda). Co więcej, w chemii istnieje * wiele * takich technik i trudno jest zastosować podejście uniwersalne dla wszystkich ze względu na różnorodność układów molekularnych i dużą liczbę stopni swobody.
@ChemInformationist Nie ma ogólnego algorytmicznego sposobu sprawdzenia, czy dane minimum lokalne jest globalne, czy nie; pierwsze dwa zdania odpowiedzi Deathbreath są nieuchronnie poprawne. Wszystko, co możesz zrobić, to wymyślić metody, które dadzą ci pewność, że zbadałeś „wystarczająco” przestrzeń konformacji. Jednym z interesujących systemów badanych tą metodą jest [periodane] (http://dx.doi.org/10.1002/qua.21322).
Deathbreath
2016-08-17 18:34:25 UTC
view on stackexchange narkive permalink

Nie jest możliwe wyprowadzenie algorytmu, który znalazłby optymalne globalne niezależnie od początkowej struktury. Jest to ogólny fakt dotyczący optymalizacji bez wypukłości. Potencjalne powierzchnie energii stabilnych związków nie są globalnie wypukłe, ponieważ energia „dysocjacji” dowolnych dwóch składników jest skończona przy nieskończonej separacji.

Chociaż istnieją tak zwane "globalne" metody optymalizacji, żadna metoda, która nie przeszukuje każdego punktu przestrzeni poszukiwań, nie może generalnie zagwarantować znalezienia optymalnego globalnego (tj. zawsze jest funkcja algorytmowi nie udaje się znaleźć optymalnego globalnego lub odwiedza każdy punkt domeny).

To powiedziawszy, metody optymalizacji globalnej przypuszczalnie zmniejszają czułość wyniku do punktu początkowego. Jest wiele tych, z których niektóre są deterministyczne (lista zdecydowanie nie jest wyczerpująca, ponieważ brakuje w niej mojej ulubionej metody uzupełniania poziomów :-).

Na koniec chciałbym zwrócić uwagę, że konkretne sformułowanie problemu wpływa na to, jak dobrze działa algorytm, np. wybrany układ współrzędnych lub sposób egzekwowania ograniczeń. Brak globalnej optymalizacji jest w rzeczywistości przydatny, jeśli ma być utrzymana pewna topologia wiązania; w przeciwnym razie t-butan zawsze optymalizowałby się na przykład do n-butanu.



To pytanie i odpowiedź zostało automatycznie przetłumaczone z języka angielskiego.Oryginalna treść jest dostępna na stackexchange, za co dziękujemy za licencję cc by-sa 3.0, w ramach której jest rozpowszechniana.
Loading...