Search

The atomistic description of the folding process of a structureless chain of amino acids to a functioning protein is still considered to be a challenge for computational biophysics. An in-depth understanding of protein-folding can be achieved through atomistic dynamics accounting for the solvent effects, combined with the theoretical description of conformational changes in shorter polypeptide chains. The present paper studies the folding transitions in short polypeptide chains consisting of 11 alanine residues in explicit solvent employing all-atom molecular dynamics. The multiple observed folding ↔ unfolding events of the peptide are interpreted as a dynamic process and rationalised through the analysis of the potential energy surface of the system. It is demonstrated that alanine polypeptide folding dynamics is governed by the backbone dihedral angles and involves the formation of spontaneously folded α -helices which emerge and live for ca. 1–20 picoseconds. The helical content within the polypeptide at different temperatures was quantified through a statistical mechanics approach, which showed a reasonable agreement with the results of molecular dynamics simulations and experiments performed for alanine-rich peptides. [ABSTRACT FROM AUTHOR]

Density-functional-theory-based chemical reactivity indicators are formulated for degenerate and near-degenerate ground states. For degenerate states, the functional derivatives of the energy with respect to the external potential do not exist, and must be replaced by the weaker concept of functional variation. The resultant reactivity indicators depend on the specific perturbation. Because it is sometimes impractical to compute reactivity indicators for a specific perturbation, we consider two special cases: point-charge perturbations and Dirac delta function perturbations. The Dirac delta function perturbations provide upper bounds on the chemical reactivity. Reactivity indicators using the common used 'average of degenerate states approximation' for degenerate states provide a lower bound on the chemical reactivity. Unfortunately, this lower bound is often extremely weak. Approximate formulas for the reactivity indicators within the frontier-molecular-orbital approximation and special cases (two or three degenerate spatial orbitals) are presented in the supplementary material. One remarkable feature that arises in the frontier molecular orbital approximation, and presumably also in the exact theory, is that removing electrons sometimes causes the electron density to increase at the location of a negative (attractive) Dirac delta function perturbation. That is, the energetic response to a reduction in the external potential can increase even when the number of electrons decreases. [ABSTRACT FROM AUTHOR]

Radon gas was allowed to accumulate in its radium source and then injected into a 36 m3 test room, resulting in an initial radon concentration of 15 kBq m–3. Filter papers were used to collect the short-lived radon progeny and thus to measure the Potential Alpha Energy Concentration (PAEC) in-situ in the year 1984 at different times and conditions according to the experimental design. The radon progeny collected on the filter papers were studied as a function of aerosol particle concentration ranging from 102–105 particles cm–3 in three different experiments. The highest aerosol particle concentration was generated by indoor cigarette smoking. Those filters were stored after the experiment, and were used after 16 years to study the activity of the radon long-lived alpha emitter progeny, 210Po (T1/2=138 days). This isotope is separated from the short-lived progeny by 210Pb beta emitter with 22.3 years half-life. After 16 years’ storage of these filters, each filter paper was sandwiched and wrapped between two CR-39 nuclear track detectors, to put the detectors in contact with the surfaces of different filters, for 337 days. Correlation between the PAEC measured using filter papers in the year 1984 and the activity of long-lived alpha emitter 210Po on the same filter papers measured in the year 2000 were studied. The results of the 210Po activity showed a very good correlation of 0.92 with the PAEC 16 years ago. The results also depict that the PAEC and 210Po activity in indoor air increased with the increase of aerosol particle concentration, which shows the attachment of short-lived radon progeny with the aerosol particles. The experiment proves that indoor cigarette smoking is a major source of aerosol particles carrying radon progeny and, thus, indoor cigarette smoking is an additional source of internal radiation hazard to the occupants whether smoker or non-smoker. [Copyright &y& Elsevier]

The Journal of Chemical Physics has released a special issue focused on software for atomistic machine learning. This issue aims to address the lack of journals dedicated to publishing scientific software papers. The collection of papers in this issue provides insight into the tools and goals of software implementations in the field of atomistic machine learning. The articles cover a range of topics, including machine-learning interatomic potentials, sampling, dataset repositories, workflows, and auxiliary tooling and analysis. The article concludes by emphasizing the importance of software implementations in the field and encourages further submissions on relevant topics. [Extracted from the article]

The potential energy surfaces of the reactions of Cl4MCH2(M = Cr, Mo, W, Ru, Re) with ethylene, models of potential chain-carrying catalysts and olefins respectively in the transition metal-catalyzed olefin metathesis reaction, have been explored using density functional theory at the B3LYP/LACVP* level of theory. In Cl4MCH2(M = Cr, Ru), the carbenoid complexes Cl3MCH2–Cl have been found to be more stable than the corresponding carbene Cl4MCH2complexes whereas in Cl4MCH2(M = Mo, W, Re) the carbene complexes are more stable than the carbenoid complexes. The carbenoid complexes have been found not to favor the formation of metallacyclobutanes, a key step in the olefin metathesis reaction according to the Herrison–Chauvin mechanism, indicating that the active species for the metathesis reaction is a carbene complex and not a carbenoid complex. Therefore, even though the formation of metallacyclobutanes through formal [2+2] cycloaddition has been found to be a low-barrier process with each of the metal carbene complexes investigated, metathesis is likely to occur only in Cl4MCH2(M = Mo, W, Re) but not in Cl4MCH2(M = Cr, Ru) since the reaction surface in the latter complexes is likely to be populated by the carbenoid complex rather than the carbene complex. In Cl4WCH2the kinetic and thermodynamic preference of the productive [2+2] pathway leading to the metallacyclobutane over the [3+2] pathway is unambiguous whereas in Cl4MoCH2the [3+2] pathway is likely to be competitive with the [2+2] pathway. The metallacyclobutane formation in W and Re has been found to have lower barriers than in Mo, suggesting that the W and Re complexes may have a greater metathesis activity than the Mo complex. [ABSTRACT FROM AUTHOR]

This work follows a companion article, which will be referred to as Paper I [Campeggio et al., J. Chem. Phys. 158, 244104 (2023)] in which a quantum-stochastic Liouville equation for the description of the quantum–classical dynamics of a molecule in a dissipative bath has been formulated in curvilinear internal coordinates. In such an approach, the coordinates of the system are separated into three subsets: the quantum coordinates, the classical relevant nuclear degrees of freedom, and the classical irrelevant (bath) coordinates. The equation has been derived in natural internal coordinates, which are bond lengths, bond angles, and dihedral angles. The resulting equation needs to be parameterized. In particular, one needs to compute the potential energy surfaces, the friction tensor, and the rate constants for the nonradiative jumps among the quantum states. While standard methods exist for the calculation of energy and dissipative properties, an efficient evaluation of the transition rates needs to be developed. In this paper, an approximated treatment is introduced, which leads to a simple explicit formula with a single adjustable parameter. Such an approximated expression is compared with the exact calculation of transition rates obtained via molecular dynamics simulations. To make such a comparison possible, a simple sandbox system has been used, with two quantum states and a single internal coordinate (together with its conjugate momentum). Results show that the adjustable parameter, which is an effective decoherence time, can be parameterized from the effective relaxation times of the autocorrelation functions of the conjugated momenta of the relevant nuclear coordinates. [ABSTRACT FROM AUTHOR]

Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances. [ABSTRACT FROM AUTHOR]

The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen–Bunker–Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems. [ABSTRACT FROM AUTHOR]

This paper proposes a combination approach termed combined approach for dynamic available transmission capacity (CDATC) for the quick and accurate computation of ATC. In this regard, the suggested model makes use of the static enlarged downhill and Newton–Raphson-Seidel algorithm. In this paper, the boundary surface of potential energy and the maximum point of potential energy are also estimated in a manner to ascertain transient stability. A CDATC calculation was developed by combining the static method, transient stability, and online state estimate. When compared to existing DATC approaches, the simulation results on a large-scale real-world system confirm the suitability of the suggested method with faster speed and tighter accuracy. Simulation results on the East Iran network, IEEE 118 bus, and Iowa State 145 bus, all of which have a high penetration of wind farms-showed that the developed CDATC is robust at all operating conditions and can be used online with fewer computational resources. The unique application of state estimation in DATC computation, which enabled the study to produce a more realistic answer even in the face of ambiguity, is an additional value. [ABSTRACT FROM AUTHOR]

The rotationally resolved infrared (IR) spectrum of the He– H 3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν2 mode of the H 3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He– H 3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He– H 3 + complex is extremely floppy, with nearly unhindered internal rotation of the H 3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He– H 3 + , higher levels of excitation, including the complex stretching mode, need further attention. [ABSTRACT FROM AUTHOR]

We propose a new collocation multi-configuration time-dependent Hartree (MCTDH) method. It reduces point-set error by using more points than basis functions. Collocation makes it possible to use MCTDH with a general potential energy surface without computing any integrals. The collocation points are associated with a basis larger than the basis used to represent wavefunctions. Both bases are obtained from a direct product basis built from single-particle functions by imposing a pruning condition. The collocation points are those on a sparse grid. Heretofore, collocation MCTDH calculations with more points than basis functions have only been possible if both the collocation grid and the basis set are direct products. In this paper, we exploit a new pseudo-inverse to use both more points than basis functions and a pruned basis and grid. We demonstrate that, for a calculation of the lowest 50 vibrational states (energy levels and wavefunctions) of CH2NH, errors can be reduced by two orders of magnitude by increasing the number of points, without increasing the basis size. This is true also when unrefined time-independent points are used. [ABSTRACT FROM AUTHOR]

We present the computational methodology, which for the first time allows rigorous twelve-dimensional (12D) quantum calculations of the coupled intramolecular and intermolecular vibrational states of hydrogen-bonded trimers of flexible diatomic molecules. Its starting point is the approach that we introduced recently for fully coupled 9D quantum calculations of the intermolecular vibrational states of noncovalently bound trimers comprised of diatomics treated as rigid. In this paper, it is extended to include the intramolecular stretching coordinates of the three diatomic monomers. The cornerstone of our 12D methodology is the partitioning of the full vibrational Hamiltonian of the trimer into two reduced-dimension Hamiltonians, one in 9D for the intermolecular degrees of freedom (DOFs) and another in 3D for the intramolecular vibrations of the trimer, and a remainder term. These two Hamiltonians are diagonalized separately, and a fraction of their respective 9D and 3D eigenstates is included in the 12D product contracted basis for both the intra- and intermolecular DOFs, in which the matrix of the full 12D vibrational Hamiltonian of the trimer is diagonalized. This methodology is implemented in the 12D quantum calculations of the coupled intra- and intermolecular vibrational states of the hydrogen-bonded HF trimer on an ab initio calculated potential energy surface (PES). The calculations encompass the one- and two-quanta intramolecular HF-stretch excited vibrational states of the trimer and low-energy intermolecular vibrational states in the intramolecular vibrational manifolds of interest. They reveal several interesting manifestations of significant coupling between the intra- and intermolecular vibrational modes of (HF)3. The 12D calculations also show that the frequencies of the v = 1, 2 HF stretching states of the HF trimer are strongly redshifted in comparison to those of the isolated HF monomer. Moreover, the magnitudes of these trimer redshifts are much larger than that of the redshift for the stretching fundamental of the donor-HF moiety in (HF)2, most likely due to the cooperative hydrogen bonding in (HF)3. The agreement between the 12D results and the limited spectroscopic data for the HF trimer, while satisfactory, leaves room for improvement and points to the need for a more accurate PES. [ABSTRACT FROM AUTHOR]

Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin–orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129–4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions. [ABSTRACT FROM AUTHOR]

The impact of the chalcogen atomic electronegativity (O, S, and Se atoms) of new organic molecules on excited-state dynamical reactions is self-evident. Inspired by this kind of distinguished photochemical characteristic, in this work, we performed a computational investigation of chalcogen-substituted 3,6-bis(4,5-dihydroxyoxazo-2-yl)benzene-1,2-diol (BDYBD) derivatives (i.e., BDYBD-O, BDYBD-S, and BDYBD-Se). In this paper, we pay close attention to characteristic BDYBD derivatives that contain intramolecular double hydrogen bonds (O1–H2···N3 and O4–H5···N6). The main goal of this study was to explore how changes in atomic electronegativity affect the way hydrogen bonds interact and how excited molecules affect transfer protons. We go into further detail in the main text of the paper. By fixing our attention to geometrical variations and infrared (IR) vibrational spectra between the S0 and S1 states, exploring hydrogen bonding behaviors using the core-valence bifurcation (CVB) index, and simulating hydrogen bonding energy (EHB) via the atom in molecule (AIM) method, we clarified the photo-induced strengthened dual hydrogen bonding interactions that facilitate the excited-state dual-proton transfer (ESDPT) behavior of BDYBD derivatives. The reorganization of charge stemming from photoexcitation further verifies the tendencies of ESDPT reactions. We relied on constructing potential energy surfaces (PESs) by adopting a restrictive optimization approach, and herein, we finally clarify the gradual ESDPT mechanism of BDYBD derivatives. Particularly, we confirm that the variation in chalcogen atomic electronegativity has a regulatory effect on the ESDPT behavior of BDYBD derivatives; that is, the lower the atomic electronegativity, the more favorable it is for the gradual ESDPT reaction. [ABSTRACT FROM AUTHOR]

This paper concerns the representation of angular momentum operators in the Born–Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born–Oppenheimer parameter κ = (m/M)1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born–Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that block-diagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born–Oppenheimer theory. Our analysis accommodates large-amplitude motions and is not dependent on small-amplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space. [ABSTRACT FROM AUTHOR]

Copyright of Momento: Revista de Física is the property of Universidad Nacional de Colombia, Departamento de Fisica and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The geometrical re parameters of trans-azobenzene (E-AB) free molecule were refined by gas electron diffraction (GED) method using available experimental data obtained previously by S. Konaka and coworkers. Structural analysis was carried out by various techniques. First of all, these included the widely used molecular orbital constrained gas electron diffraction method and regularization method. The results of the refinements using different models were also compared—a semirigid model, three variants of one-dimensional dynamic models, and a two-dimensional pseudoconformer model. Several descriptions have been used due to the fact that E-AB has a shallow potential energy surface along the rotation coordinates of phenyl groups. Despite this, it turned out that the semirigid model is suitable for use for E-AB and allows good agreement with experimental data to be achieved. According to the results of GED structural analysis, coupled with the results of DLPNO-CCSD(T0) calculations, E-AB has a planar structure. Based only on GED data, it is impossible to unambiguously determine the rotational angle of the phenyl group due to the facts that (i) with rotation over a wide range of angles, the bonded distances in the molecule change insignificantly and (ii) potential function in a structural analysis within a dynamic model is not determined with the necessary accuracy. This work also examines the sensitivity of the GED method to structural changes caused by trans-cis isomerization. The paper also analyzes the applicability of different variants of density functional theory (DFT) calculations in GED structural analysis using E-AB as an example. There are not enough similar methodological works in the literature. This experimental and methodological information is especially important and relevant for planning and implementing GED experiments and corresponding processing of the results for azobenzene derivatives, in which the conformer and isomeric diversity are even more complicated due to the presence of different substituents. [ABSTRACT FROM AUTHOR]

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Qy transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies. [ABSTRACT FROM AUTHOR]

Some diarylethene molecular switches have a low quantum yield for cycloreversion when excited by a single photon, but react more efficiently following sequential two-photon excitation. The increase in reaction efficiency depends on both the relative time delay and the wavelength of the second photon. This paper examines the wavelength-dependent mechanism for sequential excitation using excited-state resonance Raman spectroscopy to probe the ultrafast (sub-30 fs) dynamics on the upper electronic state following secondary excitation. The approach uses femtosecond stimulated Raman scattering (FSRS) to measure the time-gated, excited-state resonance Raman spectrum in resonance with two different excited-state absorption bands. The relative intensities of the Raman bands reveal the initial dynamics in the higher-lying states, Sn, by providing information on the relative gradients of the potential energy surfaces that are accessed via secondary excitation. The excited-state resonance Raman spectra reveal specific modes that become enhanced depending on the Raman excitation wavelength, 750 or 400 nm. Many of the modes that become enhanced in the 750 nm FSRS spectrum are assigned as vibrational motions localized on the central cyclohexadiene ring. Many of the modes that become enhanced in the 400 nm FSRS spectrum are assigned as motions along the conjugated backbone and peripheral phenyl rings. These observations are consistent with earlier measurements that showed higher efficiency following secondary excitation into the lower excited-state absorption band and illustrate a powerful new way to probe the ultrafast dynamics of higher-lying excited states immediately following sequential two-photon excitation. [ABSTRACT FROM AUTHOR]

The structure–property relationship is always the focus when organic light-emitting materials are used mainly in the aggregated state in some fields, such as optoelectronic devices, sensors, and bioimaging. To clarify a clear structure–property relationship between the aggregate structure and the luminescent property, the intermolecular dimer, which represents the smallest aggregate in a number, is selected as an ideal research model. However, how to achieve an intermolecular dimer by a rational molecular design is rarely found in previous papers in this field. Here, we design and synthesize three pyrene derivatives, aiming at achieving the pyrene-based dimer in a solid state, by bridging pyrene and benzo[d]thiazole (BZT) in a different position on a phenylene. The pyrene and BZT groups are connected with each other in the ortho-, meta-, and para-positions of the phenylene to construct three molecules Py-oBZT, Py-mBZT, and Py-pBZT. Experimental and theoretical investigations demonstrate that only Py-oBZT shows pyrene-based dimer stacking with a discrete feature due to the large steric hindrance provided by the BZT group. Compared to individual pyrenes in the crystal, the pyrene-based dimer in the Py-oBZT crystal exhibits a different molecular orientation, which is verified by the presence of multiple potential wells of the pyrene dimer using a theoretical simulation of the potential energy surface. This work not only examines the effect of steric hindrance on the formation of the pyrene dimer packing, but also provides an understanding of the multiplicity of the pyrene dimer structures. [ABSTRACT FROM AUTHOR]

Nuclear dynamics in the coupled electronic states of mono-, di-(ortho and meta), and pentafluorobenzene molecules is investigated here. Attempts are specifically made to understand the complexity and broadening of the recorded gas phase electronic absorption spectra of these molecules. Justification is also provided for the low quantum yield of fluorescence emission with increasing number of fluorine substitutions. The nuclear dynamics is simulated from first principles both by time-independent and time-dependent quantum mechanical methods. Potential energy surfaces of the low-lying excited electronic states of these molecules constructed in Paper I [Mondal and Mahapatra, J. Chem. Phys. 133, 084304 (2010)] are employed for the purpose. Theoretical results presented in this paper are compared with the available experimental data and the agreement between the two is found to be excellent. While structured electronic absorption bands are observed for the S1 state of mono- and difluorobenzene molecules, the same for the pentafluorobenzene is broad and structureless. Occurrence of S1-S2 conical intersections in pentafluorobenzene leads to a nonradiative internal conversion of the S1 state in ∼165 fs and contributes to the broadening of the S1←S0 absorption band and a biexponential decay of its fluorescence emission. [ABSTRACT FROM AUTHOR]

Summary: In the present paper, the static approximate Newton‐Raphson‐Sidel (ANRS) and the developed down‐hill (DDH) algorithm methods were used. To take advantage of both methods, ANRS‐DDH and the revised holomorphic embedding load flow (RHELF) were combined. An approximation of the potential energy boundary surface (APEBS) was also employed in combination with determine transient stability and correct the previous methods. A combination of static method and transient stability led to a combined approach for dynamic available transfer capability (C‐DATC) computation. It is shown that a small correction for these four methods makes the proposed method faster or more accurate. The average root mean square error of the proposed method is half the other methods and the average relative speed of calculations in the proposed methods is 10% higher than the other method. These were confirmed in the networks of Iowa State (145 bus), modified IEEE 118‐bus, and the west of Iran (848 bus) with a large number of wind farms. Finally, the present paper demonstrated that the suggested approach is ideal not only for the operation mode but also for the planning mode and the existence of wind farms. [ABSTRACT FROM AUTHOR]

The dissociating states of the triplet-excited CO-He complex are studied by means of scattering calculations on ab initio diabatic potential energy surfaces produced in the preceding paper (Paper I). With the aid of an effective transition dipole function and the bound states of the CO-He complex in the ground singlet state we obtain the photoabsorption cross section as a function of the excitation energy and generate the full UV spectrum of the singlet-triplet transition. It was found that the dominant contributions to the spectrum, in the energy range from -5 to +10 cm[SUP-1] relative to the band origin at 48 473.201 cm[SUP-1], originate from resonances that correspond to higher spin-orbit levels of the excited CO(a[SUP3]П) - He complex with approximate quantum number |Ω| = 1. Rapid predissociation, with the triplet CO fragment decaying into its lower spin-orbit levels with Ω = 0, limits the lifetime of these excited levels to, typically, 10-700 ps. We also predict the rotational and spin-orbit state distribution of the triplet CO fragment and the maximum deflection angle of the photodissociation products in a molecular beam experiment. [ABSTRACT FROM AUTHOR]

This paper continues the investigation of the exponentially repulsive EXP pair-potential system of Paper I [A. K. Bacher et al., J. Chem. Phys. 149, 114501 (2018)] with a focus on isomorphs in the low- temperature gas and liquid phases. As expected from the EXP system's strong virial potential-energy correlations, the reduced-unit structure and dynamics are isomorph invariant to a good approximation. Three methods for generating isomorphs are compared: the small-step method that is exact in the limit of small density changes and two versions of the direct-isomorph-check method that allows for much larger density changes. Results from the latter two approximate methods are compared to those of the small-step method for each of the three isomorphs generated by 230 one percent density changes, covering one decade of density variation. Both approximate methods work well [ABSTRACT FROM AUTHOR]

It is uncontested that a machine learning scheme cannot correctly reproduce physical properties that vary on a manifold in configuration space if the fingerprint, used as an input for the machine learning scheme, is constant on this manifold. In our original paper, we found several manifolds whose fingerprint variation is sufficiently small to prevent machine learning based on a standard training scheme even if some 40 configurations on the manifold are included in the training set. Response to "Comment on 'Manifolds of quasi-constant SOAP and ACSF fingerprints and the resulting failure to machine learn four-body interactions'" [J. Chem. Phys. [Extracted from the article]

In this paper, we have developed a unitary variant of a double exponential coupled cluster theory, which is capable of handling molecular strong correlation with arbitrary electronic complexity. With the Hartree–Fock determinant taken as the reference, we introduce a sequential product of parameterized unitary Ansätze. While the first unitary, containing the excitation operators, acts directly on the reference determinant, the second unitary, containing a set of rank-two, vacuum-annihilating scattering operators, has nontrivial action only on certain entangled states. We demonstrate the theoretical bottleneck of such an implementation in a classical computer, whereas the same is implemented in the hybrid quantum–classical variational quantum eigensolver framework with a reasonably shallow quantum circuit without any additional approximation. We have further introduced a number of variants of the proposed Ansatz with different degrees of sophistication by judiciously approximating the scattering operators. With a number of applications on strongly correlated molecules, we have shown that all our schemes can perform uniformly well throughout the molecular potential energy surface without significant additional implementation cost over the conventional unitary coupled cluster approach with single and double excitations. [ABSTRACT FROM AUTHOR]

This paper presents the experimental and theoretical analyses of ArICl(IP,vIP,nIP) states' population and decay at energies lower than the ArICl(E,vE = 0,nE) dissociation limit (IP = E0+, D′2, β1), vIP = 0, 1, and nIP are the quantum numbers of the van der Waals (vdW) modes. We have measured the excitation spectra of the ArICl(E,vE = 0,1,nE → X,vX,nX) and ArICl(β,0,nβ → A and/or D ′ , v D ′ , n D ′ → A ′ luminescence as well as luminescence spectra themselves. To construct potential energy surfaces (PESs) for valence (A1, A′2) and ion-pair (E, β, and D′) electronic states of the complex, we utilized the intermolecular diatomic-in-molecule perturbation theory first order method. The experimental and calculated spectroscopic characteristics of the T-shaped ArICl valence and E, β states agree well. The ArICl(D′) state PES has no vdW levels in the T-shaped configuration, and collinear ArICl(D′) binding energy is larger than that of the T-shaped ArICl(β) state. We calculated vibrational state energies and the ArICl(IP → valence states) luminescence excitation spectra, as well as luminescence spectra themselves, by using the Heidelberg MCTDH code. The comparison of the experimental and calculated excitation spectra shows that the latter describe their principal features. The bound–bound ArICl(E,0,nE → X and β,0,nβ → A) parts of experimental luminescence spectra are described adequately by the calculated spectra, whereas bound-free ArICl(E,0,nE → X, D′, 0, nD′ → A′) parts are not described since the bound-free transitions occur in repulsive parts of the ArICl(X, A ′ PESs, which we cannot describe accurately. [ABSTRACT FROM AUTHOR]

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules. We show that HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant errors. The ability of the HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems, which require a larger active space than possible with a complete active space self-consistent field. Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT). We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional. We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet–triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface. [ABSTRACT FROM AUTHOR]

The mechanism of vibrational relaxation and dissociation in the O2–O system at elevated temperatures is investigated by means of molecular dynamics. The most recent O3 potential energy surfaces (PESs), obtained from the first principles quantum mechanical calculations [Varga et al., J. Chem. Phys. 147, 154312 (2017)], are used to derive a complete set of state-specific rate coefficients of vibrational energy transfer and dissociation. Unlike most of the previous efforts that utilize only the lowest and supposedly most reactive 11A′ O3 PES [A. Varandas and A. Pais, Mol. Phys. 65, 843 (1988)], this paper demonstrates the necessity to account for a complete ensemble of all excited O3 PESs that correlate with O2(X) and O(3P) when high-temperature kinetics is of interest. At the same time, it is found that the Varandas 11A′ O3 PES adequately describes vibrational energy transfer and dissociating dynamics when compared to the most recent 11A′ O3 PES by Varga et al. [J. Chem Phys. 147, 154312 (2017)]. The differences between this new dataset and previous rate coefficients are quantified by the master equation model. [ABSTRACT FROM AUTHOR]

Context: With the union of machine learning (ML) and quantum chemistry, amid the debate between machine-learned functionals and human-designed functionals in density functional theory (DFT), this paper aims to demonstrate the generation of potential energy surfaces using computations with machine-learned density functional approximation (ML-DFA). A recent research trend is the application of ML in quantum sciences in the design of density functionals such as DeepMind's Deep Learning model (DeepMind21, DM21). Though science reported the state-of-the-art performance of DM21, the opportunity to utilize DeepMind's pretrained DM21 neural networks in computations in quantum chemistry has not yet been tapped. So far in the literature, the Deep Learning density functionals (DM21) have not been applied to generate potential energy surfaces. While the superior accuracy of DM21 has been reported, there is still a scarcity of publications that apply DM21 in calculations in the field. In this context, for the first time in literature, neural density functionals inferring 2D potential energy surfaces (ML-DFA-PES) based on machine-learned DFA-based computational method is contributed in this paper. This paper reports the ML-DFA-generated PES for C4H8, H2O, H2, and H2+ by employing a pretrained DM21m TensorFlow model with cc-pVDZ basis set. In addition, we also analyze the long-range behavior of DM21 based PES to investigate the ability to describe a system at long ranges. Furthermore, we compare PES diagrams from DM21 with popular DFT functionals (b3lyp/ PW6B95) and CCSD(T). Methods: In this method, 2D potential energy surfaces are obtained using a method that relies upon the neural network's ability to accurately learn the mapping between 3D electron density and exchange–correlation potential. By inserting Deep Learning inference in DFT with a pretrained neural network, self-consistent field (SCF) energy at different geometries along the coordinates of interest is computed, and then, potential energy surfaces are plotted. In this method, first, the electron density is computed mathematically, and this computed 3D electron density is used as a ML feature vector to predict the exchange correlation potential as a ML inference computed by a forward pass of pre-trained DM21 TensorFlow computational graph, followed by the computation of self-consistent field energy at multiple geometries, and then, SCF energies at different bond lengths/angles are plotted as 2D PES. We implement this in a python source code using frameworks such as PySCF and DM21. This paper contributes this implementation in open source. The source code and DM21-DFA-based PES are contributed at https://sites.google.com/view/MLfunctionals-DeepMind-PES. [ABSTRACT FROM AUTHOR]

In internal combustion piston engines, the formation of an oil film is completely different from that seen in industrial machines. The molecular adhesion force at the interface between the surface coating of engine parts and the lubricating oil determines the load-carrying capacity and the ability to form a lubricated film. The geometry of the lubricating wedge between the surfaces of the piston rings and the cylinder wall is created by the thickness of the oil film and the height of the ring's coverage with lubricating oil. This condition is affected by many of the parameters that characterize the engine's operation and the physical and chemical parameters of the coatings used for the cooperating pairs. For lubricant particles that reach energies that are higher than the potential energy barrier regarding adhesive attraction at the interface, slippage occurs. Therefore, the value of the contact angle of the liquid on the surface of the coating depends on the value of the intermolecular force of attraction. According to the current author, there is a strong relationship between the contact angle and the lubrication effect. The paper shows that the surface potential energy barrier is a function of the contact angle and contact angle hysteresis (CAH). The innovation of the current work consists in examining the contact angle and CAH under the conditions of thin layers of lubricating oil, in cooperation with hydrophilic and hydrophobic coatings. The thickness of the lubricant film was measured under various speed and load conditions, using optical interferometry. The study shows that CAH is a better interfacial parameter for correlation with the effect of hydrodynamic lubrication. This paper presents the mathematical relationships relating to a piston engine, various coatings, and lubricants. [ABSTRACT FROM AUTHOR]

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction. [ABSTRACT FROM AUTHOR]

Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born–Oppenheimer approximation. In this paper, we analyze the electron–photon Hamiltonian within the cavity Born–Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. In particular, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron–photon correlation and the accuracy of CBOA. [ABSTRACT FROM AUTHOR]

Accurate modeling of anionic abundances in the interstellar and circumstellar media requires calculations of collisional data with the most abundant species that are usually He atoms and H2 molecules. In this paper, we focus on smaller cyclic molecular anion, c-C3H−, an astrophysical candidate, following the detection of larger CnH− carbon chains. From a new three-dimensional potential energy surface, the rotational (de-)excitation of the c-C3H−(X1A1) anion by collision with He is investigated. The surface is obtained in the supermolecular approach at the CCSD(T)-F12/aug-cc-pVTZ level of theory. Fully quantum close-coupling calculations of inelastic integral cross sections are performed on a grid of collisional energies large enough to ensure the convergence of the state-to-state rate coefficients for the 34 first rotational levels up to j K a , K c = 77,0 of c-C3H− and temperatures ranging from 5 to 100 K. For this collisional system, rate coefficients exhibit a strong dominance in favor of 21,2 → l1,1 downward transition. This transition was previously used for the detection of the cyclic parent c-C3H. The c-C3H−–He rate coefficients (∼10−11 cm3 s−1) are of the same order of magnitude as those of the detected anions CnH− (as C2H−, C4H−, and C6H−) in collision with He and one order of magnitude smaller than those with H2. The critical densities of H2 were also estimated, and a discussion on the validity of the local thermodynamic equilibrium conditions is carried out. This work represents the contribution to understanding and modeling abundances and chemistry of hydrocarbon radicals, CnH, in astrophysical media. [ABSTRACT FROM AUTHOR]

Recent work has demonstrated the promise of using machine-learned surrogates, in particular, Gaussian process (GP) surrogates, in reducing the number of electronic structure calculations (ESCs) needed to perform surrogate model based (SMB) geometry optimization. In this paper, we study geometry meta-optimization with GP surrogates where a SMB optimizer additionally learns from its past "experience" performing geometry optimization. To validate this idea, we start with the simplest setting where a geometry meta-optimizer learns from previous optimizations of the same molecule with different initial-guess geometries. We give empirical evidence that geometry meta-optimization with GP surrogates is effective and requires less tuning compared to SMB optimization with GP surrogates on the ANI-1 dataset of off-equilibrium initial structures of small organic molecules. Unlike SMB optimization where a surrogate should be immediately useful for optimizing a given geometry, a surrogate in geometry meta-optimization has more flexibility because it can distribute its ESC savings across a set of geometries. Indeed, we find that GP surrogates that preserve rotational invariance provide increased marginal ESC savings across geometries. As a more stringent test, we also apply geometry meta-optimization to conformational search on a hand-constructed dataset of hydrocarbons and alcohols. We observe that while SMB optimization and geometry meta-optimization do save on ESCs, they also tend to miss higher energy conformers compared to standard geometry optimization. We believe that further research into characterizing the divergence between GP surrogates and potential energy surfaces is critical not only for advancing geometry meta-optimization but also for exploring the potential of machine-learned surrogates in geometry optimization in general. [ABSTRACT FROM AUTHOR]

Hybrid methods that combine molecular dynamics methods capable of analyzing dynamics with Monte Carlo (MC) methods that can efficiently treat thermodynamically stable states are valuable for understanding complex chemical processes in which an equilibrium state is reached through many elementary processes. The hybrid MC (HMC) method is one such promising method; however, it often fails to sample configurations properly from the canonical multimodal distribution due to the rugged potential energy surfaces. In this paper, we extend the HMC method to overcome this difficulty. The new method, which is termed potential scaling HMC (PS-HMC), makes use of an artificially modulated trajectory to propose a new configuration. The trajectory is generated from Hamilton's equations, but the potential energy surface is scaled to be gradually flattened and then recovered to the original surface, which facilitates barrier-crossing processes. We apply the PS-HMC method to three kinds of molecular processes: the thermal motion of argon particles, butane isomerization, and an atom transfer chemical reaction. These applications demonstrate that the PS-HMC method is capable of correctly constructing the canonical ensemble with a multimodal distribution. The sampling efficiency and accepted trajectories are examined to clarify the features of the PS-HMC method. Despite the potential scaling, many reactive atom transfer trajectories (elementary processes) pass through the vicinity of the minimum energy path. Furthermore, we demonstrate that the method can properly imitate the relaxation process owing to the inherent configurational continuity. By comparing the PS-HMC method with other relevant methods, we can conclude that the new method is a unique approach for studying both the dynamic and thermodynamic aspects of chemical processes. [ABSTRACT FROM AUTHOR]

The results of an extensive ab initio study of the cyanobutadiyne cation, initially motivated by threshold-photoelectron spectroscopy experiments [see the study by Gans et al., J. Chem. Phys. 150, 244304 (2019)], are reported in the present paper. Calculations at the internally contracted multireference configuration interaction level of theory have been performed to derive the rovibronic properties of the seven lowest electronic states of HC5N+. Equilibrium geometries, rotational constants, vibrational frequencies, electric dipole moments, and spin-orbit constants have been calculated and compared with experimental data when available. Adiabatic and vertical ionization energies from the neutral ground state as well as transition energies within the cation electronic manifold are predicted, using the convergence to the complete basis set limit. The accurate description of the complex energy landscape up to 32 000 cm−1 above the ionization potential allows us to perform Franck-Condon simulations of the photoionization spectrum to the X+2Π, A+2Π, B+2Σ+, and C+2Π states and allows us to simulate the A+2Π → X+2Π emission spectrum. The vibronic perturbations occurring on the excited potential energy surfaces are revealed and discussed, in particular, for the 3 2Π surface, which presents a double-well topography. [ABSTRACT FROM AUTHOR]

In the companion paper [T. S. Ingebrigtsen, S. Toxvaerd, O. J. Heilmann, T. B. Schro\der, and J. C. Dyre, 'NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface,' J. Chem. Phys. (in press)] an algorithm was developed for tracing out a geodesic curve on the constant-potential-energy hypersurface. Here, simulations of NVU dynamics are compared to results for four other dynamics, both deterministic and stochastic. First, NVU dynamics is compared to the standard energy-conserving Newtonian NVE dynamics by simulations of the Kob-Andersen binary Lennard-Jones liquid, its WCA version (i.e., with cut-off's at the pair potential minima), and the Lennard-Jones Gaussian liquid. We find identical results for all quantities probed: radial distribution functions, incoherent intermediate scattering functions, and mean-square displacement as function of time. Arguments are presented for the equivalence of NVU and NVE dynamics in the thermodynamic limit; in particular, to leading order in 1/N these two dynamics give identical time-autocorrelation functions. In the final part of the paper, NVU dynamics is compared to Monte Carlo dynamics, to a diffusive dynamics of small-step random walks on the constant-potential-energy hypersurface, and to Nosé-Hoover NVT dynamics. If time is scaled for the two stochastic dynamics to make single-particle diffusion constants identical to that of NVE dynamics, the simulations show that all five dynamics are equivalent at low temperatures except at short times. [ABSTRACT FROM AUTHOR]

In a previous paper [J. L. McAfee and B. Poirier, J. Chem. Phys. 130, 064701 (2009)], using spin-polarized density functional theory (DFT), the authors reported a binding energy of 0.755 eV, for a single hydrogen atom adsorbed on a pristine (unrelaxed) (5,5) single-walled carbon nanotube (SWNT) substrate. A full three-dimensional (3D) potential energy surface (PES) for the SWNT-H system was also developed, and used in a quantum dynamics calculation to compute all rovibrational bound states, and associated equatorial and longitudinal adsorbate migration rates. A highly pronounced preference for the latter migration pathway at ambient temperatures was observed. In this work, we extend the aforementioned study to include multiple H-atom adsorbates. Extensive DFT calculations are performed, in order to ascertain the most relevant dynamical pathways. For two adsorbates, the SWNT-H-H system is found to exhibit highly site-specific binding, as well as long-range correlation and pronounced binding energy enhancement. The latter effect is even more pronounced in the full-hydrogenation limit, increasing the per-adsorbate binding energy to 2.6 eV. To study migration dynamics, a single-hole model is developed, for which the binding energy drops to 2.11 eV. A global 3D PES is developed for the hole migration model, using 40 radial × 18 cylindrical ab initio geometries, fit to a Fourier basis with radially dependent expansion coefficients (rms error 4.9 meV). As compared with the single-adsorbate case, the hole migration PES does not exhibit separate chemisorption and physisorption wells. The barrier to longitudinal migration is also found to be much lower. Quantum dynamics calculations for all rovibrational states are then performed (using a mixed spectral basis/phase-space optimized discrete variable representation), and used to compute longitudinal migration rates. Ramifications for the use of SWNTs as potential hydrogen storage materials are discussed. [ABSTRACT FROM AUTHOR]

The variational approach detailed in the previous paper (Paper I) for the treatment of the Renner–Teller effect in linear tetra-atomic molecules including all degrees of freedom and couplings between angular momenta is applied for HCCH+. The accurate six-dimensional potential energy surfaces of the X 2Πu electronic state, presented in Paper I is incorporated in the variational treatment in order to obtain all rovibronic levels including the spin-orbit coupling for 1/2≤J≤7/2 and up to 2600 cm-1 above the global zero point energy. The “pure” stretching levels are calculated up to 11 000 cm-1 from the stretching zero point energy. The calculated rovibronic energies are compared with previous theoretical and experimental data. The mean agreement with the zero kinetic energy photoelectron measurements of Tang et al. [J. Chem. Phys. 125, 133201 (2006)] is of 16.7 cm-1. The Renner–Teller parameters have been determined at νtrans=690.0 cm-1, εtrans=0.30, νcis=715.0 cm-1, and εcis=-0.063. A detailed analysis of the rovibronic Hund’s cases is presented and the rotational structures of some vibronic bands recorded by Yang and Mo [J. Phys. Chem. A 110, 11001 (2006)] are given. [ABSTRACT FROM AUTHOR]

Crystal structure prediction has been a subject of topical interest but remains a substantial challenge especially for complex structures as it deals with the global minimization of the extremely rugged high-dimensional potential energy surface. In this paper, a symmetry-orientated divide-and-conquer scheme was proposed to construct a symmetry tree graph, where the entire search space is decomposed into a finite number of symmetry dependent subspaces. An artificial intelligence-based symmetry selection strategy was subsequently devised to select the low-lying subspaces with high symmetries for global exploration and in-depth exploitation. Our approach can significantly simplify the problem of crystal structure prediction by avoiding exploration of the most complex P1 subspace on the entire search space and has the advantage of preserving the crystal symmetry during structure evolution, making it well suitable for predicting the complex crystal structures. The effectiveness of the method has been validated by successful prediction of the candidate structures of binary Lennard-Jones mixtures and the high-pressure phase of ice, containing more than 100 atoms in the simulation cell. The work therefore opens up an opportunity toward achieving the long-sought goal of crystal structure prediction of complex systems. [ABSTRACT FROM AUTHOR]

Cyclodextrins (CDs) constitute a class of cyclic oligosaccharides that are well recognized and largely applied in the drug delivery field, thanks to their biocompatibility, low cost, and the possibility to be derivatized in order to tune and optimize the complexation/release of the specific drug. The conformational flexibility of these systems is one of their key properties and requires a cost-effective methodology to be studied by combining the accuracy of results with the possibility of exploring a large set of conformations. In the present paper, we have explored the conformational potential energy surface of the monomers and dimers of α-, β-, and γ-cyclodextrins (i.e., 6, 7, and 8 monomeric units, respectively) by means of fast but accurate semiempirical methods, which are then refined by state-of-the-art DFT functionals. Moreover, the crystal structure is considered for a more suitable comparison with the IR spectrum experimentally recorded. Calculations are carried out in the gas phase and in water environments, applying both implicit and explicit treatments. We show that the conformation of the studied molecules changes from the gas phase to the water, even if treated implicitly, thus modifying their complexation capability. [ABSTRACT FROM AUTHOR]

By extending an earlier study [Gianturco et al., J. Chem. Phys. 154, 054311 (2021)] on the purely rotational excitation of HeH+ by He atoms, we report in this paper integral cross sections and rate coefficients for rovibrational excitation and de-excitation processes in HeH+ due to collisions with He. The data were obtained using a new ab initio potential energy surface that includes the vibrational degree of freedom. The results are compared with those computed using the earlier potential energy surface by Panda and Sathyamurthy [J. Phys. Chem. A 107, 7125 (2003)] that additionally accounts for the proton-exchange reaction between HeH+ and He. It is shown that the exchange channel contributes nearly as much as the inelastic channel to the vibrational excitation and de-excitation processes and that the total rate constants pertaining to the purely inelastic processes are largely of the same magnitude as those obtained when both inelastic and reactive channels are included in the dynamics. The inelastic rovibrational rate coefficients involving this astrophysical cation are also found to be much larger than those obtained for anions present in similar interstellar environments. [ABSTRACT FROM AUTHOR]

The state-to-state (StS) model can accurately describe high-temperature thermochemical nonequilibrium flows. For the five-species air gas mixture, we develop a comprehensive database for the state-specific rate coefficients for temperatures 300–25 000 K in this paper. The database incorporates recent molecular dynamics simulations (based on the ab initio potential energy surfaces) in the literature, and theoretical methods, including the forced harmonic oscillator model and the Marrone–Treanor model, are employed to complement the rate coefficients that are unavailable from molecular dynamics calculations. The post-shock StS simulations using the present database agree with the experimental NO infrared radiation. Based on this updated StS kinetics database, we investigate the post-shock high-enthalpy air flows by employing both the StS and coarse-grained models (CGM). The CGM, which lumps molecular vibrational states into groups, shows results that align with the StS model, even utilizing only two groups for each molecule. However, the CGM-1G model, with only one group per molecule and belonging to the multi-temperature model (but uses StS kinetics), fails to reproduce the StS results. Analysis of vibrational energy source terms for different kinetic processes and fractions of vibrational groups reveals that the deficiency of the CGM-1G model stems from the overestimation of high-lying vibrational states, leading to higher dissociation rates and increased consumption of vibrational energy in dissociation. Furthermore, the presence of the Zeldovich-exchange processes indirectly facilitates energy transfer in N2 and O2, a phenomenon not observed in binary gas systems. These findings have important implications for developing the reduced-order model based on coarse-grained treatment. [ABSTRACT FROM AUTHOR]

A sample of NaCr(P2O7) was synthesized using the solid-state reaction method. X-ray diffraction and Rietveld refinement confirmed the formation of a monoclinic structure with the P21/a space group. Scanning electron microscopy (SEM) and energy dispersive X-ray (EDX) analysis identified the morphology and homogeneity of the chemical composition. The optical absorption spectrum revealed a direct optical band gap of 2.9 eV and an Urbach energy of 0.44 eV. Notably, the absorption spectrum showed an interference dip on the lowest energy band 4T2g(4F), attributed to spin--orbit coupling between states ²Eg(²G) and 4T2g(4F). This phenomenon was analyzed using the Neuhauser model based on coupled potential energy surfaces. From this analysis, the electronic structure of Cr3+ (3d3) ions in the NaCr(P2O7) sample was determined, enabling a reliable calculation of Racah and crystal-field parameters. The results showed good agreement between experimental and theoretical energy levels. The study demonstrates the impact of spin--orbit interactions using the coupled potential energy surface model. [ABSTRACT FROM AUTHOR]

For molecules with more than three atoms, it is difficult to fit or interpolate a potential energy surface (PES) from a small number of (usually ab initio) energies at points. Many methods have been proposed in recent decades, each claiming a set of advantages. Unfortunately, there are few comparative studies. In this paper, we compare neural networks (NNs) with Gaussian process (GP) regression. We re-fit an accurate PES of formaldehyde and compare PES errors on the entire point set used to solve the vibrational Schrödinger equation, i.e., the only error that matters in quantum dynamics calculations. We also compare the vibrational spectra computed on the underlying reference PES and the NN and GP potential surfaces. The NN and GP surfaces are constructed with exactly the same points, and the corresponding spectra are computed with the same points and the same basis. The GP fitting error is lower, and the GP spectrum is more accurate. The best NN fits to 625/1250/2500 symmetry unique potential energy points have global PES root mean square errors (RMSEs) of 6.53/2.54/0.86 cm−1, whereas the best GP surfaces have RMSE values of 3.87/1.13/0.62 cm−1, respectively. When fitting 625 symmetry unique points, the error in the first 100 vibrational levels is only 0.06 cm−1 with the best GP fit, whereas the spectrum on the best NN PES has an error of 0.22 cm−1, with respect to the spectrum computed on the reference PES. This error is reduced to about 0.01 cm−1 when fitting 2500 points with either the NN or GP. We also find that the GP surface produces a relatively accurate spectrum when obtained based on as few as 313 points. [ABSTRACT FROM AUTHOR]

Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree–Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree–Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single Ms (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree–Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis. [ABSTRACT FROM AUTHOR]

The emission of ESIPT-fluorophores is known to be sensitive to various external and internal stimuli and can be fine-tuned through substitution in the proton-donating and proton-accepting groups. The incorporation of metal ions in the molecules of ESIPT fluorophores without their deprotonation is an emerging area of research in coordination chemistry which provides chemists with a new factor affecting the ESIPT reaction and ESIPT-coupled luminescence. In this paper we present 1-hydroxy-5-methyl-4-(pyridin-2-yl)-2-(quinolin-2-yl)-1H-imidazole (HLq) as a new ESIPT-capable ligand. Due to the spatial separation of metal binding and ESIPT sites this ligand can coordinate metal ions without being deprotonated. The reactions of ZnHal2 with HLq afford ESIPT-capable [Zn(HLq)Hal2] (Hal = Cl, Br, I) complexes. In the solid state HLq and [Zn(HLq)Hal2] luminesce in the orange region (λmax = 600–650 nm). The coordination of HLq by Zn2+ ions leads to the increase in the photoluminescence quantum yield due to the chelation-enhanced fluorescence effect. The ESIPT process is barrierless in the S1 state, leading to the only possible fluorescence channel in the tautomeric form (T), S1T → S0T. The emission of [Zn(HLq)Hal2] in the solid state is blue-shifted as compared with HLq due to the stabilization of the ground state and destabilization of the excited state. In CH2Cl2 solutions, the compounds demonstrate dual emission in the UV (λmax = 358 nm) and green (λmax = 530 nm) regions. This dual emission is associated with two radiative deactivation channels in the normal (N) and tautomeric (T) forms, S1N → S0N and S1T → S0T, originating from two minima on the excited state potential energy surfaces. High energy barriers for the GSIPT process allow the trapping of molecules in the minimum of the tautomeric form, S0T, resulting in the possibility of the S0T → S1T photoexcitation and extraordinarily small Stokes shifts in the solid state. Finally, the π-system of quinolin-2-yl group facilitates the delocalization of the positive charge in the proton-accepting part of the molecule and promotes the ESIPT reaction. [ABSTRACT FROM AUTHOR]

In this paper, we analyze the influence of asymmetry on a Caldera potential energy surface. We first study the effect of asymmetry on the structure of the periodic orbit dividing surfaces associated with the unstable periodic orbits of the higher energy index-1 saddles. Then we detect a new type of dynamical matching due to the influence of the asymmetry in the Caldera potential energy surface. Finally, we present the phase space mechanism for this new type of dynamical matching using Lagrangian Descriptors. [ABSTRACT FROM AUTHOR]