Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “charge clusters”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Selective Detection of Diglycolamide Diastereomers By Metal-Ligand Cluster Formation

Introduction (120 words): Differentiation of isomers using mass spectrometry is a challenge, as these species tend to produce ions of identical mass-to-charge ratio. Therefore, separation with mass spectrometry must rely on chemical differences between isomeric species. Separation of stereoisomers is especially challenging, as these species often have very similar chemistry. However, recent investigations of the solvent extraction behavior of diastereomers of diglycolamide molecules have shown that different stereoisomers can have significantly different lanthanide extraction efficiencies. These studies showed that the R,S stereoisomer of 2,2'-oxybis(N,N-dioctylpropanamide efficiently extracts lanthanides, while the S,S stereoisomer does not. Thus, gas-phase metal complexation may provide a pathway for differentiating these species using mass spectrometry. Methods (120 words): Europium, samarium, and holmium clusters with the R,S, and the S,S diastereomers of 2,2'-oxybis(N,N-didecylpropanamide (trivially tetradecyl diglycolamide, or mTDDGA) were generated in gas-phase using the electrospray ionization source of a Bruker (Billerica, MA, USA) micrOTOF-Q II quadrupole time-of-flight mass spectrometer. Solutions of aqueous lanthanide nitrate and diglycolamide ligands in 2-propanol were added to 50%:50% by volume acetonitrile:2-propanol to generate the spray solutions. Ligand concentration was maintained at approximately 6 µM, while lanthanide concentration was maintained at either 6 µM or 30 µM . High resolution, high mass accuracy spectra were recorded using the time-of-flight. Mass accuracy was ensured using external calibration with Agilent (Santa Clara, CA, USA) ESI-L Low Concentration tuning mix. Preliminary data (300 words): Electrospray of solutions containing only R,S-mTDDGA or S,S-m-TDDGA yield abundant peaks at m/z= 721.7, corresponding to [mTDDGA+H]+1, m/z = 743.7, corresponding to [mTDDGA+Na]+1, and m/z = 1464.4, corresponding to the sodium-bound dimer, [Na(mTDDGA)2]+1. Spray of solutions containing 6 µM europium (Eu) with 3 µM R,S-mTDDGA produced an abundant cluster ion with an isotopic envelope centered at m/z=772.0, corresponding to [Eu(R,S-mTDDGA)3]3+. Europium-containing cluster ions with isotopic envelopes apexing at m/z = 1012.3 and 1189.0 correspond to the [Eu(R,S-mTDDGA)4]3+ and [Eu(R,S-mTDDGA)5]3+ cluster ions, respectively. Protonated and sodiated R,S-mTDDGA are present as well, along with the sodium-bound dimer. In contrast, spray of solutions containing 6 µM europium (Eu) with 3 µM S,S-mTDDGA diastereomer yielded only protonated and sodiated ligand ions; no Eu- S,S-mTDDGA ions were present. This mirrors the solution-phase behavior, where the S,S diestereomer extracts lanthanides with significantly less efficiency than the R,S diestereomer. Experiments with holmium and samarium yielded the same results. Increasing the lanthanide concentration to 30 µM resulted in limited formation of lanthanide-S,S-mTDDGA clusters. However, the abundance of the [Ln(S,S-mTDDGA)3]3+ clusters remained low compared to the abundance of the protonated and sodiated ligand. Additionally, the abundance of the [Ln(R,S-mTDDGA)3]3+ cluster ion was significantly higher than that of the protonated and sodiated ligand in experiments with spray solutions containing 30 µM lanthanide with R,S-mTDDA, supporting the hypothesis that R,S-mTDDGA has a higher intrinsic affinity for the lanthanides than S,S-mTDDGA. Thus, the intrinsic complexation with lanthanides appears to offer a route to differentiate diglycolamide diestereomers in the gas-phase. Novelty (20): Stereoisomers can be differentiated in the gas phase by exploitation of differences in metal complexation chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport of Light Ions in Matter

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NTJCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.

Wilson, J. W.↗

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an “optimal” value of the charge-transfer energy for the model where T c is largest. Finally, these results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural diversity and clustering of bacterial flagellar outer domains

Supercoiled flagellar filaments function as mechanical propellers within the bacterial flagellum complex, playing a crucial role in motility. Flagellin, the building block of the filament, features a conserved inner D0/D1 core domain across different bacterial species. In contrast, approximately half of the flagellins possess additional, highly divergent outer domain(s), suggesting varied functional potential. In this study, we report atomic structures of flagellar filaments from three distinct bacterial species: Cupriavidus gilardii , Stenotrophomonas maltophilia , and Geovibrio thiophilus . Our findings reveal that the flagella from the facultative anaerobic G. thiophilus possesses a significantly more negatively charged surface, potentially enabling adhesion to positively charged minerals. Furthermore, we analyze all AlphaFold predicted structures for annotated bacterial flagellins, categorizing the flagellin outer domains into 682 structural clusters. This classification provides insights into the prevalence and experimental verification of these outer domains. Remarkably, two of the flagellar structures reported herein belong to a distinct cluster, indicating additional opportunities on the study of the functional diversity of flagellar outer domains. Our findings underscore the complexity of bacterial flagellins and open up possibilities for future studies into their varied roles beyond motility.

Science & Technology - Other Topics↗

Discovery of a Peptoid-Based Nanoparticle Platform for Therapeutic mRNA Delivery via Diverse Library Clustering and Structural Parametrization

Nanoparticle-mediated mRNA delivery has emerged as a promising therapeutic modality, but its growth is still limited by the discovery and optimization of effective and well-tolerated delivery strategies. Lipid nanoparticles containing charged or ionizable lipids are an emerging standard for in vivo mRNA delivery, so creating facile, tunable strategies to synthesize these key lipid-like molecules is essential to advance the field. Here, we generate a library of N-substituted glycine oligomers, peptoids, and undertake a multistage down-selection process to identify lead candidate peptoids as the ionizable component in our Nutshell nanoparticle platform. First, we identify a promising peptoid structural motif by clustering a library of >200 molecules based on predicted physical properties and evaluate members of each cluster for reporter gene expression in vivo. Then, the lead peptoid motif is optimized using design of experiments methodology to explore variations on the charged and lipophilic portions of the peptoid, facilitating the discovery of trends between structural elements and nanoparticle properties. We further demonstrate that peptoid-based Nutshells leads to expression of therapeutically relevant levels of an anti-respiratory syncytial virus antibody in mice with minimal tolerability concerns or induced immune responses compared to benchmark ionizable lipid, DLin-MC3-DMA. Through this work, we present peptoid-based nanoparticles as a tunable delivery platform that can be optimized toward a range of therapeutic programs.

59 BASIC BIOLOGICAL SCIENCES↗

Exotic bonding interactions and coexistence of chemically distinct periodic lattice distortions in the charge density wave compound TaTe2

The superstructure of TaTe 2 is examined by total high-energy x-ray scattering and large-scale structure modeling over a broad temperature range. At room temperature, it features double zigzag chains of Ta atoms exhibiting metallic bonding. The chains are sandwiched between planes of weakly interacting Te atoms. Below 170 K, the chains appear broken to “butterflylike” clusters and single zigzag chains of covalently bonded Te atoms emerge while Ta-Te interactions remain largely ionic. This is a rare example of a charge density wave compound stabilized by bonding that is a mix of covalent, ionic, and metallic character giving rise to chemically distinct periodic lattice distortions. We argue that the formation of periodic lattice distortions in the Te sublattice in addition to those in the Ta sublattice favors charge delocalization, which may explain the unusually diminished resistivity and increased magnetic susceptibility exhibited by the low-temperature charge density wave phase of TaTe 2 . The effect may be common for 5d transition metal chalcogenides containing Te because of the extended electronic orbitals of the constituent atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rapid-charging aluminium-sulfur batteries operated at 85 °C with a quaternary molten salt electrolyte

Molten salt aluminum-sulfur batteries are based exclusively on resourcefully sustainable materials, and are promising for large-scale energy storage owed to their high-rate capability and moderate energy density; but the operating temperature is still high, prohibiting their applications. Here we report a rapid-charging aluminium-sulfur battery operated at a sub-water-boiling temperature of 85°C with a tamed quaternary molten salt electrolyte. The quaternary alkali chloroaluminate melt – possessing abundant electrochemically active high-order Al-Cl clusters and yet exhibiting a low melting point – facilitates fast Al 3+ desolvation. A nitrogen-functionalized porous carbon further mediates the sulfur reaction, enabling the battery with rapid-charging capability and excellent cycling stability with 85.4% capacity retention over 1400 cycles at a charging rate of 1 C. Importantly, we demonstrate that the asymmetric sulfur reaction mechanism that involves formation of polysulfide intermediates, as revealed by operando X-ray absorption spectroscopy, accounts for the high reaction kinetics at such temperature wherein the thermal management can be greatly simplified by using water as the heating media.

25 ENERGY STORAGE↗

Pyrotechnic hazards classification and evaluation program. Electrostatic vulnerability of the E8 and XM15/XM165 clusters, phase 2

An investigation into the electrostatic phenomena associated with the manufacturing and handling of explosives is discussed. The testing includes measurement of the severity of the primary charge generation mechanism, triboelectric effects between dissimilar surfaces; refinement of equivalent circuits of the XM15/XM165 and E8 fuse trains; evaluation of the electrostatic spark discharge characteristics predicted by an equivalent circuit analysis; and determination of the spark ignition sensitivity of materials, components, junctions, and subassemblies which compose the XM15/XM165 and E8 units. Special studies were also performed. These special tests included ignition sensitivity of the complete XM15 fuse train when subjected to discharges through its entire length, measurement of electrostatic potentials which occur during the E8 foaming operation during fabrication, and investigation of the inadvertent functioning of an XM15 cluster during manufacturing. The test results are discussed and related to the effectiveness of suggested modification to reduce the electrostatic ignition sensitivity.

Source record↗

Can Polyelectrolyte Mechanical Properties be Directly Modulated by an Electric Field? A Molecular Dynamics Study

Multifunctionality in polymers facilitates their application in emerging technologies. Electrical fields are a preferred stimulus because of the speed and ease of application to bulk polymers. While a wide range of electrically triggered actuators are developed, and electrically controlled adhesion between gels is demonstrated, modification of bulk mechanical properties via electrical stimuli remains elusive. Polymers with covalently incorporated ionic charge (polyelectrolytes) should be well suited to achieving this goal since the mechanical properties depend on electrostatic interactions and these charges are intrinsically susceptible to electric fields. Molecular dynamics simulations are utilized here to investigate whether electric fields can modulate the mechanical properties of polyelectrolytes and to understand the governing mechanisms. Mechanical property modulation by electric field is found to be sensitive to the charge distribution—charges must be tightly attached to the polymer backbone, and responsivity is greater if a single backbone contains both positive and negative charges. The dominant mechanisms are reorientation and stretching of the polymer chains, which also elongate the ionic clusters to maintain strong electrostatic interactions throughout deformation. Furthermore, these insights are critical for future experimental realization of polymers with electric field regulated mechanical properties.

36 MATERIALS SCIENCE↗

QCD hidden-color hexadiquark in the core of nuclei

Hidden-color configurations are a key prediction of QCD with important physical consequences. In this work we examine a QCD color-singlet configuration in nuclei formed by combining six scalar [ud] diquarks in a strongly bound SU (3) C channel. The resulting hexadiquark state is a charge-2, spin-0, baryon number-4, isospin-0, color-singlet state. It contributes to alpha clustering in light nuclei and to the additional binding energy not saturated by ordinary nuclear forces in as well as the alpha-nuclei sequence of interest for nuclear astrophysics. We show that the strongly bound combination of six scalar isospin-0 [ud] diquarks within the nuclear wave function - relative to free nucleons - provides a natural explanation of the EMC effect measured by the CLAS collaboration's comparison of nuclear parton distribution function ratios for a large range of nuclei. Overall, these experiments confirmed that the EMC effect; i.e., the distortion of quark distributions within nuclei, is dominantly identified with the dynamics of neutron-proton (“isophobic”) short-range correlations within the nuclear wave function rather than proton-proton or neutron-neutron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quasi-Two-Dimensional Heterostructures (K M 1 – x Te)(LaTe 3 ) ( M = Mn and Zn) with Charge Density Waves

Layered heterostructure materials with two different functional building blocks can teach us about emergent physical properties and phenomena arising from interactions between the layers. Here, we report intergrowth compounds KLaM 1 - x Te 4 (M = Mn and Zn; $x \approx$ 0.35) featuring two chemically distinct alternating layers [LaTe 3 ] and [KM 1 - x Te]. Their crystal structures are incommensurate, determined by single X-ray diffraction for the Mn compound and a transmission electron microscope study for the Zn compound. KLaMn 1 – x Te 4 crystallizes in the orthorhombic superspace group Pmnm(01/2 gamma)s00 with lattice parameters a = 4.4815(3) Å, b = 21.6649(16) Å, and c = 4.5220(3) Å. It exhibits charge density wave order at room temperature with a modulation wave vector q = 1/2b* + 0.3478c* originating from electronic instability of Te-square nets in [LaTe 3 ] layers. The Mn analog exhibits a cluster spin glass behavior with spin freezing temperature $T_f \approx$ 5 K attributed to disordered Mn vacancies and competing magnetic interactions in the [Mn 1 - x Te] layers. The Zn analog also has charge density wave order at room temperature with a similar q-vector having the c* component similar to 0.346 confirmed by selected-area electron diffraction. Electron transfer from [KM 1 - x Te] to [LaTe 3 ] layers exists in KLaM 1 – x Te 4 , leading to an enhanced electronic specific heat coefficient. The resistivities of KLaM 1 - x Te 4 (M = Mn and Zn) exhibit metallic behavior at high temperatures and an upturn at low temperatures, suggesting partial localization of carriers in the [LaTe 3 ] layers with some degree of disorder associated with the M atom vacancies in the [M 1 - x Te] layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

Completely Multipolar Model for Many-Body Water–Ion and Ion–Ion Interactions

This work constructs an advanced force field, the Completely Multipolar Model (CMM), to quantitatively reproduce each term of an energy decomposition analysis (EDA) for aqueous solvated alkali metal cations and halide anions and their ion pairings. We find that all individual EDA terms remain well-approximated in the CMM for ion-water and ion-ion interactions, except for polarization, which shows errors due to the partial covalency of ion interactions near their equilibrium. We quantify the onset of the dative bonding regime by examining the change in molecular polarizability and Mayer bond indices as a function of distance, showing that partial covalency manifests by breaking the symmetry of atomic polarizabilities while strongly damping them at short-range. This motivates an environment-dependent atomic polarizability parameter that depends on the strength of the local electric field experienced by the ions to account for strong damping, with anisotropy introduced by atomic multipoles. The resulting CMM model for ions provides accurate dimer surfaces and three-body polarization and charge transfer compared to EDA, and shows excellent performance on various ion benchmarks including vibrational frequencies and cluster geometries.

Heindel, Joseph P↗

Cation valency in water-in-salt electrolytes alters the short- and long-range structure of the electrical double layer

Highly concentrated aqueous electrolytes (termed water-in-salt electrolytes, WiSEs) at solid-liquid interfaces are ubiquitous in myriad applications including biological signaling, electrosynthesis, and energy storage. This interface, known as the electrical double layer (EDL), has a different structure in WiSEs than in dilute electrolytes. Here, we investigate how divalent salts [zinc bis(trifluoromethylsulfonyl)imide, Zn(TFSI) 2 ], as well as mixtures of mono- and divalent salts [lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) mixed with Zn(TFSI) 2 ], affect the short- and long-range structure of the EDL under confinement using a multimodal combination of scattering, spectroscopy, and surface forces measurements. Raman spectroscopy of bulk electrolytes suggests that the cation is closely associated with the anion regardless of valency. Wide-angle X-ray scattering reveals that all bulk electrolytes form ion clusters; however, the clusters are suppressed with increasing concentration of the divalent ion. To probe the EDL under confinement, we use a Surface Forces Apparatus and demonstrate that the thickness of the adsorbed layer of ions at the interface grows with increasing divalent ion concentration. Multiple interfacial layers form following this adlayer; their thicknesses appear dependent on anion size, rather than cation. Importantly, all electrolytes exhibit very long electrostatic decay lengths that are insensitive to valency. It is likely that in the WiSE regime, electrostatic screening is mediated by the formation of ion clusters rather than individual well-solvated ions. This work contributes to understanding the structure and charge-neutralization mechanism in this class of electrolytes and the interfacial behavior of mixed-electrolyte systems encountered in electrochemistry and biology.

Science & Technology - Other Topics↗

Kuramoto synchronization of quantum tunneling polarons for describing the dynamic structure in cuprate superconductors

A major open topic in cuprates is the interplay between the lattice and electronic dynamics and the importance of their coupling to the mechanism of high-temperature superconductivity (HTSC). As evidenced by extended x-ray absorption fine structure (EXAFS) experiments, anharmonic structural effects are correlated with the charge dynamics and the transition to a superconducting phase in different HTSC compounds. Here, in this paper, we describe how structural anharmonic effects can be coupled to electronic and lattice dynamics in cuprate systems by performing the exact diagonalization of a prototype anharmonic many-body Hamiltonian on a relevant 6-atom cluster and show that the EXAFS results can be understood as a Kuramoto synchronization transition between coupled internal quantum tunneling of polarons associated with the two-site distribution of the copper–apical oxygen (Cu-O ap ) pair in the dynamic structure. The transition is driven by the anharmonicity of the lattice vibrations and promotes the pumping of charge, initially stored at the apical oxygen reservoirs, into the copper-oxide plane. Simultaneously, a finite projection of the internal quantum tunneling polaron extends to the copper–planar oxygen (Cu-O pl ) pair. All these findings allow an interpretation based on an effective quantum-mechanical triple-well potential associated with the oxygen sites of the 6-atom cluster, which accurately represents the phase synchronization of apical oxygens and lattice-assisted charge transfer to the CuO 2 plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First Observation of the β3αp Decay of 13 O via β-Delayed Charged-Particle Spectroscopy

The β-delayed proton decay of 13 O has previously been studied, but the direct observation of β-delayed 3⁢α⁢p decay has not been reported. Rare 3⁢α⁢p events from the decay of excited states in 13 N* provide a sensitive probe of cluster configurations in 13 N*. To measure the low-energy products following β-delayed 3⁢αp decay, the Texas Active Target (TexAT) time projection chamber was employed using the one-at-a-time β-delayed charged-particle spectroscopy technique at the Cyclotron Institute, Texas A&M University. A total of 1.9 × 10 5 13 O implantations were made inside the TexAT time projection chamber. Furthermore, a total of 149 3⁢αp events were observed, yielding a β-delayed 3⁢αp branching ratio of 0.078(6)%. Four previously unknown α-decaying excited states were observed in 13 N at 11.3, 12.4, 13.1, and 13.7 MeV decaying via the 3⁢α + p channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structural evidence for intermediates during O 2 formation in photosystem II

In natural photosynthesis, the light-driven splitting of water into electrons, protons and molecular oxygen forms the first step of the solar-to-chemical energy conversion process. The reaction takes place in photosystem II, where the Mn 4 CaO 5 cluster first stores four oxidizing equivalents, the S 0 to S 4 intermediate states in the Kok cycle, sequentially generated by photochemical charge separations in the reaction center and then catalyzes the O–O bond formation chemistry. Here, we report room temperature snapshots by serial femtosecond X-ray crystallography to provide structural insights into the final reaction step of Kok’s photosynthetic water oxidation cycle, the S 3 →[S 4 ]→S 0 transition where O 2 is formed and Kok’s water oxidation clock is reset. Our data reveal a complex sequence of events, which occur over micro- to milliseconds, comprising changes at the Mn 4 CaO 5 cluster, its ligands and water pathways as well as controlled proton release through the hydrogen-bonding network of the Cl1 channel. Importantly, the extra O atom O x , which was introduced as a bridging ligand between Ca and Mn1 during the S 2 →S 3 transition disappears or relocates in parallel with Y z reduction starting at approximately 700 μs after the third flash. The onset of O 2 evolution, as indicated by the shortening of the Mn1–Mn4 distance, occurs at around 1,200 μs, signifying the presence of a reduced intermediate, possibly a bound peroxide.

59 BASIC BIOLOGICAL SCIENCES↗

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗