Journal of Chemical Physics最新文献

This work presents local hybrid alternatives to the orbital density approximation employed in self-interaction corrected density functional theory (SIC-DFT) and extended for use in DFT-corrected correlated wavefunction approaches (CAS-DFT). When combined with standard approximate density functionals, the orbital density approximation leaves SIC-DFT energies strongly dependent on unitary transforms among occupied orbitals and leaves CAS-DFT energies overbound. The alternatives presented here reduce both errors. The orbital density approximation and the local hybrid alternatives are shown to approximate an underlying nondiagonal exchange-correlation hole. A preliminary extension is presented to active-virtual correlation. These results motivate exploration of local hybrid concepts in SIC-DFT and CAS-DFT.

X-shaped liquid crystalline macromolecules (XLCMs) are obtained by tethering two flexible end A-blocks and two flexible side B-blocks to a semiflexible R-block. A rich array of ordered structures can be formed from XLCMs, driven by the competition between the interactions between the chemically distinct blocks and the molecular connectivity. Here, we report a theoretical study on the phase behavior of XLCMs with symmetric and asymmetric side blocks by using the self-consistent field theory (SCFT). A large number of ordered structures, including smectic phases, simple and giant polygons, are obtained as solutions of the SCFT equations. Phase diagrams of XLCMs as a function of the total length and asymmetric ratio of the side chains are constructed. For XLCMs with symmetric side blocks, the theoretically predicted phase transition sequence is in good agreement with experiments. For XLCMs with a fixed total side chain length, transitions between layered structure to polygonal phases, as well as between different polygonal phases, could be induced by varying the asymmetry of the side chains. The free energy density, domain size, side chain stretching, and molecular orientation are analyzed to elucidate mechanisms stabilizing the different ordered phases.

The computation of accurate geometric parameters at density functional theory cost for large molecules in the gas phase is addressed through a novel strategy that combines quantum chemical models with machine learning techniques. The first key step is the expansion of a database of accurate semi-experimental equilibrium structures with additional molecular geometries optimized by version 2 of the Pisa composite scheme. Then, the templating synthon approach is used to improve the accuracy of structures optimized by a hybrid density functional paired with a double zeta basis set, leveraging chemical similarity to cluster different molecular environments and refine bond lengths and valence angles. A set of prototypical biomolecular building blocks is used to demonstrate that it is possible to achieve spectroscopic accuracy for molecular systems too large to be treated by state-of-the-art composite wavefunction methods. In addition, a freely accessible web-based tool has been developed to facilitate the post-processing of geometries optimized using standard electronic structure codes, thereby providing an accurate and efficient tool for the computational study of medium- to large-sized molecules, also accessible to experiment-oriented researchers.

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange-correlation (XC) functionals spanning Jacob's ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. We further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron's absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron's spectroscopy, vertical detachment energy, or molar solvation volume.

Large-scale atomistic simulations rely on interatomic potentials, providing an efficient representation of atomic energies and forces. Modern machine-learning (ML) potentials provide the most precise representation compared to electronic structure calculations, while traditional potentials provide a less precise but computationally much faster representation and, thus, allow simulations of larger systems. We present a method to combine a traditional and a ML potential into a multi-resolution description, leading to an adaptive-precision potential with an optimum of performance and precision in large, complex atomistic systems. The required precision is determined per atom by a local structure analysis and updated automatically during simulation. We use copper as demonstrator material with an embedded atom model as classical force field and an atomic cluster expansion (ACE) as ML potential, but, in principle, a broader class of potential combinations can be coupled by this method. The approach is developed for the molecular-dynamics simulator LAMMPS and includes a load-balancer to prevent problems due to the atom dependent force-calculation times, which makes it suitable for large-scale atomistic simulations. The developed adaptive-precision copper potential represents the ACE-forces with a precision of 10 me V/Å and the ACE-energy exactly for the precisely calculated atoms in a nanoindentation of 4 × 106 atoms calculated for 100 ps and shows a speedup of 11.3 compared with a full ACE simulation.

We propose a reaction-limited evaporation model within the color-gradient lattice Boltzmann (LB) multicomponent framework to address the lack of intrinsic evaporation mechanisms. Unlike diffusion-driven approaches, our method directly enforces mass removal at the fluid interface in a reaction-limited manner while maintaining numerical stability. Using the inherent color-gradient magnitude and a single adjustable parameter, evaporation sites are chosen in a computationally efficient way with seamless mass exchange between the components, with no change to the core algorithm. Extensive validation across diverse interface geometries and evaporation flux magnitudes demonstrates high accuracy, with errors below 5% for unit density ratios. For density contrasts, the method remains robust in the limit of smaller evaporation flux magnitudes and density ratios. Our approach extends the applicability of the color-gradient LB model to scenarios involving reaction-limited evaporation, such as droplet evaporation on heated substrates, vacuum evaporation of molten metals, and drying processes in porous media.

Sparse data-driven approaches enable the approximation of governing laws of physical processes with parsimonious equations. While significant effort has been made in this field over the last decade, data-driven approaches generally rely on the paradigm of imposing a fixed base of library functions. In order to promote sparsity, finding the optimal set of basis functions is a necessary condition but a challenging task to guess in advance. Here, we propose an alternative approach that consists of optimizing the very library of functions while imposing sparsity. The robustness of our results is not only evaluated by the quality of the fit of the discovered model but also by the statistical distribution of the residuals with respect to the original noise in the data. In order to avoid choosing one metric over the other, we would rather rely on a multi-objective genetic algorithm (NSGA-II) for systematically generating a subset of optimal models sorted in a Pareto front. We illustrate how this method can be used as a tool to derive microkinetic equations from experimental data.

The nature of argentophilic interaction in the 2,2'-bipyridine-coordinated silver complex, which manifests counterintuitive cation-cation "attraction," is attributed to ligand stacking and solvation effects in the present article. While charged closed-shell transition metal complexes aggregating spontaneously to form oligomers has long been observed experimentally, the interpretation of the nature of so-called metallophilicity is still ongoing. For the dimer [(2,2'-bpy)2Ag]22+, qualitative electrostatic potential, non-covalent interaction, atoms-in-molecules analyses, and quantitative energy decomposition analysis calculations indicate that the electrostatic repulsion between two like formal charges at silver centers can be overcome by long-range dispersion attraction and short-range electronic correlation from ligands. In addition, delocalizing the net charges on silvers over the whole ligands can decrease electrostatic repulsion of metal centers to stabilize oligomers. The vital role of the screening effect of solvent has also been realized in the bound binding of the title system. Overall, this research highlights the importance of ligand stacking to argentophilicity, while d10-d10 attraction of silver centers presents quite little contribution.

A high level quantum mechanical study has been performed to explore the structural rearrangement and relative stability of the XH4+ (X = C, Si, Ge) radical cations at their X̃2T2 ground electronic states. All the stationary points located on the lowest adiabatic sheet of the Jahn-Teller (JT) split X̃2T2 state are fully optimized and characterized by performing harmonic vibrational frequency calculations. Five JT distorted stationary points with D2d(B22), C3v(A12), C2v(B22), and Cs(A'2) symmetries are located on the CH4+ ground state potential energy surface (PES), whereas four such structures are found on each of the SiH4+ and GeH4+ PESs. While the C2v(B22) isomer is found to be a global minimum and the Cs(A'2) one as a transition state for CH4+, the nature of them is reversed for SiH4+ and GeH4+. In particular, the Cs(A'2) stationary points are now global minima for the latter pair of radical cations, and C2v(B22) represents the transition state. Attempts are being made to understand such inconsistent findings via a combination of JT and epikernel principles. The barriers between equivalent C2v(B22) global minimum structures for CH4+ are found to be low, and thus CH4+ undergoes rapid interconversion along cyclic exchange of three hydrogen atoms via Cs transition state. The general features of the ground state PESs of SiH4+ and GeH4+ are similar. The pseudorotation between the Cs lowest energy structures undergoes along SiH2 and GeH2 wagging motions via C2v(B22) transition state for SiH4+ and GeH4+, respectively.

Doping boron clusters with metallic elements can tune the structural, electronic, and bonding properties. We report on the computational design of a zinc-rich D3h (1A1') B3Zn6- alloy cluster, whose global-minimum structure is a hybrid between prismatic, sandwich-like, and core-shell tubular geometries. The binary cluster features a linear B3 chain along its C3 axis, as well as three lateral Zn-Zn dimers, in which a central B atom is sandwiched by two quasi-planar BZn3 units in an eclipsed form. Chemical bonding analyses show that the B3 chain motif has Lewis-type B-B σ single bonds and a pair of orthogonal three-center two-electron (3c-2e) π bonds, collectively leading to a B-B bond order of two. Stabilizing a boron single chain is scarce in the literature, as is observing a series of double B=B bonds in a monoatomic chain fashion. The triangular pyramid BZn3 units are each in a unique triplet σ2σ*1σ*1 configuration, thus rendering σ aromaticity to the cluster according to the reversed 4n Hückel rule. It is proposed that the alloy cluster can be rationalized using the concept of electronic transmutation, wherein a close chemical analogy to the carbon dioxide (CO2) molecule is established.