Low-energy Structures Research Articles

J Coupling Constants of

Observing pyruvate metabolism in vivo has become a focal point of molecular magnetic resonance imaging. Signal amplification by reversible exchange (SABRE) has recently emerged as a versatile hyperpolarization technique. Tuning of the spin order transfer (SOT) in SABRE is challenging as the small 1H-13C J couplings, in the 13C-pyruvate case, result in SOT being not readily discernible. We demonstrate an experimental method using frequency-selective excitation of parahydrogen-derived polarization SOT sequence (SEPP-SPINEPT); its application led to up to 5700-fold 13C signal gain. In this way, we estimated the lifetime of two Ir-pyruvate SABRE complexes alongside the individual probing of eight small 1H-13C J couplings that connect the hydride protons in these complexes to 1- and 2-13C pyruvate spins, affording values between 0 and 2.69 Hz. Using electronic structure calculations, we define the low-energy structure of the corresponding complexes. Hence, this study demonstrates a novel approach to analyzing the spin topology of short-lived organometallic complexes.

Exploring the conformational landscape, hydrogen bonding, and internal dynamics in the diallyl ether and diallyl sulfide monohydrates.

The conformational spaces of the diallyl ether (DAE) and diallyl sulfide (DAS) monohydrates were explored using rotational spectroscopy from 6 to 19GHz. Calculations at the B3LYP-D3(BJ)/aug-cc-pVTZ level suggested significant differences in their conformational behavior, with DAE-w exhibiting 22 unique conformers and DAS-w featuring three stable structures within 6 kJ mol-1. However, only transitions from the lowest energy conformer of each were experimentally observed. Spectral analysis confirmed that binding with water does not alter the conformational preference for the lowest energy structure of the monomers, but it does influence the relative stabilities of all other conformers, particularly in the case of DAE. Non-covalent interaction and quantum theory of atoms in molecules analyses showed that the observed conformer for each complex is stabilized by two intermolecular hydrogen bonds (HBs), where water primarily interacts with the central oxygen or sulfur atom of the diallyl compounds, along with secondary interactions involving the allyl groups. The nature of these interactions was further elucidated using symmetry-adapted perturbation theory, which suggests that the primary HB interaction with S in DAS is weaker and more dispersive in nature compared to the primary HB in DAE. This supports the experimental observation of a tunneling splitting exclusively in the rotational spectrum of DAS-w, as the weaker contact allows water to undergo internal motions within the complex, as shown based on calculated transition state structures for possible tunneling pathways.

Flow-induced periodic chiral structures in an achiral nematic liquid crystal

Supramolecular chirality typically originates from either chiral molecular building blocks or external chiral stimuli. Generating chirality in achiral systems in the absence of a chiral input, however, is non-trivial and necessitates spontaneous mirror symmetry breaking. Achiral nematic lyotropic chromonic liquid crystals have been reported to break mirror symmetry under strong surface or geometric constraints. Here we describe a previously unrecognised mechanism for creating chiral structures by subjecting the material to a pressure-driven flow in a microfluidic cell. The chirality arises from a periodic double-twist configuration of the liquid crystal and manifests as a striking stripe pattern. We show that the mirror symmetry breaking is triggered at regions of flow-induced biaxial-splay configurations of the director field, which are unstable to small perturbations and evolve into lower energy structures. The simplicity of this unique pathway to mirror symmetry breaking can shed light on the requirements for forming macroscopic chiral structures.

Structural connectivity of 18‐Crown‐6/H+/L‐tryptophan noncovalent complexes in gas phase and in solution: Delineating host–guest–solvent interactions in solution from gas phase structures

AbstractWe investigate the structural correlation of noncovalent crown ether/H+/L‐tryptophan (CR/TrpH+) host–guest complexes in the solution phase with those in the gas phase generated through electrospray ionization/mass spectrometry (ESI/MS) techniques. We perform quantum chemical calculations to determine their structures, relative Gibbs free energies, and infrared spectra. We compare the calculated infrared (IR) spectra with the IR multiphoton dissociation (IRMPD) spectra observed for the 18‐Crown‐6/TrpH+ complex by Polfer and co‐workers [J. Phys. Chem. A 2013, 117, 1181–1188] for assigning the IR bands. We observe that the carboxyl group remains “naked,” lacking hydrogen bonding with the CR unit in the gas phase, and that this most stable conformer originates from the corresponding lowest Gibbs free energy structure in solution. Based on these findings, we propose that gas phase host–guest complexes directly correlate with those in solution, reinforcing the possibility of obtaining invaluable information about host–guest–solvent interactions in solution from the structure of the host–guest pair in the gas phase.

Revisiting a classic carbocation - DFT, coupled-cluster, and ab initio molecular dynamics computations on barbaralyl cation formation and rearrangements.

Density functional theory computations were used to model the formation and rearrangement of the barbaralyl cation (C9H+ 9). Two highly delocalized minima were located for C9H+ 9, one of C s symmetry and the other of D 3h symmetry, with the former having lower energy. Quantum chemistry-based NMR predictions affirm that the lower energy structure is the best match with experimental spectra. Partial scrambling was found to proceed through a C 2 symmetric transition structure associated with a barrier of only 2.3 kcal mol-1. The full scrambling was found to involve a C 2v symmetric transition structure associated with a 5.0 kcal mol-1 barrier. Ab initio molecular dynamics simulations initiated from the D 3h C9H+ 9 structure revealed its connection to six minima, due to the six-fold symmetry of the potential energy surface. The effects of tunneling and boron substitution on this complex reaction network were also examined.

Growth mechanism prediction for nanoparticles via structure matching polymerization.

Exploring structural and component evolution remains a challenging scientific problem for nanoscience. We propose a novel approach called principle of minimization of structure matching polymerization (SMP) change to rapidly explore the global minimum structure on the potential energy surface (PES). The new method can map low-dimensional stable structures to high-dimensional local minima, and this will make it possible for us to study the growth mechanisms of nanoparticles. Some new lowest-energy structures were found by SMP methods for sulfuric acid (SA)-dimethylamine (DMA) systems relative to previous studies. Additionally, we found that the growth process of boron clusters is mainly that the small-size boron clusters are continuously added to large-size boron clusters by structure matching for Bn (n = 2-36) systems, Bm + Bk → Bn, where m + k = n and 1 ≤ k ≤ 3. The SMP approach can greatly improve the search efficiency of other unbiased global optimization algorithms, such as basin-hopping (BH) and genetic algorithm (GA), with an enhancement of up to 19- and 7-fold relative to traditional BH and GA algorithms for searching the global minima of Bn (n = 14-22) systems. The SMP approach is general and flexible and can be applied to different kinds of problems, such as material structure design, crystal structure prediction, and new drug generation.

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach.

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

Last piece of the puzzle: lead-assisted assembly of an arsenopalladate nanostar as an antitumor metallodrug

Access to the long-sought nanostar of arsenopalladate is workable via host–guest assembly. The arsenate functions of such POP metallodrugs drive the ROS-induced apoptosis with great antitumor potential.

High-Temperature Superconductivity in La3Ni2O7

Motivated by the recent discovery of high-temperature superconductivity in bilayer La3Ni2O7 under pressure, we study its electronic properties and superconductivity due to strong electron correlation. Using the inversion symmetry, we decouple the low-energy electronic structure into block-diagonal symmetric and antisymmetric sectors. It is found that the antisymmetric sector can be reduced to a one-band system near half filling, while the symmetric bands occupied by about two electrons are heavily overdoped individually. Using the strong coupling mean field theory, we obtain strong superconducting pairing with B 1g symmetry in the antisymmetric sector. We propose that due to the spin-orbital exchange coupling between the two sectors, B 1g pairing is induced in the symmetric bands, which in turn boosts the pairing gap in the antisymmetric band and enhances the high-temperature superconductivity with a congruent d-wave symmetry in pressurized La3Ni2O7.

Beyond theory-driven discovery: introducing hot random search and datum-derived structures.

Data-driven methods have transformed the prospects of the computational chemical sciences, with machine-learned interatomic potentials (MLIPs) speeding up calculations by several orders of magnitude. I reflect on theory-driven, as opposed to data-driven, discovery based on ab initio random structure searching (AIRSS), and then introduce two new methods that exploit machine-learning acceleration. I show how long high-throughput anneals, between direct structural relaxation, enabled by ephemeral data-derived potentials (EDDPs), can be incorporated into AIRSS to bias the sampling of challenging systems towards low-energy configurations. Hot AIRSS (hot-AIRSS) preserves the parallel advantage of random search, while allowing much more complex systems to be tackled. This is demonstrated through searches for complex boron structures in large unit cells. I then show how low-energy carbon structures can be directly generated from a single, experimentally determined, diamond structure. An extension to the generation of random sensible structures, candidates are stochastically generated and then optimised to minimise the difference between the EDDP environment vector and that of the reference diamond structure. The distance-based cost function is captured in an actively learned EDDP. Graphite, small nanotubes and caged, fullerene-like, structures emerge from searches using this potential, along with a rich variety of tetrahedral framework structures. Using the same approach, the pyrope, Mg3Al2(SiO4)3, garnet structure is recovered from a low-energy AIRSS structure generated in a smaller unit cell with a different chemical composition. The relationship of this approach to modern diffusion-model-based generative methods is discussed.

Quantum computing for nuclear physics

Future quantum computers are anticipated to be able to perform simulations of quantum many-body systems and quantum field theories that lie beyond the capabilities of classical computation. This will lead to new insights and predictions for systems ranging from dense non-equilibrium matter, to low-energy nuclear structure and reactions, to high-energy collisions. I present an overview of digital quantum simulations in nuclear physics, with select examples relevant for studies of quark matter.

DLPNO-CCSD(T) and DFT study of the acetate-assisted C-H activation of benzaldimine at [RuCl2(p-cymene)]2: the relevance of ligand exchange processes at ruthenium(II) complexes in polar protic media.

To gain mechanistic insights into the acetate-assisted cyclometallations of arylimines promoted by [RuCl2(p-cymene)]2 in polar protic media, DFT geometry optimizations (with M06 and ωB97X-D3 functionals and the cc-pVDZ-PP[Ru] basis set) followed by DLPNO-CCSD(T)/CBS energy evaluations were performed using benzaldimine as a model substrate and methanol as the solvent (with CPCM or SMD models). The calculation results show that coordination of the imine to an acetate ruthenium precursor is followed by anion (chloride or acetate) dissociation as the rate-determining step of the process. H-Bonding of two explicit MeOH to the anion reduces the calculated activation energy to ca. 23 kcal mol-1, in good agreement with the experimental half-life at room temperature. Subsequent AMLA/CMD C-H activation of the intermediate cationic complexes is a faster, reversible process. Alternative reaction pathways involving neutral diacetate ruthenium complexes offer AMLA/CMD transition state structures of lower energy but are precluded due to higher energy barriers for the initial ligand exchange processes at ruthenium. Solvent assistance accelerates the final chloride/acetate exchange processes on the cycloruthenate intermediates, particularly when compression in the condensed phase is taken into consideration. The performance of six DFT functionals (with the aug-pVTZ-PP[Ru] basis set) was assessed using the DLPNO-CCSD(T)/CBS reference energies. Neutral diacetate ruthenium complexes were incorrectly predicted as being kinetically relevant when using hybrid DFT methods (PBE0-D3(BJ), M06-2X or ωB97M-V). Good agreement between the calculated barrier heights and our benchmark energy results was obtained by using double-hybrid DFT methods. PWPB95 with D3(BJ) or D4 dispersion energy corrections was found to be the most accurate (ΔG≠ MUE of ca. 1 kcal mol-1). This study may aid our understanding of and help with further experimental investigations of synthetically useful carboxylate-assisted C-H bond functionalizations involving (N,C)-cyclometallated (p-cymene)Ru(II) intermediate complexes in sustainable polar protic solvents.

Sub-picosecond photodynamics of small neutral copper oxide clusters.

The ultrafast dynamics of neutral copper oxide clusters (CunOx, n < 5) are reported using femtosecond pump probe spectroscopy in the gas phase. The transient spectra recorded for each cluster demonstrate they relax on a 100s of fs timescale followed by a long-lived (>50 ps) response. Density functional theory calculations are performed to determine the lowest energy structures and spin states. Topological descripters for the excited states are calculated (time-dependent density functional theory) to relate the measured excited state dynamics to changes in the cluster's electronic structure with increasing oxidation. Strong field ionization is demonstrated here to be a soft form of ionization and able to record transient signals for clusters previously determined to be unstable to nanosecond multiphoton ionization. The relative cluster stability is further demonstrated by signal enhancement/depreciation that is recorded through the synergy from the two laser pulses. Once the oxygen atoms exceed the number of copper atoms, a weakly bound superoxide O2 unit forms, exhibiting a higher spin state. All clusters that are not in the lowest spin configuration demonstrate fragmentation.

Infrared spectra and fragmentation dynamics of isotopologue-selective mixed-ligand complexes.

Isolated mixed-ligand complexes provide tractable model systems in which to study competitive and cooperative binding effects as well as controlled energy flow. Here, we report spectroscopic and isotopologue-selective infrared photofragmentation dynamics of mixed gas-phase Au(12/13CO)n(N2O)m+ complexes. The rich infrared action spectra, which are reproduced well using simulations of calculated lowest energy structures, clarify previous ambiguities in the assignment of vibrational bands, especially accidental coincidence of CO and N2O bands. The fragmentation dynamics exhibit the same unexpected behaviour as reported previously in which, once CO loss channels are energetically accessible, these dominate the fragmentation branching ratios, despite the much lower binding energy of N2O. We have investigated the dynamics computationally by considering anharmonic couplings between a relevant subset of normal modes involving both ligand stretch and intermolecular modes. Discrepancies between correlated and uncorrelated model fit to the ab initio potential energy curves are quantified using a Boltzmann sampled root mean squared deviation providing insight into efficiency of vibrational energy transfer between high frequency ligand stretches and the softer intermolecular modes which break during fragmentation.

Binuclear methylphosphinidine complexes of cyclopentadienylruthenium carbonyls: effects of the higher ligand field strength of ruthenium derivatives relative to iron derivatives.

The structures and energetics of the binuclear methylphosphinidene complexes of cyclopentadienylruthenium carbonyls of the type [MePRu2(CO)nCp2] (n = 4, 3, 2, 1) have been investigated for comparison with their previously studied iron analogues. For the tetracarbonyls and tricarbonyls [MePM2(CO)nCp2] (n = 4, 3) substituting ruthenium for iron has relatively little effect on the energetically preferred structures. Thus such structures have two-electron donor bridging MeP groups with no metal-metal bond for the tetracarbonyls and a metal-metal single bond for the tricarbonyls. This leads to favored 18-electron configurations for both ruthenium atoms in all cases. However, the higher ligand field strengths of ruthenium complexes relative to analogous iron complexes have major effects on the energetically preferred structures for the dicarbonyls and monocarbonyls [MePM2(CO)nCp2] (M = Fe, Ru; n = 2, 1). Thus the 11 lowest energy structures for the dicarbonyl [MePFe2(CO)2Cp2] are triplet or quintet spin state structures whereas the 6 lowest energy structures for the ruthenium analogue [MePRu2(CO)2Cp2] are all singlet structures. These low-energy singlet [MePRu2(CO)2Cp2] structures include species in which both ruthenium atoms attain the favored 18-electron configurations in different ways: either by a Ru-Ru single bond and an agostic C-H-Ru interaction from the methyl group, a Ru-Ru single bond and a four-electron donor bridging MeP ligand with PRu double bonds, or a formal RuRu double bond with a two-electron donor bridging MeP ligand. The 8 lowest energy structures for the diiron monocarbonyl [MePFe2(CO)Cp2] are all triplet or quintet spin structures whereas the lowest energy structure for the diruthenium monocarbonyl [MePRu2(CO)Cp2] by more than 20 kcal mol-1 is a singlet structure with a formal RuRu double bond and bridging CO and four-electron donor MeP groups. Thermochemical information predicts such monocarbonyl derivatives to be the dominant binuclear decarbonylation products of the tricarbonyls [RPRu2(CO)3Cp2] by thermal or photochemical methods.

Oligomers of diphenylalanine examined using cold ion spectroscopy and neural network-based conformational search.

Diphenylalanine (Phe2) is the primary building block of many self-assembling nanostructures that are important in biology and materials science. Understanding the detailed mechanism of their formation requires knowledge of the structural motives that the smallest oligomers attain at the very first steps of the process. Herein, we first employed high-resolution mass spectrometry to assign protonated Phe2 and its 2-13-unit oligomers formed in the gas phase from solution via electrospray ionization and then used cold ion spectroscopy to record UV and IR spectra for the monomer, dimer and hexamer. UV spectroscopy suggests the likely lack of specific strong proton-π interactions in oligomers larger than octamers, implying their certain structural stabilization. IR spectroscopy and quantum chemical calculations, enhanced by neural network-based conformational search, jointly determined the lowest-energy structures of the Phe2 monomer and dimer.

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics.

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N'-dimethyl, N,N'-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Hierarchical structures and magnetism of Co clusters: a perspective from integration of deep learning and a hybrid differential evolution algorithm.

Theoretically determining the lowest-energy structure of a cluster has been a persistent challenge due to the inherent difficulty in accurate description of its potential energy surface (PES) and the exponentially increasing number of local minima on the PES with the cluster size. In this work, density-functional theory (DFT) calculations of Co clusters were performed to construct a dataset for training deep neural networks to deduce a deep potential (DP) model with near-DFT accuracy while significantly reducing computational consumption comparable to classic empirical potentials. Leveraging the DP model, a high-efficiency hybrid differential evolution (HDE) algorithm was employed to search for the lowest-energy structures of CoN (N = 11-50) clusters. Our results revealed 38 of these clusters superior to those recorded in the Cambridge Cluster Database and identified diverse architectures of the clusters, evolving from layered structures for N = 11-27 to Marks decahedron-like structures for N = 28-42 and to icosahedron-like structures for N = 43-50. Subsequent analyses of the atomic arrangement, structural similarity, and growth pattern further verified their hierarchical structures. Meanwhile, several highly stable clusters, i.e., Co13, Co19, Co22, Co39, and Co43, were discovered by the energetic analyses. Furthermore, the magnetic stability of the clusters was verified, and a competition between the coordination number and bond length in affecting the magnetic moment was observed. Our study provides high-accuracy and high-efficiency prediction of the optimal structures of clusters and sheds light on the growth trend of Co clusters containing tens of atoms, contributing to advancing the global optimization algorithms for effective determination of cluster structures.

Edge-sharing water prisms.

Regular-shaped water clusters and nanostructures can be of particular interest for the self-assembly of complex structures and functional materials involving water molecules. Polyhedral water clusters are the most well studied. The low-energy structures of water hexamer, octamer and decamer are shaped like right prisms. The cavities of gas hydrates also have a polyhedral shape. Ice nanotubes are of no less interest both as configurations of global energy minima and in terms of possible applications. This article presents the results of the first systematic study of the structure and properties of a special class of water clusters. It reveals different combinations of three right prisms sharing an edge. According to modern calculations, the configurations of global energy minima of water clusters with the number of molecules being eighteen and twenty have this structure precisely. The topological characteristics of the edge-sharing water prisms are studied and the formulas for estimating the residual entropy are obtained. For these clusters (3-prisms), the usefulness of using a discrete model of intermolecular interaction [strong-weak-effective-bond model], previously developed for polyhedral water clusters, is shown. Calculations of the binding energy were performed using the non-additive Amoeba potential. To understand the relative stability of 3-prisms among other cluster forms, we use the structures recently reported in a published database containing over 3 × 106 unique water cluster networks (H2O)N of size N = 3-25 [Rakshit et al., J. Chem. Phys., 2019, 151(21), 214307, DOI: 10.1063/1.5128378].

H2 formation from the E2-E4 states of nitrogenase.

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 ×102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).