Implementação e desenvolvimento de algoritmo eficiente para deformação intramolecular com o método Monte Carlo (2018)
- Authors:
- Autor USP: CEZAR, HENRIQUE MUSSELI - IF
- Unidade: IF
- Sigla do Departamento: FGE
- Subjects: FÍSICA MOLECULAR; FÍSICA COMPUTACIONAL; SOLVATAÇÃO; ESPECTROSCOPIA
- Agências de fomento:
- Language: Português
- Abstract: Com o avanço do poder computacional nas últimas décadas, a modelagem molecular de problemas em diversas áreas se tornou mais acessível, sendo hoje uma ferramenta fundamental para o entendimento de diversos processos. Em especial, simulações moleculares com campos de força clássicos vem sendo importante para a amostragem de propriedades termodinâmicas, para a determinação de estruturas e população de confôrmeros, e seleção de configurações para utilização com métodos que combinam mecânica quântica com mecânica molecular. Os principais métodos de simulação atualmente utilizados são a dinâmica molecular (MD, do inglês molecular dynamics) e o Monte Carlo (MC). Ambos os métodos são a princípio equivalentes quando o objetivo é a amostragem configuracional, tendo a MD a vantagem de permitir a análise da dinâmica e evolução temporal, e o MC a vantagem de poder gerar configurações de forma probabilística, sem a necessidade de seguir um caminho sobre a superfície de energia potencial, o que pode resultar em uma amostragem mais eficiente. Contudo, não há ainda uma metodologia de MC que possa ser considerada eficiente e bem estabelecida para a amostragem dos graus de liberdade internos de moléculas com complexidade arbitrária. Visando avançar no desenvolvimento de métodos que trabalham nesse sentido, neste trabalho o método apresentado por Shah e Maginn [1] foi implementado e aprimorado. No método, a molécula é fragmentada em partes menores, formadas de graus de liberdade rígidos, que não variam drasticamente durante a simulação. Esses fragmentos por sua vez são conectadas por graus de liberdade maleáveis, os diedros e termos não ligados da energia. Durante a simulação a molécula tem fragmentos apagados, e então é reconstruída utilizando um esquema de Configurational Bias Monte Carlo (CBMC).A contribuição deste trabalho para a metodologia consiste em generalizar os tipos de fragmentos possíveis dentro do método, a simplificação do critério de aceitação e extensões com vieses adicionais, como é o caso do viés no potencial eletrostático e no critério de aceitação da amostragem preferencial com o CBMC. A validade da implementação do método foi avaliada através de simulações em dois sistemas simples: o octano e o 1,2-dicloroetano. Comparando a amostragem e os resultados obtidos para as populações de confôrmeros com resultados de MD e experimentais ou obtidos com outros métodos de MC da literatura foi possível verificar que a implementação reproduz os resultados esperados. Além disso, o equilíbrio conformacional da molécula de óxido mesityl (MOx) que possui duas conformações: syn e anti, foi estudado. Graças a barreira de potencial entre as duas conformações de cerca de 10 kcal/mol, a MD não é capaz de realizar uma amostragem ergódica, enquanto o CBMC realiza a amostragem sem problemas. Na simulação CBMC, a inversão da população dominante em fase gasosa e em água foi observada. Simulações do MOx em acetonitrila e metanol mostraram que a população de confôrmeros anti de fato aumenta conforme a polaridade do solvente. Entretanto, devido ao estiramento da ligação C = O do MOx em metanol devido as ligações de hidrogênio, a contribuição conformacional à posição do máximo da banda de absorção não segue a tendência da polaridade, tendo um deslocamento para o azul maior na acetonitrila do que em metanol. O estiramento da ligação C = O só pode ser observado graças a introdução da amostragem de fragmentos com deslocamentos Cartesianos e algoritmo de Metropolis ao método CBMC original. Esse efeito devido ao estiramento é compensado pelo efeito solvente, de modo que a contribuição total à energia de transição segue a polaridade do solvente.
- Imprenta:
- Data da defesa: 27.09.2018
-
ABNT
CEZAR, Henrique Musseli. Implementação e desenvolvimento de algoritmo eficiente para deformação intramolecular com o método Monte Carlo. 2018. Tese (Doutorado) – Universidade de São Paulo, São Paulo, 2018. Disponível em: http://www.teses.usp.br/teses/disponiveis/43/43134/tde-19102018-165856/. Acesso em: 26 abr. 2024. -
APA
Cezar, H. M. (2018). Implementação e desenvolvimento de algoritmo eficiente para deformação intramolecular com o método Monte Carlo (Tese (Doutorado). Universidade de São Paulo, São Paulo. Recuperado de http://www.teses.usp.br/teses/disponiveis/43/43134/tde-19102018-165856/ -
NLM
Cezar HM. Implementação e desenvolvimento de algoritmo eficiente para deformação intramolecular com o método Monte Carlo [Internet]. 2018 ;[citado 2024 abr. 26 ] Available from: http://www.teses.usp.br/teses/disponiveis/43/43134/tde-19102018-165856/ -
Vancouver
Cezar HM. Implementação e desenvolvimento de algoritmo eficiente para deformação intramolecular com o método Monte Carlo [Internet]. 2018 ;[citado 2024 abr. 26 ] Available from: http://www.teses.usp.br/teses/disponiveis/43/43134/tde-19102018-165856/ - Implementação do método Parallel Tempering Monte Carlo para o estudo de propriedades termodinâmicas de nanoclusters
- Thermodynamic properties of 55-atom Pt-based nanoalloys: Phase changes and structural effects on the electronic properties
- Understanding the absorption spectrum of mesityl oxide dye in solvents of different polarities
- Solvent effect on the syn/anti conformational stability: A comparison between conformational bias Monte Carlo and molecular dynamics methods
- N2 adsorption on carbon nanomaterials: Validity of Force Fields and Effects on Electronic Structure
- DICE: A Monte Carlo Code for Molecular Simulation Including the Configurational Bias Monte Carlo Method
- A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution
- Unraveling the acid–base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand
How to cite
A citação é gerada automaticamente e pode não estar totalmente de acordo com as normas