{"title":"Series expansion of a scalable Hermitian excitonic renormalization method.","authors":"Marco Bauer, Andreas Dreuw, Anthony D Dutoi","doi":"10.1063/5.0238707","DOIUrl":"https://doi.org/10.1063/5.0238707","url":null,"abstract":"<p><p>Utilizing the sparsity of the electronic structure problem, fragmentation methods have been researched for decades with great success, pushing the limits of ab initio quantum chemistry ever further. Recently, this set of methods has been expanded to include a fundamentally different approach called excitonic renormalization, providing promising initial results. It builds a supersystem Hamiltonian in a second-quantized-like representation from transition-density tensors of isolated fragments, contracted with biorthogonalized molecular integrals. This makes the method fully modular in terms of the quantum chemical methods applied to each fragment and enables massive truncation of the state-space required. Proof-of-principle tests have previously shown that an excitonically renormalized Hamiltonian can efficiently scale to hundreds of fragments, but the ad hoc approach to building the Hamiltonian was not scalable to larger fragments. On the other hand, initial tests of the originally proposed modular Hamiltonian build, presented here, show the accuracy to be poor on account of its non-Hermitian character. In this study, we bridge the gap between these with an operator expansion that is shown to converge rapidly, tending toward a Hermitian Hamiltonian while retaining the modularity, yielding an accurate, scalable method. The accuracy is tested here for a beryllium dimer. At distances near equilibrium and longer, the zeroth-order method is comparable to coupled-cluster singles, doubles, and perturbative triples and the first-order method is comparable to full configuration interaction (FCI). The second-order method agrees with FCI for distances well up the inner repulsive wall of the potential. Deviations occurring at shorter bond distances are discussed along with approaches to scaling to larger fragments.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142949430","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spectroscopic study of energy transfer in collisions between vibrational excited H2 and CO2.","authors":"Chu Qin, Xiaofang Zhao, Shuying Wang, Yifan Shen","doi":"10.1063/5.0239602","DOIUrl":"https://doi.org/10.1063/5.0239602","url":null,"abstract":"<p><p>The collisional energy transfer between vibrational excited H2(1, 7) and CO2 was investigated by exciting H2 to a vibrational excited state of v = 1, J = 7 by the stimulated Raman scattering technique. The coherent anti-Stokes Raman spectroscopy (CARS) technique determined that H2 was excited to the H2(1, 7) state. Varying the cuvette temperature, the number of H2(1, 7) particles was found to increase with the increase in H2 molar ratio α by scanning the intensity of the CARS spectrum, with peaks at different α at a temperature of 363 ± 15 K, but the peak temperature was not sensitive to α. Scanning CARS spectra after H2 collisions yielded that the energies of the excited states of H2 were mainly distributed in the vibrational and translational states, proving that the collisions between the excited states of H2 were linear collisions. The collisional transfer rate coefficients of H2(1, 7) and CO2 were obtained by fitting the Stern-Volmer equation as kv(H2) = (2.89 ± 0.30) × 10-13 cm3 s-1 and kv(CO2) = (8.23 ± 0.42) × 10-13 cm3 s-1. Exciting H2 to different states, it was found that the collisional transfer rate coefficient of CO2 was less affected by the energy of the vibrational excited H2. The rotational temperature was obtained from the Boltzmann distribution of the rotational dynamics, and it was found that the rotational temperature of CO2(0000, J) was about 3.4 times higher than that of CO2(0001, J), which proved that the energy of the vibrational excited H2 was mainly allocated to the higher rotational state of CO2(0000).</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142949432","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Photon-mediated energy transfer between molecules and atoms in a cavity: A numerical study.","authors":"Jun Zhang, Shaohong Wang, Mengdi Guo, Xin-Ke Li, Yong-Chen Xiong, Wanghuai Zhou","doi":"10.1063/5.0242420","DOIUrl":"https://doi.org/10.1063/5.0242420","url":null,"abstract":"<p><p>The molecular energy transfer is crucial for many different physicochemical processes. The efficiency of traditional resonance energy transfer relies on dipole-dipole distance between molecules and becomes negligible when the distance is larger than ∼10 nm, which is difficult to overcome. Cavity polariton, formed when placing molecules inside the cavity, is a promising way to surmount the distance limit. By hybridizing a two-level atom (TLA) and a lithium fluoride (LiF) molecule with a cavity, we numerically simulate the reaction process and the energy transfer between them. Our results show that the TLA can induce a deep potential well, which can be seen as a replica of the potential energy surface of bare LiF, acting as a reservoir to absorb/release the molecular kinetic energy. In addition, the energy transfer shows a molecular nuclear kinetic energy dependent behavior, namely, more nuclear kinetic energy igniting more energy transfer. These findings show us a promising way to manipulate the energy transfer process within the cavity using an intentional TLA, which can also serve as a knob to control the reaction process.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142949429","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Thomas Schnappinger, Cyril Falvo, Markus Kowalewski
{"title":"Disentangling collective coupling in vibrational polaritons with double quantum coherence spectroscopy.","authors":"Thomas Schnappinger, Cyril Falvo, Markus Kowalewski","doi":"10.1063/5.0239877","DOIUrl":"10.1063/5.0239877","url":null,"abstract":"<p><p>Vibrational polaritons are formed by strong coupling of molecular vibrations and photon modes in an optical cavity. Experiments have demonstrated that vibrational strong coupling can change molecular properties and even affect chemical reactivity. However, the interactions in a molecular ensemble are complex, and the exact mechanisms that lead to modifications are not fully understood yet. We simulate two-dimensional infrared spectra of molecular vibrational polaritons based on the double quantum coherence technique to gain further insight into the complex many-body structure of these hybrid light-matter states. Double quantum coherence uniquely resolves the excitation of hybrid light-matter polaritons and allows one to directly probe the anharmonicities of the resulting states. By combining the cavity Born-Oppenheimer Hartree-Fock ansatz with a full quantum dynamics simulation of the corresponding eigenstates, we go beyond simplified model systems. This allows us to study the influence of self-polarization and the response of the electronic structure to the cavity interaction on the spectral features even beyond the single-molecule case.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7617315/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142894796","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
T Oppelstrup, L G Stanton, J O B Tempkin, T N Ozturk, H I Ingólfsson, T S Carpenter
{"title":"Anisotropic interactions for continuum modeling of protein-membrane systems.","authors":"T Oppelstrup, L G Stanton, J O B Tempkin, T N Ozturk, H I Ingólfsson, T S Carpenter","doi":"10.1063/5.0237408","DOIUrl":"https://doi.org/10.1063/5.0237408","url":null,"abstract":"<p><p>In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein-lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS-RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein-membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142948448","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Chaima Hammami, Conner Penson, Mehdi Adrien Ayouz, Viatcheslav Kokoouline
{"title":"Electronic states of the N2+ ion dissociating to the four lowest dissociation limits: Energies, transition dipole moments, and Einstein coefficients.","authors":"Chaima Hammami, Conner Penson, Mehdi Adrien Ayouz, Viatcheslav Kokoouline","doi":"10.1063/5.0245622","DOIUrl":"https://doi.org/10.1063/5.0245622","url":null,"abstract":"<p><p>This study presents Born-Oppenheimer energies and transition dipole moments of the 36 lowest electronic states of the N2+ ion as a function of internuclear distance in the interval between 1.5 and 10 bohrs obtained in first-principles calculations. The electronic states are of the total electronic spin S = 1/2, 3/2, and 5/2, dissociating toward to the lowest four N(4S0) + N+(3P), N(2P0) + N+(3P), N(2D0) + N+(3P), and N(4S0) + N+(1D) dissociation limits. Energies of the lowest states, dissociating toward to the N(4S0) + N+(3P) limit, are computed accounting for relativistic corrections. The obtained potential energy curves and the transition dipole moments are employed to compute vibrational energies in these states, vibronic transition dipole moments, and the Einstein coefficients for radiative transitions between the vibronic levels.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142949359","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Generalized Langevin equation for a tagged monomer in a Gaussian semiflexible polymer.","authors":"Xavier Durang, Chan Lim, Jae-Hyung Jeon","doi":"10.1063/5.0229919","DOIUrl":"https://doi.org/10.1063/5.0229919","url":null,"abstract":"<p><p>In this study, we present a comprehensive analysis of the motion of a tagged monomer within a Gaussian semiflexible polymer model. We carefully derived the generalized Langevin equation (GLE) that governs the motion of a tagged central monomer. This derivation involves integrating out all the other degrees of freedom within the polymer chain, thereby yielding an effective description of the viscoelastic motion of the tagged monomer. A critical component of our analysis is the memory kernel that appears in the GLE. By examining this kernel, we characterized the impact of bending rigidity on the non-Markovian diffusion dynamics of the tagged monomer. Furthermore, we calculated the mean-squared displacement of the tagged monomer using the derived GLE. Our theoretical findings were corroborated by the Langevin dynamics simulation and scaling theory. Our results not only show remarkable agreement with previously known results in certain limiting cases but also provide dynamic features over the entire timescale.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142894815","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Variance-reduced random batch Langevin dynamics.","authors":"Zhenli Xu, Yue Zhao, Qi Zhou","doi":"10.1063/5.0246661","DOIUrl":"https://doi.org/10.1063/5.0246661","url":null,"abstract":"<p><p>The random batch method is advantageous in accelerating force calculations in particle simulations, but it poses a challenge of removing the artificial heating effect in application to the Langevin dynamics. We develop an approach to solve this issue by estimating the force variance, resulting in a variance-reduced random batch Langevin dynamics. Theoretical analysis shows the high-order local truncation error of the time step in the numerical discretization scheme, consistent with the fluctuation-dissipation theorem. The numerical results indicate that the method can achieve a significant variance reduction since a smaller batch size provides accurate approximation, demonstrating the attractive feature of the variance-reduced random batch method for Langevin dynamics.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142949437","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A brief introduction to the diffusion Monte Carlo method and the fixed-node approximation.","authors":"Alfonso Annarelli, Dario Alfè, Andrea Zen","doi":"10.1063/5.0232424","DOIUrl":"https://doi.org/10.1063/5.0232424","url":null,"abstract":"<p><p>Quantum Monte Carlo (QMC) methods represent a powerful family of computational techniques for tackling complex quantum many-body problems and performing calculations of stationary state properties. QMC is among the most accurate and powerful approaches to the study of electronic structure, but its application is often hindered by a steep learning curve; hence it is rarely addressed in undergraduate and postgraduate classes. This tutorial is a step toward filling this gap. We offer an introduction to the diffusion Monte Carlo (DMC) method, which aims to solve the imaginary time Schrödinger equation through stochastic sampling of the configuration space. Starting from the theoretical foundations, the discussion leads naturally to the formulation of a step-by-step algorithm. To illustrate how the method works in simplified scenarios, examples such as the harmonic oscillator and the hydrogen atom are provided. The discussion extends to the fixed-node approximation, a crucial approach for addressing the fermionic sign problem in multi-electron systems. In particular, we examine the influence of trial wave function nodal surfaces on the accuracy of DMC energy by evaluating results from a non-interacting two-fermion system. Extending the method to excited states is feasible in principle, but some additional considerations are needed, supported by practical insights. By addressing the fundamental concepts from a hands-on perspective, we hope this tutorial will serve as a valuable guide for researchers and students approaching DMC for the first time.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 24","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142948386","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rikkert J Nap, Paola Carillo Gonzalez, Aria E Coraor, Ranya K A Virk, Juan J de Pablo, Vadim Backman, Igal Szleifer
{"title":"The impact of charge regulation and ionic intranuclear environment on the nucleosome core particle.","authors":"Rikkert J Nap, Paola Carillo Gonzalez, Aria E Coraor, Ranya K A Virk, Juan J de Pablo, Vadim Backman, Igal Szleifer","doi":"10.1063/5.0241529","DOIUrl":"10.1063/5.0241529","url":null,"abstract":"<p><p>We theoretically investigate how the intranuclear environment influences the charge of a nucleosome core particle (NCP)-the fundamental unit of chromatin consisting of DNA wrapped around a core of histone proteins. The molecular-based theory explicitly considers the size, shape, conformation, charge, and chemical state of all molecular species-thereby linking the structural state with the chemical/charged state of the system. We investigate how variations in monovalent and divalent salt concentrations, as well as pH, affect the charge distribution across different regions of an NCP and quantify the impact of charge regulation. The effective charge of an NCP emerges from a delicate and complex balance involving the chemical dissociation equilibrium of the amino acids and the DNA-phosphates, the electrostatic interaction between them, and the translational entropy of the mobile solution ions, i.e., counter ion release and ion condensation. From our results, we note the significant effect of divalent magnesium ions on the charge and electrostatic energy as well as the counterion cloud that surrounds an NCP. As a function of magnesium concentration, charge neutralization, and even charge inversion is predicted-in line with experimental observation of NCPs. The strong Mg-dependence of the nucleosome charge state arises from ion bridges between two DNA-phosphates and one Mg2+ ion. We demonstrate that to describe and predict the charged state of an NCP properly, it is essential to consider molecular details, such as DNA-phosphate ion condensation and the acid-base equilibrium of the amino acids that comprise the core histone proteins.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 23","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142864343","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}