Photochemical Implications of Changes in the Spectral Properties of Chromophoric Dissolved Organic Matter: A Model Assessment for Surface Waters

Chromophoric dissolved organic matter (CDOM) is the main sunlight absorber in surface waters and a very important photosensitiser towards the generation of photochemically produced reactive intermediates (PPRIs), which take part in pollutant degradation. The absorption spectrum of CDOM (ACDOM(λ), unitless) can be described by an exponential function that decays with increasing wavelength: ACDOM(λ) = 100 d DOC Ao e− Sλ, where d [m] is water depth, DOC [mgC L−1] is dissolved organic carbon, Ao [L mgC−1 cm−1] is a pre-exponential factor, and S [nm−1] is the spectral slope. Sunlight absorption by CDOM is higher when Ao and DOC are higher and S is lower, and vice versa. By the use of models, here we investigate the impact of changes in CDOM spectral parameters (Ao and S) on the steady-state concentrations of three PPRIs: the hydroxyl radical (•OH), the carbonate radical (CO3•−), and CDOM excited triplet states (3CDOM*). A first finding is that variations in both Ao and S have impacts comparable to DOC variations on the photochemistry of CDOM, when reasonable parameter values are considered. Therefore, natural variability of the spectral parameters or their modifications cannot be neglected. In the natural environment, spectral parameters could, for instance, change because of photobleaching (prolonged exposure of CDOM to sunlight, which decreases Ao and increases S) or of the complex and still poorly predictable effects of climate change. A second finding is that, while the steady-state [3CDOM*] would increase with increasing ACDOM (increasing Ao, decreasing S), the effect of spectral parameters on [•OH] and [CO3•−] depends on the relative roles of CDOM vs. NO3− and NO2− as photochemical •OH sources.

Computation of entropy values for non-electrolyte solute molecules in solution based on semi-empirical corrections to a polarized continuum model

A simple heuristic model was developed for estimating the entropy of a solute molecule in an ideal solution based on quantum mechanical calculations with polarizable continuum models (QM/PCMs). A translational term was incorporated that included free-volume compensation for the Sackur-Tetrode equation and a rotational term was modeled based on the restricted rotation of a dipole in an electrostatic field. The configuration term for the solute at a given concentration was calculated using a simple lattice model that considered the number of configurations of the solute within the lattice. The configurational entropy was ascertained from this number based on Boltzmann's principle. Standard entropy values were determined for 41 combinations of solutes and solvents at a set concentration of 1 mol dm^-3 using the proposed model, and the computational values were compared with experimental data. QM/PCM calculations were conducted at the ωB97X-D/6-311++G(d,p)/IEF-PCM level using universal force field van der Waals radii scaled by 1.2. The proposed model accurately reproduced the entropy values reported for solutes in non-aqueous solvents within a mean absolute deviation of 9.2 J mol^-1 K^-1 for 33 solutions. This performance represents a considerable improvement relative to that obtained using the method based on the ideal gas treatment that is widely utilized in commercially available computation packages. In contrast, computations for aqueous molecules overestimated the entropies because hydrophobic effects that decrease the entropy of aqueous solutions were not included in the present model.

Accurate determination of reaction energetics and kinetics of the HO(2)˙ + O(3) → OH˙ + 2O(2) reaction

In the present work, we have studied the HO₂˙ + O₃ → HO˙ + 2O₂ reaction using chemical kinetics and quantum chemical calculations. We have employed the post-CCSD(T) method to estimate the barrier height and reaction energy for the title reaction. In the post-CCSD(T) method, we have included zero point energy corrections, contributions from full triple excitations and partial quadratic excitations at the coupled-cluster level, and core corrections. We have also computed the reaction rate in the temperature range of 197-450 K and found good agreement with all the available experimental results. In addition, we have also fitted the computed rate constants with the Arrhenius expression and obtained an activation energy of 1.0 ± 0.1 kcal mol^-1, almost identical to the value recommended by IUPAC and JPL.

1-Azahomocubane

Highly strained cage hydrocarbons have long stood as fundamental molecules to explore the limits of chemical stability and reactivity, probe physical properties, and more recently as bioactive molecules and in materials discovery. Interestingly, the nitrogenous congeners have attracted much less attention. Previously absent from the literature, azahomocubanes, offer an opportunity to investigate the effects of a nitrogen atom when incorporated into a highly constrained polycyclic environment. Herein disclosed is the synthesis of 1-azahomocubane, accompanied by comprehensive structural characterization, physical property analysis and chemical reactivity. These data support the conclusion that nitrogen is remarkably well tolerated in a highly strained environment.

Chemical similarity of dialkyl carbonates and carbon dioxide opens an avenue for novel greenhouse gas scavengers: cheap recycling and low volatility via experiments and simulations

Global warming linked to the industrial emissions of greenhouse gases may be the end of mankind unless it is adequately and timely handled. To prevent irreversible changes to the climate of the Earth, numerous research groups are striving to develop robust CO₂ sorbents. Dialkyl carbonates (DACs) and CO₂ exhibit obvious chemical similarities in their structure and properties. The degrees of oxidation of all atoms composing DACs and CO₂ are identical resulting in very similar nucleophilicities and electrophilicities of all interaction centers. While both compounds possess relatively high partial atomic charges on their polar moieties, the molecular geometries prevent tight binding of the head groups. The computed DAC-DAC binding energies are ∼40 kJ mol^-1, whereas the effect of the alkyl chain length is marginal. The phase transition points and shear viscosities of DACs are very low. We herein hypothesize and numerically rationalize that DACs represent noteworthy physical sorbents for CO₂ thanks to the similar sorbent-CO₂ and sorbent-sorbent interaction energies. By reporting in silico-derived sorption thermodynamics at various conditions, spectral and structural properties, and experimentally derived CO₂ capacities and recyclabilities, we highlight the mutual affinity of DACs and CO₂. Indeed, the experimentally determined CO₂ sorption capacity of 0.88 mol% (diethyl carbonate) at 278.15 K and 30 bar is competitive. The unprecedentedly low DAC-CO₂ binding energies, ∼14 kJ mol^-1, suggest a low-cost desorption process and outstanding recyclability of the sorbent. We also note that DACs possessing long alkyl chains (butyl, hexyl, octyl) exhibit negligible volatilities, while preserving the liquid aggregate state over a practically important temperature range. The reported results may foster the development of a new class of CO₂ scavengers with possibly quite peculiar characteristics.

Study on Benzylamine(BZA) and Aminoethylpiperazine(AEP) Mixed Absorbent on Ship-Based Carbon Capture

To find suitable absorbents for ship-based carbon capture, the absorption and desorption properties of four mixed aqueous amines based on BZA were investigated, and the results indicated that BZA-AEP had the best absorption and desorption performance. Then, the absorption and desorption properties of different mole ratios of BZA-AEP were tested. The results showed that the average CO2 absorption rate had the highest value at the mole ratio of BZA to AEP of three. The average CO2 desorption rate had the maximum value at the mole ratio of BZA to AEP of one. Three fitted models of the absorption and desorption performance of BZA-AEP based on the test data were obtained. The p-values of all three models were less than 0.0001. Considering the performance and material cost, the BZA-AEP mole ratio of 1.5 is more appropriate for ship carbon capture. Compared with MEA, the average CO2 absorption rate increased by 48%, the CO2 desorption capacity increased by 120%, and the average CO2 desorption rate increased by 161%.

Kinetics and Mechanisms of Pressure-Induced Ice Amorphization and Polyamorphic Transitions in a Machine-Learned Coarse-Grained Water Model

Water glasses have attracted considerable attention due to their potential connection to a liquid-liquid transition in supercooled water. Here we use molecular simulations to investigate the formation and phase behavior of water glasses using the machine-learned bond-order parameter (ML-BOP) water model. We produce glasses through hyperquenching of water, pressure-induced amorphization (PIA) of ice, and pressure-induced polyamorphic transformations. We find that PIA of polycrystalline ice occurs at a lower pressure than that of monocrystalline ice and through a different mechanism. The temperature dependence of the amorphization pressure of polycrystalline ice for ML-BOP agrees with that in experiments. We also find that ML-BOP accurately reproduces the density, coordination number, and structural features of low-density (LDA), high-density (HDA), and very high-density (VHDA) amorphous water glasses. ML-BOP accurately reproduces the experimental radial distribution function of LDA but overpredicts the minimum between the first two shells in high-density glasses. We examine the kinetics and mechanism of the transformation between low-density and high-density glasses and find that the sharp nature of these transitions in ML-BOP is similar to that in experiments and all-atom water models with a liquid-liquid transition. Transitions between ML-BOP glasses occur through a spinodal-like mechanism, similar to ice crystallization from LDA. Both glass-to-glass and glass-to-ice transformations have Avrami-Kolmogorov kinetics with exponent n = 1.5 ± 0.2 in experiments and simulations. Importantly, ML-BOP reproduces the competition between crystallization and HDA→LDA transition above the glass transition temperature T_g, and separation of their time scales below T_g, observed also in experiments. These findings demonstrate the ability of ML-BOP to accurately reproduce water properties across various regimes, making it a promising model for addressing the competition between polyamorphic transitions and crystallization in water and solutions.

Heart-Type Fatty Acid Binding Protein Binds Long-Chain Acylcarnitines and Protects against Lipotoxicity

Heart-type fatty-acid binding protein (FABP3) is an essential cytosolic lipid transport protein found in cardiomyocytes. FABP3 binds fatty acids (FAs) reversibly and with high affinity. Acylcarnitines (ACs) are an esterified form of FAs that play an important role in cellular energy metabolism. However, an increased concentration of ACs can exert detrimental effects on cardiac mitochondria and lead to severe cardiac damage. In the present study, we evaluated the ability of FABP3 to bind long-chain ACs (LCACs) and protect cells from their harmful effects. We characterized the novel binding mechanism between FABP3 and LCACs by a cytotoxicity assay, nuclear magnetic resonance, and isothermal titration calorimetry. Our data demonstrate that FABP3 is capable of binding both FAs and LCACs as well as decreasing the cytotoxicity of LCACs. Our findings reveal that LCACs and FAs compete for the binding site of FABP3. Thus, the protective mechanism of FABP3 is found to be concentration dependent.

Infrared spectra of isoquinolinium (iso-C(9)H(7)NH(+)) and isoquinolinyl radicals (iso-C(9)H(7)NH and 1-, 3-, 4-, 5-, 6-, 7- and 8-iso-HC(9)H(7)N) isolated in solid para-hydrogen

Protonated polycyclic aromatic nitrogen heterocycles (H⁺PANH) are prospective candidates that may contribute to interstellar unidentified infrared (UIR) emission bands because protonation enhances the relative intensities of the bands near 6.2, 7.7 and 8.6 μm, and the presence of the N atom induces a blue shift of the ring-stretching modes so that the spectra of H⁺PANH match better with the 6.2 μm feature in class-A UIR spectra. We report the infrared (IR) spectra of protonated isoquinoline (the 2-isoquinolinium cation, iso-C₉H₇NH⁺), its neutral counterpart (the 2-isoquinolinyl radical, iso-C₉H₇NH), and another mono-hydrogenated product (the 6-isoquinolinyl radical, 6-iso-HC₉H₇N), produced on the electron-bombardment of a mixture of isoquinoline (iso-C₉H₇N) with excess para-hydrogen (p-H₂) during matrix deposition at 3.2 K. To generate additional isomers of hydrogenated isoquinoline, we irradiated iso-C₉H₇N/Cl₂/p-H₂ matrices at 365 nm to generate Cl atoms, followed by IR irradiation to generate H atoms via Cl + H₂ (v = 1) → HCl + H; the H atoms thus generated reacted with iso-C₉H₇N. In addition to iso-C₉H₇NH and 6-iso-HC₉H₇N observed in the electron-bombardment experiments, we identified six additional hydrogenated isoquinoline species, 1-, 3-, 4-, 5-, 7- and 8-iso-HC₉H₇N, via their IR spectra; hydrogenation on the N atom and all available carbon atoms except for the two sharing carbon atoms on the fused ring was observed. Spectral groupings were achieved according to their behaviors after maintenance of the matrix in darkness and on secondary photolysis at various wavelengths. The assignments were supported via comparison of the experimental results with the vibrational wavenumbers and IR intensities of possible isomers predicted using the B3LYP/6-311++G(d,p) method. The implications in the identification of the UIR band are discussed.

Classical density functional theory for interfacial properties of hydrogen, helium, deuterium, neon, and their mixtures

We present a classical density functional theory (DFT) for fluid mixtures that is based on a third-order thermodynamic perturbation theory of Feynman-Hibbs-corrected Mie potentials. The DFT is developed to study the interfacial properties of hydrogen, helium, neon, deuterium, and their mixtures, i.e., fluids that are strongly influenced by quantum effects at low temperatures. White Bear fundamental measure theory is used for the hard-sphere contribution of the Helmholtz energy functional, and a weighted density approximation is used for the dispersion contribution. For mixtures, a contribution is included to account for non-additivity in the Lorentz-Berthelot combination rule. Predictions of the radial distribution function from DFT are in excellent agreement with results from molecular simulations, both for pure components and mixtures. Above the normal boiling point and 5% below the critical temperature, the DFT yields surface tensions of neon, hydrogen, and deuterium with average deviations from experiments of 7.5%, 4.4%, and 1.8%, respectively. The surface tensions of hydrogen/deuterium, para-hydrogen/helium, deuterium/helium, and hydrogen/neon mixtures are reproduced with a mean absolute error of 5.4%, 8.1%, 1.3%, and 7.5%, respectively. The surface tensions are predicted with an excellent accuracy at temperatures above 20 K. The poor accuracy below 20 K is due to the inability of Feynman-Hibbs-corrected Mie potentials to represent the real fluid behavior at these conditions, motivating the development of new intermolecular potentials. This DFT can be leveraged in the future to study confined fluids and assess the performance of porous materials for hydrogen storage and transport.

Selective Oxidation of Cyclohexane to Cyclohexanol/Cyclohexanone by Surface Peroxo Species on Cu-Mesoporous TiO(2)

Selective oxidation of cyclohexane to cyclohexanol/cyclohexanone (KA-oil) is an important chemical process, which is still constrained by low conversion and selectivity and high energy consumption. In this study, Cu-doped mesoporous TiO₂ (Cu-MT) has been successfully synthesized via calcinating MIL-125(Ti) doped with copper acetylacetonate, which shows high reactivity in selective oxidation of cyclohexane to KA-oil by persulfate (PS) with the desirable cyclohexane conversion of 16.8% and a selectivity of 98.0% under mild conditions and the low ratio of PS/cyclohexane of 1:1. A series of characterizations and density functional theory calculations reveal that the doped Cu(I,II) on Cu-MT is the reactive site for non-radical activation of PS with the moderate elongation of the O-O bond in PS, which then abstracts 1H (1H⁺ + 1e^-) from cyclohexane to form Cy^• and eventually KA-oil. This study gives new insight on the importance of moderately activated PS in selective oxidation of C-H.

Machine Learning Models for Predicting Molecular UV-Vis Spectra with Quantum Mechanical Properties

Accurate understanding of ultraviolet-visible (UV-vis) spectra is critical for the high-throughput synthesis of compounds for drug discovery. Experimentally determining UV-vis spectra can become expensive when dealing with a large quantity of novel compounds. This provides us an opportunity to drive computational advances in molecular property predictions using quantum mechanics and machine learning methods. In this work, we use both quantum mechanically (QM) predicted and experimentally measured UV-vis spectra as input to devise four different machine learning architectures, UVvis-SchNet, UVvis-DTNN, UVvis-Transformer, and UVvis-MPNN, and assess the performance of each method. We find that the UVvis-MPNN model outperforms the other models when using optimized 3D coordinates and QM predicted spectra as input features. This model has the highest performance for predicting UV-vis spectra with a training RMSE of 0.06 and validation RMSE of 0.08. Most importantly, our model can be used for the challenging task of predicting differences in the UV-vis spectral signatures of regioisomers.

Modeling of scavenging systems in water radiolysis with Geant4-DNA

PURPOSE

This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H₂O₂ following a ⁶⁰Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons.

METHODS

To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO₂^- and NO₃^- as scavengers and a ⁶⁰Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules.

RESULTS

The production of H₂O₂ under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature.

CONCLUSIONS

The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.

Metalla-Claisen Rearrangement in Gold-Catalyzed [4+2] Reaction: A New Elementary Reaction Suggested for Future Reaction Design

We report here computational evidence for a metalla-Claisen rearrangement (MCR) in the case of gold-catalyzed [4+2] cycloaddition reaction of yne-dienes. The [4+2] reaction starts from exo cyclopropanation, followed by MCR and reductive elimination. The cyclopropane moiety formed in the first step is crucial for a low barrier of the MCR step. In addition, the importance of an appropriate combination of the tether group and the terminal substituent on alkyne in the yne-diene substrates was studied. The mechanism of rhodium-catalyzed [4+2] reaction of yne-dienes was also investigated to see whether an MCR mechanism is involved or not. The findings and new understanding hereby reported represent an important advance in the catalysis field.

L-Shaped Heterobidentate Imidazo[1,5-a]pyridin-3-ylidene (N,C)-Ligands for Oxidant-Free Au(I) /Au(III) Catalysis

In the last decade, major advances have been made in homogeneous gold catalysis. However, Au^I /Au^III catalytic cycle remains much less explored due to the reluctance of Au^I to undergo oxidative addition and the stability of the Au^III intermediate. Herein, we report activation of aryl halides at gold(I) enabled by NHC (NHC=N-heterocyclic carbene) ligands through the development of a new class of L-shaped heterobidentate ImPy (ImPy=imidazo[1,5-a]pyridin-3-ylidene) N,C ligands that feature hemilabile character of the amino group in combination with strong σ-donation of the carbene center in a rigid conformation, imposed by the ligand architecture. Detailed characterization and control studies reveal key ligand features for Au^I /Au^III redox cycle, wherein the hemilabile nitrogen is placed at the coordinating position of a rigid framework. Given the tremendous significance of homogeneous gold catalysis, we anticipate that this ligand platform will find widespread application.

Theoretical Study on the Gas-Phase and Aqueous Interface Reaction Mechanism of Criegee Intermediates with 2-Methylglyceric Acid and the Nucleation of Products

Criegee intermediates (CIs) are important in the sink of many atmospheric substances, including alcohols, organic acids, amines, etc. In this work, the density functional theory (DFT) method was used to calculate the energy barriers for the reactions of CH3CHOO with 2-methyl glyceric acid (MGA) and to evaluate the interaction of the three functional groups of MGA. The results show that the reactions involving the COOH group of MGA are negligibly affected, and that hydrogen bonding can affect the reactions involving α-OH and β-OH groups. The water molecule has a negative effect on the reactions of the COOH group. It decreases the energy barriers of reactions involving the α-OH and β-OH groups as a catalyst. The Born-Oppenheimer molecular dynamic (BOMD) was applied to simulate the reactions of CH3CHOO with MGA at the gas-liquid interface. Water molecule plays the role of proton transfer in the reaction. Gas-phase calculations and gas-liquid interface simulations demonstrate that the reaction of CH3CHOO with the COOH group is the main pathway in the atmosphere. The molecular dynamic (MD) simulations suggest that the reaction products can form clusters in the atmosphere to participate in the formation of particles.

Color-variable dual-dyed photodynamic antimicrobial polyethylene terephthalate (PET)/cotton blended fabrics

The urgent demand for scalable, potent, color variable, and comfortable antimicrobial textiles as personal protection equipment (PPE) to help reduce infection transmission in hospitals and healthcare facilities has significantly increased since the start of the COVID-19 pandemic. Here, we explored photodynamic antimicrobial polyethylene terephthalate/cotton (TC) blended fabrics comprised of photosensitizer-conjugated cotton fibers and polyethylene terephthalate (PET) fibers dyed with disperse dyes. A small library of TC blended fabrics was constructed wherein the PET fibers were embedded with traditional disperse dyes dominating the fabric color, thereby enabling variable color expression, while the cotton fibers were covalently coupled with the photosensitizer thionine acetate as the microbicidal agent. Physical (SEM, CLSM, TGA, XPS and mechanical strength) and colorimetric (K/S and CIELab values) characterization methods were employed to investigate the resultant fabrics, and photooxidation studies with DPBF demonstrated the ability of these materials to generate reactive oxygen species (i.e., singlet oxygen) upon visible light illumination. The best results demonstrated a photodynamic inactivation of 99.985% (~ 3.82 log unit reduction, P = 0.0021) against Gram-positive S. aureus, and detection limit inactivation (99.99%, 4 log unit reduction, P ≤ 0.0001) against Gram-negative E. coli upon illumination with visible light (60 min; ~ 300 mW/cm²; λ ≥ 420 nm). Enveloped human coronavirus 229E showed a photodynamic susceptibility of ~ 99.99% inactivation after 60 min illumination (400-700 nm, 65 ± 5 mW/cm²). The presence of the disperse dyes on the fabrics showed no significant effects on the aPDI results, and furthermore, appeared to provide the photosensitizer with some measure of protection from photobleaching, thus improving the photostability of the dual-dyed fabrics. Taken together, these results suggest the feasibility of low cost, scalable and color variable thionine-conjugated TC blended fabrics as potent self-disinfecting textiles.

Positron Scattering from Pyrimidine

The positron impact cross-sections of pyrimidine molecules are reported from 1 eV to 5000 eV. These cross-sections include differential elastic, integral elastic, and direct ionisation. The elastic cross-sections are computed using the single-centre expansion scheme whereas the direct ionisation cross-sections are obtained using the binary-encounter-Bethe formula. The integral and differential cross-sections exhibit consistency with the experimental and other theoretical results. The direct ionisation cross-sections, which are reported for the first time, are compared with the experimental inelastic cross-sections (the sum of excitation and ionisation) to assess the trends in theoretically computed ionisation cross-sections and with the corresponding results for the electrons. The incoherently summed elastic and ionisation cross-sections match very well with the total cross-sections after 40 eV indicating the minimal impact of the positronium formation and electronic excitation processes. Based on this study, we recommend that the experimental data of the inelastic cross-sections reported by Palihawadana et al. be revisited.

Two Opposing Effects of Monovalent Cations on the Stability of i-Motif Structure

At acidic pH, cytosine-rich single-stranded DNA can be folded into a tetraplex structure called i-motif (iM). In recent studies, the effect of monovalent cations on the stability of iM structure has been addressed, but a consensus about the issue has not been reached yet. Thus, we investigated the effects of various factors on the stability of iM structure using fluorescence resonance energy transfer (FRET)-based analysis for three types of iM derived from human telomere sequences. We confirmed that the protonated cytosine-cytosine (C:C⁺) base pair is destabilized as the concentration of monovalent cations (Li⁺, Na⁺, K⁺) increases and that Li⁺ has the greatest tendency of destabilization. Intriguingly, monovalent cations would play an ambivalent role in iM formation by making single-stranded DNA flexible and pliant for an iM structure. In particular, we found that Li⁺ has a notably greater flexibilizing effect than Na⁺ and K⁺. All taken together, we conclude that the stability of iM structure is controlled by the subtle balance of the two counteractive effects of monovalent cations: electrostatic screening and disruption of cytosine base pairing.

Generalizable Synthesis of Highly Fluorinated Ionic Liquids

The unique chemistry of fluorocarbons (in particular, their weak intermolecular interactions and high degree of intrinsic free volume) makes them promising building blocks for ionic liquids with high gas capacities. Here, we report a generalizable method for the synthesis of fluorinated ionic liquids, which relies on the evolution of gaseous byproducts to drive product formation. This synthetic strategy overcomes solubility challenges that can hinder the synthesis of highly fluorinated ionic liquids via conventional methods and enables a systematic investigation of the effect of fluorination on ionic liquid viscosity and gas solubility.

Machine Learning-Boosted Design of Ionic Liquids for CO(2) Absorption and Experimental Verification

Efficient CO₂ capture is indispensable for achieving a carbon-neutral society while maintaining a high quality of life. Since the discovery that ionic liquids (ILs; room-temperature molten salts) can absorb CO₂, various solvents composed of molecular ions have been studied. However, it is challenging to observe the properties of each isolated ion component to control the function of ILs as they are mixtures of ions. Finding the optimal cation-anion combination for the CO₂ absorbent from their enormous chemical space had been impossible in a practical sense. This study applied electronic structure informatics to explore ILs with high CO₂ solubility from 402,114 IL candidates. The feature variables were determined by a set of cheap quantum chemistry calculations for isolated small-ion fragments, and the importance of molecular geometries and electronic states governing molecular interactions was identified via the wrapper method. As a result, it was clearly shown that the electronic states of ionic species must have essential roles in the CO₂ physisorption capacity of ILs. Considering synthetic easiness for the candidates narrowed by the machine learning model, trihexyl(tetradecyl)phosphonium perfluorooctanesulfonate was synthesized. Using a magnetic suspension balance, it was experimentally confirmed that this IL has higher CO₂ solubility than trihexyl(tetradecyl)phosphonium bis(trifluoromethanesulfonyl)amide, which is the previous best IL for CO₂ absorption.

Glycosylated BODIPY- Incorporated Pt(II) Metallacycles for Targeted and Synergistic Chemo-Photodynamic Therapy

Pt(II)-BODIPY complexes combine the chemotherapeutic activity of Pt(II) with the photocytotoxicity of BODIPYs. Additional conjugation with targeting ligands can boost the uptake by cancer cells that overexpress the corresponding receptors. We describe two Pt(II) triangles, 1 and 2, built with pyridyl BODIPYs functionalized with glucose (3) or triethylene glycol methyl ether (4), respectively. Both 1 and 2 showed higher singlet oxygen quantum yields than 3 and 4, due to the enhanced singlet-to-triplet intersystem crossing. To evaluate the targeting effect of the glycosylated derivative, in vitro experiments were performed using glucose transporter 1 (GLUT1)-positive HT29 and A549 cancer cells, and noncancerous HEK293 cells as control. Both 1 and 2 showed higher cellular uptake than 3 and 4. Specifically, 1 was selective and highly cytotoxic toward HT29 and A549 cells. The synergistic chemo- and photodynamic behavior of the metallacycles was also confirmed. Notably, 1 exhibited superior efficacy toward the cisplatin-resistant R-HepG2 cells.

Confinement effect on hydrolysis in small lipid vesicles

In living organisms most chemical reactions take place within the confines of lipid-membrane bound compartments, while confinement within the bounds of a lipid membrane is thought to be a key step in abiogenesis. In previous work we demonstrated that confinement in the aqueous cavity of a lipid vesicle affords protection against hydrolysis, a phenomenon that we term here confinement effect (C _e) and that we attributed to the interaction with the lipid membrane. Here, we show that both the size and the shape of the cavity of the vesicle modulate the C _e. We link this observation to the packing of the lipid following changes in membrane curvature, and formulate a mathematical model that relates the C _e to the radius of a spherical vesicle and the packing parameter of the lipids. These results suggest that the shape of the compartment where a molecule is located plays a major role in controlling the chemical reactivity of non-enzymatic reactions. Moreover, the mathematical treatment we propose offers a useful tool for the design of vesicles with predictable reaction rates of the confined molecules, e.g., drug delivery vesicles with confined prodrugs. The results also show that a crude form of signal transduction, devoid of complex biological machinery, can be achieved by any external stimuli that drastically changes the structure of the membrane, like the osmotic shocks used in the present work.

Ferromagnetic order controlled by the magnetic interface of LaNiO(3)/La(2/3)Ca(1/3)MnO(3) superlattices

Interface engineering in complex oxide superlattices is a growing field, enabling manipulation of the exceptional properties of these materials, and also providing access to new phases and emergent physical phenomena. Here we demonstrate how interfacial interactions can induce a complex charge and spin structure in a bulk paramagnetic material. We investigate a superlattice (SLs) consisting of paramagnetic LaNiO₃ (LNO) and highly spin-polarized ferromagnetic La_2/3Ca_1/3MnO₃ (LCMO), grown on SrTiO₃ (001) substrate. We observed emerging magnetism in LNO through an exchange bias mechanism at the interfaces in X-ray resonant magnetic reflectivity. We find non-symmetric interface induced magnetization profiles in LNO and LCMO which we relate to a periodic complex charge and spin superstructure. High resolution scanning transmission electron microscopy images reveal that the upper and lower interfaces exhibit no significant structural variations. The different long range magnetic order emerging in LNO layers demonstrates the enormous potential of interfacial reconstruction as a tool for tailored electronic properties.

Solid–Liquid Equilibrium in Co-Amorphous Systems: Experiment and Prediction

In this work, the solid–liquid equilibrium (SLE) of four binary systems combining two active pharmaceutical ingredients (APIs) capable of forming co-amorphous systems (CAMs) was investigated. The binary systems studied were naproxen-indomethacin, naproxen-ibuprofen, naproxen-probucol, and indomethacin-paracetamol. The SLE was experimentally determined by differential scanning calorimetry. The thermograms obtained revealed that all binary mixtures investigated form eutectic systems. Melting of the initial binary crystalline mixtures and subsequent quenching lead to the formation of CAM for all binary systems and most of the compositions studied. The experimentally obtained liquidus and eutectic temperatures were compared to theoretical predictions using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state and conductor-like screening model for real solvents (COSMO-RS), as implemented in the Amsterdam Modeling Suite (COSMO-RS-AMS). On the basis of the obtained results, the ability of these models to predict the phase diagrams for the investigated API–API binary systems was evaluated. Furthermore, the glass transition temperature (Tg) of naproxen (NAP), a compound with a high tendency to recrystallize, whose literature values are considerably scattered, was newly determined by measuring and modeling the Tg values of binary mixtures in which amorphous NAP was stabilized. Based on this analysis, erroneous literature values were identified.

A configuration interaction correction on top of pair coupled cluster doubles

Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N₂ and F₂ dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.

Optical trapping and light scattering in atmospheric aerosol science

Aerosol particles are ubiquitous in the atmosphere, and currently contribute a large uncertainty to climate models. Part of the endeavour to reduce this uncertainty takes the form of improving our understanding of aerosol at the microphysical level, thus enabling chemical and physical processes to be more accurately represented in larger scale models. In addition to modeling efforts, there is a need to develop new instruments and methodologies to interrogate the physicochemical properties of aerosol. This perspective presents the development, theory, and application of optical trapping, a powerful tool for single particle investigations of aerosol. After providing an overview of the role of aerosol in Earth's atmosphere and the microphysics of these particles, we present a brief history of optical trapping and a more detailed look at its application to aerosol particles. We also compare optical trapping to other single particle techniques. Understanding the interaction of light with single particles is essential for interpreting experimental measurements. In the final part of this perspective, we provide the relevant formalism for understanding both elastic and inelastic light scattering for single particles. The developments discussed here go beyond Mie theory and include both how particle and beam shape affect spectra. Throughout the entirety of this work, we highlight numerous references and examples, mostly from the last decade, of the application of optical trapping to systems that are relevant to the atmospheric aerosol.

Structural Features Promoting Photocatalytic Degradation of Contaminants of Emerging Concern: Insights into Degradation Mechanism Employing QSA/PR Modeling

Although heterogeneous photocatalysis has shown promising results in degradation of contaminants of emerging concern (CECs), the mechanistic implications related to structural diversity of chemicals, affecting oxidative (by HO•) or reductive (by O2•−) degradation pathways are still scarce. In this study, the degradation extents and rates of selected organics in the absence and presence of common scavengers for reactive oxygen species (ROS) generated during photocatalytic treatment were determined. The obtained values were then brought into correlation as K coefficients (MHO•/MO2•−), denoting the ratio of organics degraded by two occurring mechanisms: oxidation and reduction via HO• and O2•−. The compounds possessing K >> 1 favor oxidative degradation over HO•, and vice versa for reductive degradation (i.e., if K << 1 compounds undergo reductive reactions driven by O2•−). Such empirical values were brought into correlation with structural features of CECs, represented by molecular descriptors, employing a quantitative structure activity/property relationship (QSA/PR) modeling. The functional stability and predictive power of the resulting QSA/PR model was confirmed by internal and external cross-validation. The most influential descriptors were found to be the size of the molecule and presence/absence of particular molecular fragments such as C − O and C − Cl bonds; the latter favors HO•-driven reaction, while the former the reductive pathway. The developed QSA/PR models can be considered robust predictive tools for evaluating distribution between degradation mechanisms occurring in photocatalytic treatment.

Shear-induced phase transition in the aqueous solution of an imidazolium-based ionic liquid

An ionic liquid (IL) is a salt in the liquid state that consists of a cation and an anion, one of which possesses an organic component. Because of their non-volatile property, these solvents have a high recovery rate, and, hence, they are considered as environment-friendly green solvents. It is necessary to study the detailed physicochemical properties of these liquids for designing and processing techniques and find suitable operating conditions for IL-based systems. In the present work, the flow behavior of aqueous solutions of an imidazolium-based IL, 1-methyl-3-octylimidazolium chloride, is investigated, where the dynamic viscosity measurements indicate non-Newtonian shear thickening behavior in the solutions. Polarizing optical microscopy shows that the pristine samples are isotropic and transform into anisotropic after shear. These shear thickened liquid crystalline samples change into an isotropic phase upon heating, which is quantified by the differential scanning calorimetry. The small angle x-ray scattering study revealed that the pristine isotropic cubic phase of spherical micelles distort into non-spherical micelles. This has provided the detailed structural evolution of mesoscopic aggregates of the IL in an aqueous solution and the corresponding viscoelastic property of the solution.

Influence of Disorder on the Bad Metal Behavior in Polar Amalgams

The two new ternary amalgams K_1-xRb_xHg₁₁ [x = 0.472(7)] and Cs_3-xCa_xHg₂₀ [x = 0.20(3)] represent two different examples of how to create ternary compounds from binaries by statistical atom substitution. K_1-xRb_xHg₁₁ is a Vegard-type mixed crystal of the isostructural binaries KHg₁₁ and RbHg₁₁ [cubic, BaHg₁₁ structure type, space group Pm3̅m, a = 9.69143(3) Å, Rietveld refinement], whereas Cs_3-xCa_xHg₂₀ is a substitution variant of the Rb₃Hg₂₀ structure type [cubic, space group Pm3̅n, a = 10.89553(14) Å, Rietveld refinement] for which a fully substituted isostructural binary Ca phase is unknown. In K_1-xRb_xHg₁₁, the valence electron concentration (VEC) is not changed by the substitution, whereas in Cs_3-xCa_xHg₂₀, the VEC increases with the Ca content. Amalgams of electropositive metals form polar metal bonds and show "bad metal" properties. By thermal analysis, magnetic susceptibility and resistivity measurements, and density functional theory calculations of the electronic structures, we investigate the effect of the structural disorder introduced by creating mixed-atom occupation on the physical properties of the two new polar amalgam systems.

Unveiling the theoretical aspects of superelectrophilic activation in an inverse demand Diels-Alder reaction

The present computational study using B3LYP functional and 6-31+G(d) basis set has been accomplished to investigate the mechanism of the inverse demand Diels-Alder reaction between pyridyl imine and propene. The highly charged dicationic superelectrophilic diene with exceptionally low-lying LUMO makes the cycloaddition reaction with propene more favorable by significantly lowering the activation energy. The Wiberg bond indices are calculated in accordance with the formation and breaking processes of bonds. The synchronicity concept is also utilized to explain the global nature of the reaction. A potential outcome of this investigation is the utilization of propene as a C₂ building block in the industry.

Hybrid Molecular Magnets with Lanthanide- and Countercation-Mediated Interfacial Electron Transfer between Phthalocyanine and Polyoxovanadate

A series of {V₁₂}-nuclearity polyoxovanadate cages covalently functionalized with one or sandwiched by two phthalocyaninato (Pc) lanthanide (Ln) moieties via V-O-Ln bonds were prepared and fully characterized for paramagnetic Ln = Sm^III-Er^III and diamagnetic Ln = Lu^III, including Y^III. The LnPc-functionalized {V₁₂O₃₂} cages with fully oxidized vanadium centers in the ground state were isolated as (nBu₄N)₃[HV₁₂O₃₂Cl(LnPc)] and (nBu₄N)₂[HV₁₂O₃₂Cl(LnPc)₂] compounds. As corroborated by a combined experimental (EPR, DC and AC SQUID, laser photolysis transient absorption spectroscopy, and electrochemistry) and computational (DFT, MD, and model Hamiltonian approach) methods, the compounds feature intra- and intermolecular electron transfer that is responsible for a partial reduction at V(3d) centers from V^V to V^IV in the solid state and at high sample concentrations. The effects are generally Ln dependent and are clearly demonstrated for the (nBu₄N)₃[HV₁₂O₃₂Cl(LnPc)] representative with Ln = Lu^III or Dy^III. Intramolecular charge transfer takes place for Ln = Lu^III and occurs from a Pc ligand via the Ln center to the {V₁₂O₃₂} core of the same molecule, whereas for Ln = Dy^III, only intermolecular charge transfer is allowed, which is realized from Pc in one molecule to the {V₁₂O₃₂} core of another molecule usually via the nBu₄N⁺ countercation. For all Ln but Dy^III, two of these phenomena may be present in different proportions. Besides, it is demonstrated that (nBu₄N)₃[HV₁₂O₃₂Cl(DyPc)] is a field-induced single molecule magnet with a maximal relaxation time of the order 10^-3 s. The obtained results open up the way to further exploration and fine-tuning of these three modular molecular nanocomposites regarding tailoring and control of their Ln-dependent charge-separated states (induced by intramolecular transfer) and relaxation dynamics as well as of electron hopping between molecules. This should enable us to realize ultra-sensitive polyoxometalate powered quasi-superconductors, sensors, and data storage/processing materials for quantum technologies and neuromorphic computing.

3-D molecular stars with covalent axial bonding

In designing three-dimensional (3-D) molecular stars, it is very difficult to enhance the molecular rigidity through forming the covalent bonds between the axial and equatorial groups because corresponding axial groups will generally break the delocalized π bond over equatorial frameworks and thus break their star-like arrangement. In this work, exemplified by designing the 3-D stars Be₂ ©Be₅ E₅ ⁺ (E = Au, Cl, Br, I) with three delocalized σ bonds and delocalized π bond over the central Be₂ ©Be₅ moiety, we propose that the desired covalent bonding can be achieved by forming the delocalized σ bond(s) and delocalized π bond(s) simultaneously between the axial groups and equatorial framework. The covalency and rigidity of axial bonding can be demonstrated by the total Wiberg bond indices of 1.46-1.65 for axial Be atoms and ultrashort Be-Be distances of 1.834-1.841 Å, respectively. Beneficial also from the σ and π double aromaticity, these mono-cationic 3-D molecular stars are dynamically viable global energy minima with well-defined electronic structures, as reflected by wide HOMO-LUMO gaps (4.68-5.06 eV) and low electron affinities (4.70-4.82 eV), so they are the promising targets in the gas phase generation, mass-separation, and spectroscopic characterization.

Radiative Recombination Plasma Rate Coefficients for Multiply Charged Ions

Radiative recombination (RR) plasma rate coefficients are often applied to estimate electron densities and temperatures under quite different plasma conditions. Despite their frequent use, however, these rate coefficients are available only for selected (few-electron) ions and isoelectronic sequences, mainly because of the computational efforts required. To overcome this limitation, we report here a (relativistic) cascade model which helps compute fine-structure and shell-resolved as well as total RR plasma rate coefficients for many, if not most, elements of the periodic table. This model is based on Jac, the Jena Atomic Calculator, and supports studies on how the electron is captured in selected levels of the recombined ion, a relativistic (Maxwellian) electron distribution, or how the multipoles beyond the electric-dipole field in the electron-photon interaction affect the RR rate coefficients and, hence, the ionization and recombination dynamics of hot plasma. As a demonstration of this model, we compute, compare, and discuss different RR plasma rate coefficients for initially helium-like ions, with an emphasis especially on Fe24+ ions.

Direct Oral FXa Inhibitors Binding to Human Serum Albumin: Spectroscopic, Calorimetric, and Computational Studies

Direct FXa inhibitors are an important class of bioactive molecules (rivaroxaban, apixaban, edoxaban, and betrixaban) applied for thromboprophylaxis in diverse cardiovascular pathologies. The interaction of active compounds with human serum albumin (HSA), the most abundant protein in blood plasma, is a key research area and provides crucial information about drugs' pharmacokinetics and pharmacodynamic properties. This research focuses on the study of the interactions between HSA and four commercially available direct oral FXa inhibitors, applying methodologies including steady-state and time-resolved fluorescence, isothermal titration calorimetry (ITC), and molecular dynamics. The HSA complexation of FXa inhibitors was found to occur via static quenching, and the complex formation in the ground states affects the fluorescence of HSA, with a moderate binding constant of 10⁴ M^-1. However, the ITC studies reported significantly different binding constants (10³ M^-1) compared with the results obtained through spectrophotometric methods. The suspected binding mode is supported by molecular dynamics simulations, where the predominant interactions were hydrogen bonds and hydrophobic interactions (mainly π-π stacking interactions between the phenyl ring of FXa inhibitors and the indole moiety of Trp214). Finally, the possible implications of the obtained results regarding pathologies such as hypoalbuminemia are briefly discussed.

Insights into the Absorption of Hydrocarbon Gases in Phosphorus-Containing Ionic Liquids

The solubility of ethane, ethylene, propane, and propylene was measured in two phosphorus-containing ionic liquids, trihexyltetradecylphosphonium bis(2,4,4-trimethylpentyl)phosphinate, [P_6,6,6,14][DiOP], and 1-butyl-3-methylimidazolium dimethylphosphate, [C₄C₁Im][DMP], using an isochoric saturation method. The ionic liquid [C₄C₁Im][DMP] absorbed between 1 and 20 molecules of gas per 1000 ion pairs, at 313 K and 0.1 MPa, while [P_6,6,6,14][DiOP] absorbed up to 169 molecules of propane per 1000 ion pairs under the same conditions. [C₄C₁Im][DMP] had a higher capacity to absorb olefins than paraffins, while the opposite was true for [P_6,6,6,14][DiOP], with the former being slightly more selective than the later. From the analysis of the thermodynamic properties of solvation, we concluded that in both ionic liquids and for all of the studied gases the solvation is ruled by the entropy, even if its contribution is unfavorable. These results, together with density measurements, 2D NMR studies, and self-diffusion coefficients suggest that the gases' solubility is ruled mostly by nonspecific interactions with the ionic liquids and that the looser ion packing in [P_6,6,6,14][DiOP] makes it easier to accommodate the gases compared to [C₄C₁Im][DMP].

Rigorous Negative Ion Binding Energies in Low-Energy Electron Elastic Collisions with Heavy Multi-Electron Atoms and Fullerene Molecules: Validation of Electron Affinities

Dramatically sharp resonances manifesting stable negative ion formation characterize Regge pole-calculated low-energy electron elastic total cross sections (TCSs) of heavy multi-electron systems. The novelty of the Regge pole analysis is in the extraction of rigorous and unambiguous negative ion binding energies (BEs), corresponding to the measured electron affinities (EAs) of the investigated multi-electron systems. The measured EAs have engendered the crucial question: is the EA of multi-electron atoms and fullerene molecules identified with the BE of the attached electron in the ground, metastable or excited state of the formed negative ion during a collision? Inconsistencies in the meaning of the measured EAs are elucidated and new EA values for Bk, Cf, Fm, and Lr are presented.

Exploring the Effect of Different Reactivity Promoters on the Oxidation of Ammonia in a Jet-Stirred Reactor

The low reactivity of ammonia (NH₃) is the main barrier to applying neat NH₃ as fuel in technical applications, such as internal combustion engines and gas turbines. Introducing combustion promoters as additives in NH₃-based fuel can be a feasible solution. In this work, the oxidation of ammonia by adding different reactivity promoters, i.e., hydrogen (H₂), methane (CH₄), and methanol (CH₃OH), was investigated in a jet-stirred reactor (JSR) at temperatures between 700 and 1200 K and at a pressure of 1 bar. The effect of ozone (O₃) was also studied, starting from an extremely low temperature (450 K). Species mole fraction profiles as a function of the temperature were measured by molecular-beam mass spectrometry (MBMS). With the help of the promoters, NH₃ consumption can be triggered at lower temperatures than in the neat NH₃ case. CH₃OH has the most prominent effect on enhancing the reactivity, followed by H₂ and CH₄. Furthermore, two-stage NH₃ consumption was observed in NH₃/CH₃OH blends, whereas no such phenomenon was found by adding H₂ or CH₄. The mechanism constructed in this work can reasonably reproduce the promoting effect of the additives on NH₃ oxidation. The cyanide chemistry is validated by the measurement of HCN and HNCO. The reaction CH₂O + NH₂ ⇄ HCO + NH₃ is responsible for the underestimation of CH₂O in NH₃/CH₄ fuel blends. The discrepancies observed in the modeling of NH₃ fuel blends are mainly due to the deviations in the neat NH₃ case. The total rate coefficient and the branching ratio of NH₂ + HO₂ are still controversial. The high branching fraction of the chain-propagating channel NH₂ + HO₂ ⇄ H₂NO + OH improves the model performance under low-pressure JSR conditions for neat NH₃ but overestimates the reactivity for NH₃ fuel blends. Based on this mechanism, the reaction pathway and rate of production analyses were conducted. The HONO-related reaction routine was found to be activated uniquely by adding CH₃OH, which enhances the reactivity most significantly. It was observed from the experiment that adding ozone to the oxidant can effectively initiate NH₃ consumption at temperatures below 450 K but unexpectedly inhibit the NH₃ consumption at temperatures higher than 900 K. The preliminary mechanism reveals that adding the elementary reactions between NH₃-related species and O₃ is effective for improving the model performance, but their rate coefficients have to be refined.

Quantum Chemical Modeling of Propellant Degradation

An ab initio quantum chemical approach for the modeling of propellant degradation is presented. Using state-of-the-art bonding analysis techniques and composite methods, a series of potential degradation reactions are devised for a sample hydroxyl-terminated-polybutadiene (HTPB) type solid fuel. By applying thermochemical procedures and isodesmic reactions, accurate thermochemical quantities are obtained using a modified G3 composite method based on the resolution of the identity. The calculated heats of formation for the different structures produced presents an ∼2 kcal/mol average error when compared against experimental values.

R2022: A DFT/MRCI Ansatz with Improved Performance for Double Excitations

A reformulation of the combined density functional theory and multireference configuration interaction method (DFT/MRCI) is presented. Expressions for ab initio matrix elements are used to derive correction terms for a new effective Hamiltonian. On the example of diatomic carbon, the correction terms are derived, focusing on the doubly excited ¹Δ_g state, which was problematic in previous formulations of the method, as were double excitations in general. The derivation shows that a splitting of the parameters for intra- and interorbital interactions is necessary for a concise description of the underlying physics. Results for ¹L_a and ¹L_b states in polyacenes and ¹A_u and ¹A_g states in mini-β-carotenoids suggest that the presented formulation is superior to former effective Hamiltonians. Furthermore, statistical analysis reveals that all the benefits of the previous DFT/MRCI Hamiltonians are retained. Consequently, the here presented formulation should be considered as the new standard for DFT/MRCI calculations.

EUV Beam–Foil Spectra of Germanium and a Blind-Spot Problem in Spectroscopy

Beam–foil extreme-ultraviolet survey spectra of Ge (Z=32) are presented. The data have been garnered at the performance limit of the heavy-ion accelerator available, with a correspondingly limited statistical and calibrational reliability. However, the Ge spectra have been recorded at various delays after excitation, and this technique points to a possible blind spot in some other spectroscopic techniques, and thus in the literature coverage. A similarly patchy coverage can be noted in various atomic structure computations. The experimental and theoretical gaps seem to be correlated.