Sensitivity of Intra- and Intermolecular Interactions of Benzo[h]quinoline from Car–Parrinello Molecular Dynamics and Electronic Structure Inspection

Jarosław J. Panek,Joanna Zasada,Bartłomiej M. Szyja,Beata Kizior,Aneta Jezierska
DOI: https://doi.org/10.3390/ijms22105220
IF: 5.6
2021-05-14
International Journal of Molecular Sciences
Abstract:The O-H...N and O-H...O hydrogen bonds were investigated in 10-hydroxybenzo[h]quinoline (HBQ) and benzo[h]quinoline-2-methylresorcinol complex in vacuo, solvent and crystalline phases. The chosen systems contain analogous donor and acceptor moieties but differently coupled (intra- versus intermolecularly). Car–Parrinello molecular dynamics (CPMD) was employed to shed light onto principle components of interactions responsible for the self-assembly. It was applied to study the dynamics of the hydrogen bonds and vibrational features as well as to provide initial geometries for incorporation of quantum effects and electronic structure studies. The vibrational features were revealed using Fourier transformation of the autocorrelation function of atomic velocity and by inclusion of nuclear quantum effects on the O-H stretching solving vibrational Schrödinger equation a posteriori. The potential of mean force (Pmf) was computed for the whole trajectory to derive the probability density distribution and for the O-H stretching mode from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling and nuclear quantum effects. The electronic structure changes of the benzo[h]quinoline-2-methylresorcinol dimer and trimers were studied based on Constrained Density Functional Theory (CDFT) whereas the Electron Localization Function (ELF) method was applied for all systems. It was found that the bridged proton is localized on the donor side in both investigated systems in vacuo. The crystalline phase simulations indicated bridged proton-sharing and transfer events in HBQ. These effects are even more pronounced when nuclear quantization is taken into account, and the quantized Pmf allows the proton to sample the acceptor area more efficiently. The CDFT indicated the charge depletion at the bridged proton for the analyzed dimer and trimers in solvent. The ELF analysis showed the presence of the isolated proton (a signature of the strongest hydrogen bonds) only in some parts of the HBQ crystal simulation. The collected data underline the importance of the intramolecular coupling between the donor and acceptor moieties.
What problem does this paper attempt to address?
This paper aims to solve the following problems: 1. **Sensitivity of intramolecular and intermolecular hydrogen bonds**: Investigate the O - H...N and O - H...O hydrogen bonds in 10 - hydroxybenzo[h]quinoline (HBQ) and the benzo[h]quinoline - 2 - methylresorcinol complex in different phases (gas phase, solvent phase, and crystal phase). These hydrogen bonds play a crucial role in the molecular self - assembly process, but their sensitivity to environmental changes has not been fully understood. 2. **Hydrogen - bond dynamics and its influence on spectral properties**: Through Car - Parrinello molecular dynamics (CPMD) simulations, study the dynamic behavior of hydrogen bonds and explore how these behaviors affect infrared (IR) spectral features. In particular, the influence of hydrogen - bond formation on the O - H stretching vibration is studied, including frequency shift, bandwidth change, etc. 3. **Influence of quantum effects**: Introduce quantum effects (such as by solving the nuclear Schrödinger equation) to more accurately describe the dynamic behavior of hydrogen bonds. This includes calculating the potential of mean force (Pmf) to obtain the probability density distribution of protons in hydrogen bonds. 4. **Electron - structure changes**: Use the constrained density functional theory (CDFT) and the electron localization function (ELF) methods to study the changes in the electron structure during the formation of hydrogen bonds, especially charge transfer and electron - density redistribution. ### Specific objectives - **Analyze the static and dynamic properties of hydrogen bonds**: Through CPMD simulations, study the changes in geometric parameters and dynamic behaviors of hydrogen bonds in different phases. - **Study spectral properties**: Analyze the influence of hydrogen bonds on infrared spectra through atomic - velocity power - spectrum analysis and correct it in combination with quantum effects. - **Calculate the free - energy curve**: Through Pmf calculations, reveal the energy landscape of hydrogen - bond formation, especially the migration behavior of protons in hydrogen bonds. - **Electron - structure analysis**: Use CDFT and ELF methods to study the changes in electron density during the formation of hydrogen bonds and evaluate the strength of hydrogen bonds. ### Specialties of the research system The paper selects two systems with similar donor and receptor groups but different coupling methods: - **10 - hydroxybenzo[h]quinoline (HBQ)**: It has an intramolecular hydrogen bond. - **Benzo[h]quinoline - 2 - methylresorcinol complex**: Forms hydrogen bonds through intermolecular interactions. This selection enables the study to compare the behavioral differences between intramolecular and intermolecular hydrogen bonds in different environments, thereby better understanding the mechanism of action of hydrogen bonds in molecular self - assembly. ### Expected results Through the above research, the author hopes to: - Reveal the dynamic behavior of hydrogen bonds in different phases and its influence on spectral properties. - Evaluate the importance of quantum effects in describing the behavior of hydrogen bonds. - Understand the influence of electron - structure changes on the strength of hydrogen bonds. - Provide new insights for designing new benzo[h]quinoline compounds with specific photochemical/photophysical properties. In summary, this paper aims to deeply understand the mechanism of action of hydrogen bonds in complex molecular systems through multi - scale simulations and theoretical analysis, and provide theoretical support for designing new functional materials.