RESUMO
Nonadiabatic excited state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, and charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamics (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Furthermore, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well to those obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.
RESUMO
The crystal nucleation from liquid in most cases is too rare to be accessed within the limited time scales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. Thus, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.
Assuntos
Cristalização , Modelos Químicos , Simulação de Dinâmica MolecularRESUMO
The temperature dependence of the solid-liquid interfacial free energy, γ, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed γ shows a linear dependence on the temperature. The values of γ extrapolated to the melting temperature agree well with previous data obtained by the capillary fluctuation method. Using the temperature dependence of γ, we estimate the nucleation free energy barrier in a wide temperature range from the classical nucleation theory. The obtained data agree very well with the results from the brute-force molecular dynamics simulations.
RESUMO
Molecular dynamics simulations are used to provide strong evidence for barrier-free nucleation events in a heterogeneous solid-solid system. The barrier-free events are characterized by an absence of an incubation time and a growth rate of the emerging phase that is independent of the system size. Furthermore, an analysis of the size and shape of the critical nucleus using the Winterbottom construction indicates that no solution exists for these barrier-free cases. We propose that barrier-free nucleation, which will have a profound effect on phase transformation kinetics, may be a general phenomenon for any polycrystalline material.
RESUMO
To investigate the effect of rosuvastatin on gait and balance disorder progression and elucidate the role of cerebrovascular reactivity (CVR) on this effect. From April 2008 to November 2010, 943 hypertensive patients aged ≥60 years were enrolled from the Shandong area of China. Patients were randomized into rosuvastatin and placebo groups. Gait, balance, CVR, fall and stroke were assessed. During an average 72 months of follow-up, the decreasing trends for step length, step speed, and Berg balance scale scores and the increasing trends for step width and chair rising test were slower in the rosuvastatin group when compared to the placebo group. The hazard ratio of incident balance impairment and falls was 0.542 [95% confidence interval (CI) 0.442-0.663] and 0.532 (95% CI 0.408-0.694), respectively, in the rosuvastatin group compared with placebo group. For CVR progression, the cerebrovascular reserve capacity and breath-holding index were increased and the pulsatility index decreased in the rosuvastatin group, while the cerebrovascular reserve capacity and breath-holding index were decreased, and pulsatility index increased in the placebo group. The changes in gait stability and balance function were independently associated with the changes in the CVR. The odds risks of balance impairment and falls were 2.178 (95% CI: 1.491-3.181) and 3.227 (95% CI: 1.634-6.373), respectively, in the patients with CVR impairment and patients without CVR impairment. Rosuvastatin ameliorated gait and balance disorder progression in older patients with hypertension. This effect might result from the improvement in the CVR. This double-blind clinical trial recruited 943 hypertensive patients aged ≥60 years who were randomly administered rosuvastatin and placebo interventions. The data indicates that rosuvastatin significantly ameliorated the progressions of gait and balance disorders in older hypertensive patients. The cerebrovascular reactivity might play an important mediating role in this amelioration.
Assuntos
Inibidores de Hidroximetilglutaril-CoA Redutases , Hipertensão , Equilíbrio Postural , Rosuvastatina Cálcica , Humanos , Rosuvastatina Cálcica/uso terapêutico , Idoso , Masculino , Feminino , Equilíbrio Postural/efeitos dos fármacos , Hipertensão/tratamento farmacológico , Hipertensão/fisiopatologia , Pessoa de Meia-Idade , Inibidores de Hidroximetilglutaril-CoA Redutases/uso terapêutico , Método Duplo-Cego , Circulação Cerebrovascular/efeitos dos fármacos , Marcha/efeitos dos fármacos , Sulfonamidas/uso terapêutico , Idoso de 80 Anos ou mais , Pirimidinas/uso terapêutico , Acidentes por QuedasRESUMO
We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.
RESUMO
The effects of tripolyphosphate (TPP) on organic pollutant degradation in saline wastewater using Fe0/H2O2 were systematically investigated to elucidate its mechanism and the main reactive oxygen species (ROS). Organic pollutant degradation was dependent on the Fe0 and H2O2 concentration, Fe0/TPP molar ratio, and pH value. The apparent rate constant (kobs) of TPP-Fe0/H2O2 was 5.35 times higher than that of Fe0/H2O2 when orange II (OGII) and NaCl were used as the target pollutant and model salt, respectively. The electron paramagnetic resonance (EPR) and quenching test results showed that â¢OH, O2â¢-, and 1O2 participated in OGII removal, and the dominant ROS were influenced by the Fe0/TPP molar ratio. The presence of TPP accelerates Fe3+/Fe2+ recycling and forms Fe-TPP complexes, which ensures sufficient soluble Fe for H2O2 activation, prevents excessive Fe0 corrosion, and thereby inhibits Fe sludge formation. Additionally, TPP-Fe0/H2O2/NaCl maintained a performance similar to those of other saline systems and effectively removed various organic pollutants. The OGII degradation intermediates were identified using high-performance liquid chromatography-mass spectrometry (HPLC-MS) and density functional theory (DFT), and possible degradation pathways for OGII were proposed. These findings provide a facile and cost-effective Fe-based AOP method for removing organic pollutants from saline wastewater.
RESUMO
Understanding conical intersection (CI) dynamics and subsequent conformational changes is key for exploring and controlling photo-reactions in aromatic molecules. Monitoring of their time-resolved dynamics remains a formidable experimental challenge. In this study, we simulate the photoinduced S3 to S1 non-adiabatic dynamics of cyclooctatetraene (COT), involving multiple CIs with relaxation times in good agreement with experiment. We further investigate the possibility to directly probe the CI passages in COT by off-resonant X-ray Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). We find that these signals sensitively monitor key chemical features during the ultrafast dynamics. First, we distinguish two CIs by using TRUECARS signals with their appearances at different Raman shifts. Second, we demonstrate that TRXD, where X-ray photons scatter off electron densities, can resolve ultrafast changes in the aromaticity of COT. It can further distinguish between planar and non-planar geometries explored during the dynamics, as e.g. two different tetraradical-type CIs. The knowledge gained from these measurements can give unique insight into fundamental chemical properties that dynamically change during non-adiabatic passages.
RESUMO
We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.
RESUMO
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
RESUMO
BACKGROUND AND AIMS: Mitophagy plays a crucial role in mitochondrial homeostasis and is closely associated with endothelial function. However, the mechanism underlying low blood flow shear stress (SS), detrimental cellular stress, regulating endothelial mitophagy is unclear. This study aimed to investigate whether low SS inhibits endothelial mitophagy via caveolin-1 (Cav-1)/miR-7-5p/Sequestosome 1 (SQSTM1) signaling pathway. METHODS: Low SS in vivo modeling was induced using a perivascular SS modifier implanted in the carotid artery of mice. In vitro modeling, low and physiological SS (4 and 15 dyn/cm2, respectively) were exerted on human aortic endothelial cells using a parallel plate chamber system. RESULTS: Compared with physiological SS, low SS significantly inhibited endothelial mitophagy shown by down-regulation of SQSTM1, PINK1, Parkin, and LC 3II expressions. Deficient mitophagy deteriorated mitochondrial dynamics shown by up-regulation of Mfn1 and Fis1 expression and led to decreases in mitochondrial membrane potential. Cav-1 plays a bridging role in the process of low SS inhibiting mitophagy. The up-regulated miR-7-5p expression induced by low SS was suppressed after Cav-1 was silenced. The results of dual-luciferase reporter assays showed that miR-7-5p targeted inhibiting the SQSTM1 gene. Oxidative stress reaction shown by the elevation of reactive oxygen species and O2â- and endothelial dysfunction by the decrease in nitric oxide and the increase in macrophage chemoattractant protein-1 were associated with the low SS inhibiting endothelial mitophagy process. CONCLUSIONS: Cav-1/miR-7-5p/SQSTM1 signaling pathway was the mechanism underlying low SS inhibiting endothelial mitophagy that involves mitochondrial homeostasis impairment and endothelial dysfunction.
Assuntos
MicroRNAs , Mitofagia , Animais , Caveolina 1/genética , Caveolina 1/metabolismo , Células Endoteliais/metabolismo , Humanos , Camundongos , MicroRNAs/genética , MicroRNAs/metabolismo , Mitofagia/genética , Proteínas Quinases/metabolismo , Proteína Sequestossoma-1/genética , Proteína Sequestossoma-1/metabolismo , Transdução de SinaisRESUMO
Background: The heart is an important target organ for the harmful effects of high dietary salt intake. However, the effects and associations of high salt intake on myocardial viability, cardiac function changes, and myocardial remodeling are unclear. Methods: A total of 3,810 participants aged 60 years and older were eligible and enrolled from April 2008 to November 2010 and from August 2019 to November 2019 in the Shandong area of China. Salt intake was estimated using 24-h urine collection consecutively for 7 days. Myocardial strain rates, cardiac function and structure, and serum high-sensitivity C-reactive protein (hsCRP) levels were assessed. Participants were classified into low (n = 643), mild (n = 989), moderate (n = 1,245), and high (n = 933) groups, corresponding to < 6, 6-9, 9-12, and >12 g/day of salt intake, respectively, depending on the salt intake estimation. Results: The global early diastolic strain rate (SRe) and late diastolic strain rate (SRa) in the high group were 1.58 ± 0.26, 1.38 ± 0.24. respectively, and significantly lower compared with the low, mild, and moderate groups (all P < 0.05). The global systolic strain rate (SRs) in the high group was -1.24 ± 0.24, and it was higher than those in the low, mild, and moderate groups (all P < 0.05). Salt intake was independently and positively correlated with global SRs, Tei index, and the parameters of left ventricular structure separately; negatively correlated with global SRe and SRa, left ventricular short axis shortening rate, left ventricular ejection fraction after adjusting for confounders (all P adjusted < 0.001). Hayes process analyses demonstrated that the mediating effects of hsCRP on global SRe, SRa, and SRs; Tei index; and left ventricular remodeling index were -0.013 (95% CI: -0.015 to -0.010), -0.010 (-0.012 to -0.008), 0.008 (0.006-0.010), 0.005 (0.003-0.006), and 0.010 (0.009-0.012), respectively (all P adjusted < 0.001). Conclusion: Our data indicate that excess salt intake is independently associated with the impairment in myocardial viability and cardiac function, as well as myocardial remodeling. Chronic inflammation might play a mediating role in the association between high salt intake and cardiac function damage and myocardial remodeling.
RESUMO
Objective: Small and dense low-density lipoprotein (sdLDL) elevation may be among the most sensitive early biomarkers for nascent cardiovascular disease. This study, therefore, investigated the association between visit-to-visit changes in sdLDL and cerebral small vessel disease (CSVD) progression in older individuals, and the influence of Apolipoprotein E (APOE) genotype on this association. Methods: Between April 2007 and July 2009, 1,143 participants ≥60 years old were recruited from the Shandong region of China, and sdLDL was measured at baseline and at each follow-up visit. White matter hyperintensities (WMHs), lacunes, microbleeds, and enlarged perivascular spaces (EPVSs) were assessed by magnetic resonance imaging. The APOE genotype was determined and participants were stratified as ε4-positive or ε4-negative. Results: During an average follow-up of 86.0 months, 225 participants (19.7%) developed WMH progression, 193 (16.9%) lacune progression, 170 (14.9%) microbleed progression, and 185 (16.2%) EPVS progression. Compared with patients in the first (lowest) tertile of visit-to-visit mean sdLDL, those in the second and third tertiles demonstrated significantly greater risks of WMH progression (53.5 and 105.3% higher), lacune progression (53.3 and 60.8%), microbleed progression (47.2 and 127.6%), and EPVS progression (54.0 and 135.0%) after adjustment for confounders (all adjusted P values for trends <0.001). Compared with patients in the first tertile of visit-to-visit sdLDL SD, those in the second and third tertiles also demonstrated significantly greater risks of WMH progression (49.9% and 143.6%), lacune progression (75.3 and 178.0%), microbleed progression (12.7 and 64.7%), and EPVS progression (41.7 and 114.6%) after adjustment (all P < 0.001). There were significant and positive visit-to-visit mean sdLDL × visit-to-visit sdLDL SD, visit-to-visit mean sdLD×ε4-positive, visit-to-visit sdLDL SD×ε4-positive, and visit-to-visit mean sdLDL×visit-to-visit sdLDL SD×ε4-positive interactions influencing CSVD progression after confounder adjustment (all P < 0.05). Conclusion: Large and variable visit-to-visit changes in sdLDL are independent predictors of aggressive CSVD progression, and this association is strongly influenced by APOE ε4 allele genotype.
RESUMO
In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S2/S3(D2d) â S1(D8h) â S0(D4h) â S0(D2d) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. The obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.
RESUMO
The role of quantum-mechanical coherences in the elementary photophysics of functional optoelectronic molecular materials is currently under active study. Designing and controlling stable coherences arising from concerted vibronic dynamics in organic chromophores is the key for numerous applications. Here, we present fundamental insight into the energy transfer properties of a rigid synthetic heterodimer that has been experimentally engineered to study coherences. Quantum non-adiabatic excited state simulations are used to compute X-ray Raman signals, which are able to sensitively monitor the coherence evolution. Our results verify their vibronic nature, that survives multiple conical intersection passages for several hundred femtoseconds at room temperature. Despite the contributions of highly heterogeneous evolution pathways, the coherences are unambiguously visualized by the experimentally accessible X-ray signals. They offer direct information on the dynamics of electronic and structural degrees of freedom, paving the way for detailed coherence measurements in functional organic materials.
RESUMO
Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum-classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.
RESUMO
The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
RESUMO
Computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. We report an implementation of the fewest-switches surface-hopping algorithm in the NWChem computational chemistry program. The surface-hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited-state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born-Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative nonadiabatic couplings. A numerical approach for the time-derivative nonadiabatic couplings together with an analytical method for calculating nonadiabatic coupling vectors is an efficient combination for surface-hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay timescales, and vibrational coordinates coupled to the excitation dynamics. We also study the photoinduced dynamics in trans-distyrylbenzene. This study provides a baseline for future implementations of higher-level frameworks for simulating nonadiabatic molecular dynamics in NWChem.
RESUMO
We present a versatile new code released for open community use, the nonadiabatic excited state molecular dynamics (NEXMD) package. This software aims to simulate nonadiabatic excited state molecular dynamics using several semiempirical Hamiltonian models. To model such dynamics of a molecular system, the NEXMD uses the fewest-switches surface hopping algorithm, where the probability of transition from one state to another depends on the strength of the derivative nonadiabatic coupling. In addition, there are a number of algorithmic improvements such as empirical decoherence corrections and tracking trivial crossings of electronic states. While the primary intent behind the NEXMD was to simulate nonadiabatic molecular dynamics, the code can also perform geometry optimizations, adiabatic excited state dynamics, and single-point calculations all in vacuum or in a simulated solvent. In this report, first, we lay out the basic theoretical framework underlying the code. Then we present the code's structure and workflow. To demonstrate the functionality of NEXMD in detail, we analyze the photoexcited dynamics of a polyphenylene ethynylene dendrimer (PPE, C30H18) in vacuum and in a continuum solvent. Furthermore, the PPE molecule example serves to highlight the utility of the getexcited.py helper script to form a streamlined workflow. This script, provided with the package, can both set up NEXMD calculations and analyze the results, including, but not limited to, collecting populations, generating an average optical spectrum, and restarting unfinished calculations.