Search

Unidirectional photochemical molecular motors can act as a power source for molecular machines. The motors operate by successive excited state isomerization and ground state helix inversion reactions, attaining unidirectionality from an interplay of steric strain and stereochemistry. Optimizing the yield of the excited state isomerization reaction is an important goal that requires detailed knowledge of excited state dynamics. Here, we investigate the effect of electron withdrawing and donating substituents on excited state structure and ultrafast dynamics in a series of newly synthesized first generation photochemical molecular motors. All substituents red-shift the absorption spectra, while some modify the Stokes shift and render the fluorescence quantum yield solvent polarity dependent. Raman spectra and density functional theory calculations reveal that the stretching mode of the C=C "axle" in the electronic ground state shows a small red-shift when conjugated with electron withdrawing substituents. Ultrafast fluorescence measurements reveal substituent and solvent polarity effects, with the excited state decay being accelerated by both polar solvent environment and electron withdrawing substituents. Excited state structural dynamics are investigated by fluorescence coherence spectroscopy and femtosecond stimulated Raman spectroscopy. The time resolved Raman measurements are shown to provide structural data specifically on the Franck–Condon excited state. The C=C localized modes have a different substituent dependence compared to the ground state, with the unsubstituted motor having the most red-shifted mode. Such measurements provide valuable new insights into pathways to optimize photochemical molecular motor performance, especially if they can be coupled with high-quality quantum molecular dynamics calculations. [ABSTRACT FROM AUTHOR]

Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems. [ABSTRACT FROM AUTHOR]

Quantum dynamics simulations are becoming a standard tool for simulating photo-excited molecular systems involving a manifold of coupled states, known as non-adiabatic dynamics. While these simulations have had many successes in explaining experiments and giving details of non-adiabatic transitions, the question remains as to their predictive power. In this work, we present a set of quantum dynamics simulations on cyclobutanone using both grid-based multi-configuration time-dependent Hartree and direct dynamics variational multi-configuration Gaussian methods. The former used a parameterized vibronic coupling model Hamiltonian, and the latter generated the potential energy surfaces on the fly. The results give a picture of the non-adiabatic behavior of this molecule and were used to calculate the signal from a gas-phase ultrafast electron diffraction (GUED) experiment. Corresponding experimental results will be obtained and presented at a later stage for comparison to test the predictive power of the methods. The results show that over the first 500 fs after photo-excitation to the S2 state, cyclobutanone relaxes quickly to the S1 state, but only a small population relaxes further to the S0 state. No significant transfer of population to the triplet manifold is found. It is predicted that the GUED experiments over this time scale will see signals related mostly to the C–O stretch motion and elongation of the molecular ring along the C–C–O axis. [ABSTRACT FROM AUTHOR]

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born–Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born–Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S 2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S 2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S 1 , where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function— Δ PDF (r) ] for each (interpolated) pathway as well as a time-resolved signal [ Δ PDF (r , t) and Δ I / I (s , t) ] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both Δ PDF (r , t) and Δ I / I (s , t) ] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals. [ABSTRACT FROM AUTHOR]

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [ r L i + < r N a + < r K + ] and higher concentration disturb the hydrogen-bond network of water, resulting in more disordered cationic hydration shell. Our ab initio simulation data for solute concentrations ranging between 1 and 10M can be used to guide future elastic and inelastic neutron-scattering experiments. [ABSTRACT FROM AUTHOR]

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. We show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory. [ABSTRACT FROM AUTHOR]

Microcanonical unimolecular rate theory is applied to Shirhatti and Wodtke's recent supersonic molecular beam experiments examining the activated dissociative chemisorption of HCl on Au(111). A precursor mediated microcanonical trapping (PMMT) model (where the surface vibrates and HCl rotations, vibration, and translation directed along the surface normal are treated as active degrees of freedom) gave dissociative sticking coefficient predictions that are several orders of magnitude higher than experimental values but in good accord with prior quantum and molecular dynamics simulations. Density functional theory (DFT) electronic structure calculations using the Perdew–Burke–Ernzerhof (PBE) functional served to fix the vibrational frequencies of the reactive transition state and the threshold energy for dissociation, E0 = 72.9 kJ/mol. To explore the possibilities of varying threshold energy, coupling to phonons, and dynamics, a three-parameter [E0, s, ɛn] dynamically biased (d-) PMMT model was fit to the experiments. A dynamical bias was introduced using an efficiency, ɛn, of normal translational energy to contribute to the active exchangeable energy capable of promoting reactivity. To achieve the low sticking probabilities observed in experiment, severe normal translational energy dampening (ɛn → 0.26) was imposed, leading to a large vibrational efficacy of ηv = εv/εn = 3.85. The optimal threshold energy for dissociation was E0 = 30.88 kJ/mol, some 40 kJ/mol below the PBE-DFT prediction, and the optimal number of Au surface oscillators was s = 1. The d-PMMT modeling indicates that HCl/Au(111) reactivity can be consistent with electronically adiabatic passage across a relatively low and late transition state that dynamically disfavors normal translational energy. [ABSTRACT FROM AUTHOR]

The molecular dynamics with electronic friction (MDEF) approach can accurately describe nonadiabatic effects at metal surfaces in the weakly nonadiabatic limit. That being said, the MDEF approach treats nuclear motion classically such that the nuclear quantum effects are completely missing in the approach. To address this limitation, we combine Electronic Friction with Ring Polymer Molecular Dynamics (EF-RPMD). In particular, we apply the averaged electronic friction from the metal surface to the centroid mode of the ring polymer. We benchmark our approach against quantum dynamics to show that EF-RPMD can accurately capture zero-point energy as well as transition dynamics. In addition, we show that EF-RPMD can correctly predict the electronic transfer rate near metal surfaces in the tunneling limit as well as the barrier crossing limit. We expect that our approach will be very useful to study nonadiabatic dynamics near metal surfaces when nuclear quantum effects become essential. [ABSTRACT FROM AUTHOR]

We present a novel theoretical formulation for performing quantum dynamics in terms of moments within the single-particle description. By expressing the quantum dynamics in terms of increasing orders of moments, instead of single-particle wave functions as generally done in time-dependent density functional theory, we describe an approach for reducing the high computational cost of simulating the quantum dynamics. The equation of motion is given for the moments by deriving analytical expressions for the first-order and second-order time derivatives of the moments, and a numerical scheme is developed for performing quantum dynamics by expanding the moments in the Taylor series as done in classical molecular dynamics simulations. We propose a few numerical approaches using this theoretical formalism on a simple one-dimensional model system, for which an analytically exact solution can be derived. The application of the approaches to an anharmonic system is also discussed to illustrate their generality. We also discuss the use of an artificial neural network model to circumvent the numerical evaluation of the second-order time derivatives of the moments, as analogously done in the context of classical molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B12H122−·nCH3CN (n = 1–4). The electron binding energies of B12H122−·nCH3CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B12H122− via a threefold dihydrogen bond (DHB) B3–H3 ⁝⁝⁝ H3C–CN unit, in which three adjacent nucleophilic H atoms in B12H122− interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B12H122− solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B12H122− and CH3CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B12H122−·CH3CN cluster. The binding motif in B12H122−·CH3CN is largely retained for the whole halogenated series B12X122−·CH3CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions. [ABSTRACT FROM AUTHOR]

The manipulation of the ultrafast quantum dynamics of a molecular system can be achieved through the application of tailored light fields. This has been done in many ways in the past. In our present investigation, we show that it is possible to exert specific control over the nonadiabatic dynamics of a generic model system describing ultrafast charge-transfer within a condensed dissipative environment by using frequency-chirped pulses. By adjusting the external photoexcitation conditions, such as the chirp parameter, we show that the final population of the excitonic and charge-transfer states can be significantly altered, thereby influencing the elementary steps controlling the transfer process. In addition, we introduce an excitation scheme based on double-pump time-resolved fluorescence spectroscopy using chirped-pulse excitations. Here, our findings reveal that chirped excitations enhance the vibrational system dynamics as evidenced by the simulated spectra, where a substantial signal intensity dependence on the chirp is observed. Our simulations show that chirped pulses are a promising tool for steering the dynamics of the charge-transfer process toward a desired target outcome. [ABSTRACT FROM AUTHOR]

Disordered molecular systems, such as amorphous catalysts, organic thin films, electrolyte solutions, and water, are at the cutting edge of computational exploration at present. Traditional simulations of such systems at length scales relevant to experiments in practice require a compromise between model accuracy and quality of sampling. To address this problem, we have developed an approach based on generative machine learning called the Morphological Autoregressive Protocol (MAP), which provides computational access to mesoscale disordered molecular configurations at linear cost at generation for materials in which structural correlations decay sufficiently rapidly. The algorithm is implemented using an augmented PixelCNN deep learning architecture that, as we previously demonstrated, produces excellent results in 2 dimensions (2D) for mono-elemental molecular systems. Here, we extend our implementation to multi-elemental 3D and demonstrate performance using water as our test system in two scenarios: (1) liquid water and (2) samples conditioned on the presence of pre-selected motifs. We trained the model on small-scale samples of liquid water produced using path-integral molecular dynamics simulations, including nuclear quantum effects under ambient conditions. MAP-generated water configurations are shown to accurately reproduce the properties of the training set and to produce stable trajectories when used as initial conditions in quantum dynamics simulations. We expect our approach to perform equally well on other disordered molecular systems in which structural correlations decay sufficiently fast while offering unique advantages in situations when the disorder is quenched rather than equilibrated. [ABSTRACT FROM AUTHOR]

C4F7N/CO2 gas mixtures have attracted extensive attention because of their excellent insulating properties and environmental friendliness. High electrical and thermal stability is an important indicator for evaluating their performance, but there have been few molecular dynamics studies of their decomposition mechanisms. In this study, using ReaxFF molecular dynamics simulations and quantum chemistry theory, the decomposition mechanism of a C4F7N/CO2 gas mixture and the effect of the O2 content on the decomposition of the mixture were simulated on the microscopic level. It was found that there are three main decomposition pathways of C4F7N molecules, of which the generation of C3F4N⋅ and CF3⋅ free radicals is the most likely to occur. COF2 is the main oxygen-containing product of the C4F7N/CO2 gas mixture, and its generation is significantly affected by the simulation time and temperature. COF2 can be regarded as the characteristic decomposition product of the C4F7N/CO2 gas mixture. The addition of O2 slightly promotes the decomposition of C4F7N, whereas the maximum decomposition rate of CO2 decreases by 0.3% and 1% after the addition of 2% and 8% O2, respectively. Relevant results of this research can provide a theoretical basis and guidance for research into the performance of C4F7N/CO2 gas mixtures and practical engineering applications of these mixtures in the future. [ABSTRACT FROM AUTHOR]

This work presents systematic comparisons between classical molecular dynamics (cMD) and quantum dynamics (QD) simulations of 15-dimensional and 75-dimensional models in their description of H atom scattering from graphene. We use an experimentally validated full-dimensional neural network potential energy surface of a hydrogen atom interacting with a large cell of graphene containing 24 carbon atoms. For quantum dynamics simulations, we apply Monte Carlo canonical polyadic decomposition to transform the original potential energy surface (PES) into a sum of products form and use the multi-layer multi-configuration time-dependent Hartree method to simulate the quantum scattering of a hydrogen or deuterium atom with an initial kinetic energy of 1.96 or 0.96 eV and an incident angle of 0°, i.e., perpendicular to the graphene surface. The cMD and QD initial conditions have been carefully chosen in order to be as close as possible. Our results show little differences between cMD and QD simulations when the incident energy of the H atom is equal to 1.96 eV. However, a large difference in sticking probability is observed when the incident energy of the H atom is equal to 0.96 eV, indicating the predominance of quantum effects. To the best of our knowledge, our work provides the first benchmark of quantum against classical simulations for a system of this size with a realistic PES. Additionally, new projectors are implemented in the Heidelberg multi-configuration time-dependent Hartree package for the calculation of the atom scattering energy transfer distribution as a function of outgoing angles. [ABSTRACT FROM AUTHOR]

We develop a deep learning-based algorithm, called DeepForce, to link ab initio physics with the continuum theory to predict concentration profiles of confined water. We show that the deep-learned forces can be used to predict the structural properties of water confined in a nanochannel with quantum scale accuracy by solving the continuum theory given by Nernst–Planck equation. The DeepForce model has an excellent predictive performance with a relative error less than 7.6% not only for confined water in small channel systems (L < 6 nm) but also for confined water in large channel systems (L = 20 nm) which are computationally inaccessible through the high accuracy ab initio molecular dynamics simulations. Finally, we note that classical Molecular dynamics simulations can be inaccurate in capturing the interfacial physics of water in confinement (L < 4.0 nm) when quantum scale physics are neglected. [ABSTRACT FROM AUTHOR]

The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures. [ABSTRACT FROM AUTHOR]

To enhance the performance of polyimides, novel polyimides (P1, P2, P3) were synthesized through the reaction of isophthalic dianhydride with a series of furan-based diamines. The structures and properties of these polyimides (P1, P2, P3) were analyzed through molecular dynamics and quantum chemical simulation methods, including glass transition temperature, dielectric properties, free volume, and cohesive energy density. Compared with traditional polyimides, the results showed that P1, P2, and P3 exhibited a significant increase in glass transition temperature and cohesive energy density, while their dielectric properties decreased compared to traditional PI (3.4–3.6), all falling within the range of (2.5-3.0). The free volume also decreased. However, there was no significant impact on the materials' radial distribution functions at the atomic level. [ABSTRACT FROM AUTHOR]

Among azobenzene derivatives, azothiophenes represent a relatively recent family of compounds that exhibit similar characteristics as dyes and photoreactive systems. Their technological applications are extensive thanks to the additional design flexibility conferred by the heteroaromatic ring. In this study, we present a comprehensive investigation of the structural and electronic properties of novel dyes derived from 3-thiophenamine, utilizing a multilevel approach. We thoroughly examined the potential energy surfaces of the E and Z isomers for three molecules, each bearing different substituents on the phenyl ring at the para position relative to the diazo group. This exploration was conducted through quantum chemistry calculations at various levels of theory, employing a continuum solvent model. Subsequently, we incorporated an explicit solvent (a dimethyl sulfoxide–water mixture) to simulate the most stable isomers using classical molecular dynamics, delivering a clear picture of the local solvation structure and intermolecular interactions. Finally, a hybrid quantum mechanics/molecular mechanics (QM/MM) approach was employed to accurately describe the evolution of the solutes' properties within their environment, accounting for finite temperature effects. [ABSTRACT FROM AUTHOR]

Molecular imprinting is a promising approach for developing polymeric materials as artificial receptors. However, only a few types of molecularly imprinted polymers (MIPs) are commercially available, and most research on MIPS is still in the experimental phase. The significant limitation has been a challenge for screening imprinting systems, particularly for weak functional target molecules. Herein, a combined method of quantum mechanics (QM) computations and molecular dynamics (MD) simulations was employed to screen an appropriate 2,4-dichlorophenoxyacetic acid (2,4-D) imprinting system. QM calculations were performed using the Gaussian 09 software. MD simulations were conducted using the Gromacs2018.8 software suite. The QM computation results were consistent with those of the MD simulations. In the MD simulations, a realistic model of the 'actual' pre-polymerisation mixture was obtained by introducing numerous components in the simulations to thoroughly investigate all non-covalent interactions during imprinting. This study systematically examined MIP systems using computer simulations and established a theoretical prediction model for the affinity and selectivity of MIPs. The combined method of QM computations and MD simulations provides a robust foundation for the rational design of MIPs. [ABSTRACT FROM AUTHOR]

The results of molecular dynamics simulation of the crystal structure of the currently most studied CQDs, based on cadmium selenide nanoparticles, are presented, and issues concerning the structure of the ligand shell of these nanoparticles, primarily made of trioctylphosphine (TOP) and trioctylphosphine oxide (TOPO) molecules, which act as a medium and as precursors in the chemical reaction that results in the formation of CQDs, are discussed in detail. [ABSTRACT FROM AUTHOR]

• A good accumulation and accommodation of interface is towards acetylacetone (AcAc). • AcAc isomerization is favorable on acidic interface as "water-bridge" is destroyed. • Carbonyl/hydroxyl O-atom of AcAc displays an energetical preference to hydration. • Unfeasible hydration on acidic interface leads to the accumulation of product. Acetylacetone (AcAc) is a typical class of β-diketones with broad industrial applications due to the property of the keto-enol isomers, but its isomerization and chemical reactions at the air-droplet interface are still unclear. Hence, using combined molecular dynamics and quantum chemistry methods, the heterogeneous chemistry of AcAc at the air-droplet interface was investigated, including the attraction of AcAc isomers by the droplets, the distribution of isomers at the air-droplet interface, and the hydration reactions of isomers at the air-droplet interface. The results reveal that the preferential orientation of two AcAc isomers (keto - and enol -AcAc) to accumulate and accommodate at the acidic air-droplet interface. The isomerization of two AcAc isomers at the acidic air-droplet interface is more favorable than that at the neutral air-droplet interface because the "water bridge" structure is destroyed by H 3 O+, especially for the isomerization from keto -AcAc to enol -AcAc. At the acidic air-droplet interface, the carbonyl or hydroxyl O-atoms of two AcAc isomers display an energetical preference to hydration. Keto -diol is the dominant products to accumulate at the air-droplet interface, and excessive keto -diol can enter the droplet interior to engage in the oligomerization. The photooxidation reaction of AcAc will increase the acidity of the air-droplet interface, which indirectly facilitate the uptake and formation of more keto -diol. Our results provide an insight into the heterogeneous chemistry of β-diketones and their influence on the environment. [Display omitted] [ABSTRACT FROM AUTHOR]

Piperine (a naturally occurring alkaloid), is a biologically active natural product belonging to the alkaloid series of compounds. In this study, the spectral characterization of piperine isolated from black pepper has been performed using FT-IR, 1H-NMR, and 13C-NMR. To explain spectral features, we have calculated the vibrational frequencies of the optimized structure of piperine using the B3LYP/6-311 + G(d,p) level of theory and NMR spectra using the GIAO method. All vibrational modes of the title molecule have been assigned based on potential energy distribution. The calculated scaled wavenumbers and NMR chemical shifts show good agreement with the experimental FT-IR and NMR data, respectively. To assess the bioactivity of piperine, we have first performed molecular docking using the ITK receptor and subsequently, the molecular dynamics simulation of the resulting complex for 100 ns. The results suggest the potential of piperine for possible anti-inflammatory action. [ABSTRACT FROM AUTHOR]

We report a theoretical investigation of a new carbon nanobelt called Hexa-Penta-belt with molecular formula C $_{42}$ 42 H $_{18}$ 18 . Hexa-Penta-belt is a variation of the synthesised Methylene-bridged[6]cycloparaphenylene (M-b[6]CPP) with molecular formula C $_{42}$ 42 H $_{24}$ 24 . Hexa-Penta-belt consists of carbon pentagon and hexagon rings and was inspired by the Penta-belt (theoretically derived from the penta-graphene sheet), the [12]cyclophenacene, and the M-b[6]CPP synthesised by Itami and Segawa [1]. Using Density Functional Theory (DFT), we calculated Hexa-Penta-belt electronic, optical, and thermodynamic properties and their vibrational frequencies. Hexa-Penta-belt is a structurally and thermally stable molecule, as Quantum Molecular Dynamics (QMD) showed. The results exhibited that the Hexa-Penta-belt has excellent electronic stability comparable to C $_{60}$ 60 , and it behaves as a semiconductor with a gap estimated at 0.89 eV. It absorbs in the ultraviolet and the visible region (corresponding to the violet, green, and red colours) presenting an absorption more consistent in the visible region than M-b[6]CPP, which could indicate application, for example, as an optical nanosensor. We verified through the carbon atom binding energy that Hexa-Penta-belt is an intermediate molecule between Penta-belt and [12]cyclophenacene. [ABSTRACT FROM AUTHOR]

In this work, a composite material between chitosan and silver oxide nanoparticles (Ag2NPs) was developed. The universal force field option in the Materials Studio® software was applied to the modelled Ag2NPs within a sphere of variable radius, being 11 å the most stable. In addition, each of these models was taken to their state of minimum energy to verify their stability, and then, a simulation was carried out using a canonical ensemble at 298 K. The X-ray diffraction pattern of these models was simulated and is according to those reported for Ag2O. Molecular and quantum dynamics simulations for the models of the chitosan chains in interaction with Ag2ONPs were made. Each proposed model was based on the Schrödinger equation, considering the quantum fields for the interaction of the nanoparticle in the chitosan substrate. It was observed that the fields were deformed, demonstrating that the presence of the nanoparticle on the surface produces deformation in the polymeric matrix, which suggests a phenomenon of surface adsorption that was confirmed with the RDF where a chemical interaction between the polymer and the silver oxide nanoparticles was shown. This result supports that is possible to obtain the chitosan–Ag2ONPs experimentally. [ABSTRACT FROM AUTHOR]

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. The proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich–Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer–Miller–Stock–Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately. [ABSTRACT FROM AUTHOR]

In this work, we prepare a simulation framework for a high‐accuracy numerical study of electron–ion temperature relaxation in nonideal (strongly coupled) plasmas. The existing relaxation rate theories require either parameter selection or some pre‐knowledge of the electron–ion correlation functions and effective interaction potentials. This makes non‐equilibrium classical and quantum molecular dynamics simulations a crucial stage in the study of energy transfer rates. We begin by revisiting the classical molecular dynamics simulations of a system of equally charged particles with different masses on a neutralizing background. We accurately simulate this simple ab‐initio (parameterless) system with controlled precision in terms of number of particles, mass ratio, and energy convergence. The predictions for the equally charged system are compared to the previous simulations and theories, which are reproduced with higher accuracy. We also perform a series of classical molecular dynamics simulations of the system of oppositely charged particles with the corrected Kelbg potential based on the quantum statistical approach. We analyze the differences and similarities between the same‐charge and opposite‐charge systems. Some remarks about the forthcoming application of quantum simulations with the help of WPMD or WPMD‐DFT methods are given. [ABSTRACT FROM AUTHOR]

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter. [ABSTRACT FROM AUTHOR]

Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs. Through neutron scattering experiments coupled with machine learning, the authors uncover the strong role of nuclear quantum effects in the dynamics of ammonia in both its solid and technologically relevant liquid phase. [ABSTRACT FROM AUTHOR]

We consider the problem of numerically computing the quantum dynamics of an electron in twisted bilayer graphene. The challenge is that atomic-scale models of the dynamics are aperiodic for generic twist angles because of the incommensurability of the layers. The Bistritzer--MacDonald PDE model, which is periodic with respect to the bilayer's moir\'e pattern, has recently been shown to rigorously describe these dynamics in a parameter regime. In this work, we first prove that the dynamics of the tight-binding model of incommensurate twisted bilayer graphene can be approximated by computations on finite domains. The main ingredient of this proof is a speed of propagation estimate proved using Combes--Thomas estimates. We then provide extensive numerical computations, which clarify the range of validity of the Bistritzer--MacDonald model. [ABSTRACT FROM AUTHOR]

In most cases, the ultrafast dynamics of resonantly excited molecules are considered and almost always computed in the molecular frame, while experiments are carried out in the laboratory frame. Here, we provide a formalism in terms of a lab frame density matrix, which connects quantum dynamics in the molecular frame to those in the laboratory frame, providing a transparent link between computation and measurement. The formalism reveals that in any such experiment, the molecular frame dynamics vary for molecules in different orientations and that certain coherences, which are potentially experimentally accessible, are rejected by the orientation-averaged reduced vibronic density matrix. Instead, molecular angular distribution moments are introduced as a more accurate representation of experimentally accessible information. Furthermore, the formalism provides a clear definition of a molecular frame quantum tomography and specifies the requirements to perform such a measurement enabling the experimental imaging of molecular frame vibronic dynamics. Successful completion of such a measurement fully characterizes the molecular frame quantum dynamics for a molecule at any orientation in the laboratory frame. [ABSTRACT FROM AUTHOR]

We study theoretically the quantum dynamics and spectroscopy of rovibrational polaritons formed in a model system composed of a single rovibrating diatomic molecule, which interacts with two degenerate, orthogonally polarized modes of an optical Fabry–Pérot cavity. We employ an effective rovibrational Pauli–Fierz Hamiltonian in length gauge representation and identify three-state vibro-polaritonic conical intersections (VPCIs) between singly excited vibro-polaritonic states in a two-dimensional angular coordinate branching space. The lower and upper vibrational polaritons are of mixed light–matter hybrid character, whereas the intermediate state is purely photonic in nature. The VPCIs provide effective population transfer channels between singly excited vibrational polaritons, which manifest in rich interference patterns in rotational densities. Spectroscopically, three bright singly excited states are identified when an external infrared laser field couples to both a molecular and a cavity mode. The non-trivial VPCI topology manifests as pronounced multi-peak progression in the spectral region of the upper vibrational polariton, which is traced back to the emergence of rovibro-polaritonic light–matter hybrid states. Experimentally, ubiquitous spontaneous emission from cavity modes induces a dissipative reduction of intensity and peak broadening, which mainly influences the purely photonic intermediate state peak as well as the rovibro-polaritonic progression. [ABSTRACT FROM AUTHOR]

Methylation of amines inside an introverted resorcinarene‐based deep methyl ester cavitand is investigated by means of molecular dynamics simulations and quantum chemical calculations. Experimentally, the cavitand has been shown to bind a number of amines and accelerate the methylation reaction by more than four orders of magnitude for some of them. Eight different amines are considered in the present study, and the geometries and energies of their binding to the cavitand are first characterized and analyzed. Next, the methyl transfer reactions are investigated and the calculated barriers are found to be in generally good agreement with experimental results. In particular, the experimentally‐observed rate acceleration in the cavitand as compared to the solution reaction is well reproduced by the calculations. The origins of this rate acceleration are analyzed by computational modifications made to the structure of the cavitand, and the role of the solvent is discussed. [ABSTRACT FROM AUTHOR]

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid–liquid phase transition and liquid–liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., T g LDL (P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., T g HDL (P) increases upon compression. (iii) NQE shift both T g LDL (P) and T g HDL (P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA–HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA–HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA–HDA transformation. [ABSTRACT FROM AUTHOR]

Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum–classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum–classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems. [ABSTRACT FROM AUTHOR]

Spin–orbit (SO) coupling can have significant effects on the quantum dynamics of molecular systems, but it is still difficult to account for accurately. One promising way to do this is to devise a diabatic SO model combined with the molecular potential energy. Few such models have been developed utilizing spatial and time-reversal symmetry. These models are tedious to derive and are specific for the molecular symmetry and included spin states. Here, we present a relatively simple approach to construct such models for various spin states with S ≠ 1 2 from a basic one-electron SO case with S = 1 2 . The multi-electron fine structure states are expressed in terms of Slater determinants of single-electron spin functions (spinors). The properties of all single-electron matrix elements over the SO operator are derived and expressed as Taylor expansions in terms of symmetry-adapted nuclear coordinates. The SO matrix elements for the multi-electron case are then obtained from these single-electron matrix elements using the Slater–Condon rules. This yields the full SO matrix and symmetry properties of the multi-electron matrix elements in a straightforward way. The matrix elements are expressed as symmetry-adapted polynomials up to arbitrary order. This approach is demonstrated first for an abstract model of two electrons in a set of p orbitals in a C3v symmetric environment and then applied to set up a diabatic model for the photodissociation of methyl iodide (CH3I). The high accuracy of this new approach is demonstrated in comparison to an available analytic SO model for CH3I. [ABSTRACT FROM AUTHOR]

Quasicentroid molecular dynamics (QCMD) is a path-integral method for approximating nuclear quantum effects in dynamics simulations, which has given promising results for gas- and condensed-phase water. In this work, by simulating the infrared spectrum of gas-phase ammonia, we test the feasibility of extending QCMD beyond water. Overall, QCMD works as well for ammonia as for water, reducing or eliminating blue shifts from the classical spectrum without introducing the artificial red shifts or broadening associated with other imaginary-time path-integral methods. However, QCMD gives only a modest improvement over the classical spectrum for the position of the symmetric bend mode, which is highly anharmonic (since it correlates with the inversion pathway). We expect QCMD to have similar problems with large-amplitude degrees of freedom in other molecules but otherwise to work as well as for water. [ABSTRACT FROM AUTHOR]

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su–Schrieffer–Heeger (Holstein–Peierls) model Hamiltonian. [ABSTRACT FROM AUTHOR]

The method of correlation femtoscopy makes it possible to estimate the parameters of the particle-emitting region (radius of emission region, , and correlation strength, ). Measurement of femtoscopic radius dependence on transverse momentum of particle pairs, , is an important tool for studying the dynamics of the emission process. This work is devoted to the study of momentum correlations of identical pions produced in Au Au collisions at GeV using the UrQMD (Ultrarelativistic Quantum Molecular Dynamics) model. Three-dimensional femtoscopic analysis was performed as a function of , rapidity and collision centrality. Physical implications will be discussed. [ABSTRACT FROM AUTHOR]

SF6 is widely used in gas switchgear but has a strong greenhouse effect. The development of an environment-friendly switchgear that can replace SF6 is a current research hotspot. As a SF6 alternative technology, vacuum arc-extinguishing chamber plus clean air insulation has shown high application prospects. The clean air inside the switchgear decomposes under the effect of high temperature and discharge; however, only a few studies focused on its decomposition mechanism. In this work, the decomposition mechanism of clean air and the effect of temperature on the decomposition are simulated at the atomic level based on ReaxFF (Reactive force field) molecular dynamics and quantum chemistry theory. Results showed that the decomposition of clean air mainly generates NO, NO2, and N2O. NO is the main product at high temperatures and thus can be the characteristic decomposition product of clean air. The clean air has good self-recovery characteristics, and its decomposition can be substantially promoted by increasing the temperature. The decomposition rates of N2 and O2 under 3000 K can reach 7.00% and 8.00%, respectively, which are twice and four times those under 2000 K. These results can provide theoretical basis and engineering guidance for the development of environment-friendly switchgear with vacuum arc-extinguishing chamber plus clean air insulation. [ABSTRACT FROM AUTHOR]

We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques. [ABSTRACT FROM AUTHOR]

We present a new data-driven potential energy function (PEF) describing chloride–water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride–water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride–water and water–water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments. [ABSTRACT FROM AUTHOR]

Expanding the set of stable, accurate, and scalable methods for simulating molecular quantum dynamics is important for accelerating the computational exploration of molecular processes. In this paper, we adapt the signed particles Monte Carlo algorithm for solving the transient Wigner equation to scenarios of chemical interest. This approach was used in the past to study electronic processes in semi-conductors, but to the best of our knowledge, it had never been applied to molecular modeling. We present the algorithm and demonstrate its excellent performance on harmonic and double well potentials for electronic and nuclear systems. We explore the stability of the algorithm, discuss the choice of hyper-parameters, and cautiously speculate that it may be used in quantum molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann–Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems. [ABSTRACT FROM AUTHOR]

In this paper, we report on the performance of various quantum molecular dynamics simulation methods in describing the photo-induced nonadiabatic dynamics underlying the isomerization process of the retinal chromophore in rhodopsin. We focus on purely quantum vibronic wavepacket techniques and on various trajectory-based schemes, discussing their capability of accurately capture the isomerization process using a two-dimensional two-state model system coupled to an environment of secondary harmonic modes. Numerical results of various algorithms and time-independent grid schemes for the purely quantum approaches are presented, which also serve as benchmark for the trajectory-based calculations. Independent-trajectory and coupled-trajectory methods are compared as well, devoting particular attention to the scaling of their computational cost when increasing the number of degrees of freedom. [ABSTRACT FROM AUTHOR]

We study the dynamical fluctuation ν dyn in the particle yield ratio K / π produced in Au + Au interactions at the laboratory energy E lab = 30 A GeV using the ultra-relativistic quantum molecular dynamics (UrQMD) and a multi-phase transport (AMPT) models. Additionally, we investigate some single particle properties of K ± and π ± and bulk properties of the matter created in the interactions under consideration. The single particle properties like the centrality dependence of normalized particle yield and mean transverse momentum measured from the models are found to differ from the corresponding experimental values. The dynamical fluctuation in K / π , denoted by ν K π , dyn , are positive valued, while the fluctuation measured for the identified particle pairs, like K + / π + , K + / π - , etc. are negative valued at all centralities. Correlations between the oppositely charged particle pairs are found to be larger than that between the same charged particle pairs. Also we study the effect of hadronic rescattering in the dynamical fluctuations by using both the UrQMD and AMPT models. [ABSTRACT FROM AUTHOR]

A novel CaCO3 scale inhibitor IK-CD) was successfully prepared by amine modification of kaolinite and loading carboxymethyl cyclodextrin onto the surface of kaolinite. The scale inhibition effect of IK-CD was deeply analyzed by Ca2e titration experiment scale inhibition experiment molecular dynamics simulation and quantum chemistry calculation. The results show that IK-CD can decrease the precipitation rate and delay the occurrence of CaCO3 precipitation. Its scale inhibition efficiency is up to 90 which is 9.3 higher than that of carboxymethyl cyclodextrin. Molecular dynamics simulations show that the main adsorption surface of amine modified kaolinite is (001) surface which can attract calcium carbonate crystal nucleus. The nucleation process of Ca2e and carbonate attached onto the active site of calcium carbonate crystal nucleus caused lattice distortion inhibited the aggregation process of calcium carbonate and achieved scale inhibition. [ABSTRACT FROM AUTHOR]

To study the structure and dynamics of nitrogen-containing compounds, NMR parameters with directly involved nitrogen can provide valuable structure information. However, this information can only be obtained using 15N-enriched compounds due to low natural abundance of 15N and extremely short relaxation time of 14N. In the synthesis of benzamides from their 15N-ammonium salts, 15N-enriched benzamides are often used as intermediates. In the present work, we studied the dynamic structure of benzamide, which is controlled by two independent factors: hindered internal rotation of the NH2 group around the C(O)–N bond and of the amide group as a whole relative to the benzene ring. Deeper knowledge of the mechanism and parameters of these processes in amides is important for meaningful interpretation and prediction of the biological activity of aromatic amides in living systems and the strength and conformation of their supramolecular complexes with lanthanide and actinide ions. A double enriched [2H5,15N]benzamide was synthesized to avoid undesirable superposition of the strong aromatic multiplet on the amide signals in the 1H NMR spectra. The 1H NMR spectrum of this compound contained only strong signals of amide protons, which allowed accurate determination of the quantitative characteristics of the studied dynamic processes. The obtained experimental data are in good agreement with the results of quantum molecular dynamics simulation. [ABSTRACT FROM AUTHOR]