Anharmonic Vibrational Frequencies of Water Borane and Associated Molecules

Water borane (BH3OH2) and borinic acid (BH2OH) have been proposed as intermediates along the pathway of hydrogen generation from simple reactants: water and borane. However, the vibrational spectra for neither water borane nor borinic acid has been investigaged experimentally due to the difficulty of isolating them in the gas phase, making accurate quantum chemical predictions for such properties the most viable means of their determination. This work presents theoretical predictions of the full rotational and fundamental vibrational spectra of these two potentially application-rich molecules using quartic force fields at the CCSD(T)-F12b/cc-pCVTZ-F12 level with additional corrections included for the effects of scalar relativity. This computational scheme is further benchmarked against the available gas-phase experimental data for the related borane and HBO molecules. The differences are found to be within 3 cm−1 for the fundamental vibrational frequencies and as close as 15 MHz in the B0 and C0 principal rotational constants. Both BH2OH and BH3OH2 have multiple vibrational modes with intensities greater than 100 km mol−1, namely ν2 and ν4 in BH2OH, and ν1, ν3, ν4, ν9, and ν13 in BH3OH2. Finally, BH3OH2 has a large dipole moment of 4.24 D, which should enable it to be observable by rotational spectroscopy, as well.     

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.     

Metastable Kitaev Magnets

Nearly two decades ago, Alexei Kitaev proposed a model for spin-1/2 particles with bond-directional interactions on a two-dimensional honeycomb lattice which had the potential to host a quantum spin-liquid ground state. This work initiated numerous investigations to design and synthesize materials that would physically realize the Kitaev Hamiltonian. The first generation of such materials, such as Na2IrO3, α-Li2IrO3, and α-RuCl3, revealed the presence of non-Kitaev interactions such as the Heisenberg and off-diagonal exchange. Both physical pressure and chemical doping were used to tune the relative strength of the Kitaev and competing interactions; however, little progress was made towards achieving a purely Kitaev system. Here, we review the recent breakthrough in modifying Kitaev magnets via topochemical methods that has led to the second generation of Kitaev materials. We show how structural modifications due to the topotactic exchange reactions can alter the magnetic interactions in favor of a quantum spin-liquid phase.     

An Exchange Mechanism for the Magnetic Behavior of Er3+ Complexes

We study the magnetic properties of the erbium based compounds, Na9[Er(W5O18)2] and [(Pc)Er{Pc{N(C4H9)2}8}]·/−, in the framework of an effective spin exchange model involving delocalized electrons occupying molecular orbitals. The calculations successfully reproduce the experimental data available in the literature for the magnetic spectrum, magnetization and molar susceptibility in dc and ac fields. Owing to their similar molecular geometry, the compounds’ magnetic behaviors are interpreted in terms of the same set of active orbitals and thus the same effective spin coupling scheme. For all three complexes, the model predicts a prompt change in the ground state from a Kramer’s doublet at zero fields to a fully polarized quartet one brought about by the action of an external magnetic field without Zeeman splitting. This alteration is attributed to the enhancement of the effect of orbital interactions over the spin exchange as the magnitude of the external magnetic field increases.     

Fluorescent Liquid Tetrazines

Tetrazines with branched alkoxy substituents are liquids at ambient temperature that despite the high chromophore density retain the bright orange fluorescence that is characteristic of this exceptional fluorophore. Here, we study the photophysical properties of a series of alkoxy-tetrazines in solution and as neat liquids. We also correlate the size of the alkoxy substituents with the viscosity of the liquids. We show using time-resolved spectroscopy that intersystem crossing is an important decay pathway competing with fluorescence, and that its rate is higher for 3,6-dialkoxy derivatives than for 3-chloro-6-alkoxytetrazines, explaining the higher fluorescence quantum yields for the latter. Quantum chemical calculations suggest that the difference in rate is due to the activation energy required to distort the tetrazine core such that the nπ*S1 and the higher-lying ππ*T2 states cross, at which point the spin-orbit coupling exceeding 10 cm−1 allows for efficient intersystem crossing to occur. Femtosecond time-resolved anisotropy studies in solution allow us to measure a positive relationship between the alkoxy chain lengths and their rotational correlation times, and studies in the neat liquids show a fast decay of the anisotropy consistent with fast exciton migration in the neat liquid films.     

Galvano-Fenton Engineering Solution with Spontaneous Catalyst’s Generation from Waste: Experimental Efficiency,Parametric Analysis and Modeling Interpretation Applied to a Clean Technology for Dyes Degradation in Water

In this paper, the degradation of the diazo dye naphthol blue black (NBB) using the Galvano-Fenton process is studied experimentally and numerically. The simulations are carried out based on the anodic, cathodic, and 34 elementary reactions evolving in the electrolyte, in addition to the oxidative attack of NBB by HO• at a constant rate of 3.35×107 mol−1·m3·s−1 during the initiation stage of the chain reactions. The selection of the operating conditions including the pH of the electrolyte, the stirring speed, and the electrodes disposition is performed by assessing the kinetics of NBB degradation; these parameters are set to 3, 350 rpm and a parallel disposition with a 3 cm inter-electrode distance, respectively. The kinetics of Fe(III) in the electrolyte were monitored using the principles of Fricke dosimetry and simulated numerically. The model showed more than a 96% correlation with the experimental results in both the blank test and the presence of the dye. The effects of H2O2 and NBB concentrations on the degradation of the dye were examined jointly with the evolution of the simulated H2O2, Fe2+, and HO• concentrations in the electrolyte. The model demonstrated a good correlation with the experimental results in terms of the initial degradation rates, with correlation coefficients exceeding 98%.     

Operando Raman Shift Replaces Current in Electrochemical Analysis of Li-ion Batteries: A Comparative Study

Li-rich and catalytically active γ-LixV2O5 (x = 1.48) was investigated as a cathode for its heterogeneous charge transfer kinetics. Using a specially designed two-electrode system lithium half cell, Butler–Volmer analysis was performed, and Raman spectra were acquired in 18 mV intervals. A direct correlation was observed between the Raman shift of the active modes Ag,Bg, Au, and Bu, and the development of the Faraday current at the working electrode. The Raman intensity and the Raman shift were implemented to replace the current in a Tafel plot used for the analysis of Butler–Volmer kinetics. Striking similarities in the charge transfer proportionality constants α were found for current and Raman-based analysis. The potential of this new method of Raman-aided electrochemical detection at the diffraction limit is discussed.     

Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt

In the hydrochloride of a pyrazolyl-substituted acetylacetone, the chloride anion is hydrogen-bonded to the protonated pyrazolyl moiety. Equimolar co-crystallization with tetrafluorodiiodobenzene (TFDIB) leads to a supramolecular aggregate in which TFDIB is situated on a crystallographic center of inversion. The iodine atom in the asymmetric unit acts as halogen bond donor, and the chloride acceptor approaches the σ-hole of this TFDIB iodine subtending an almost linear halogen bond, with Cl···I = 3.1653(11) Å and Cl···I–C = 179.32(6)°. This contact is roughly orthogonal to the N–H···Cl hydrogen bond. An analysis of the electron density according to Bader’s Quantum Theory of Atoms in Molecules confirms bond critical points (bcps) for both short contacts, with ρbcp = 0.129 for the halogen and 0.321 eÅ−3 for the hydrogen bond. Our halogen-bonded adduct represents the prototype for a future class of co-crystals with tunable electron density distribution about the σ-hole contact.     

The Relationship between the IC50 Values and the Apparent Inhibition Constant in the Study of Inhibitors of Tyrosinase Diphenolase Activity Helps Confirm the Mechanism of Inhibition

Tyrosinase is the enzyme involved in melanization and is also responsible for the browning of fruits and vegetables. Control of its activity can be carried out using inhibitors, which is interesting in terms of quantitatively understanding the action of these regulators. In the study of the inhibition of the diphenolase activity of tyrosinase, it is intriguing to know the strength and type of inhibition. The strength is indicated by the value of the inhibition constant(s), and the type can be, in a first approximation: competitive, non-competitive, uncompetitive and mixed. In this work, it is proposed to calculate the degree of inhibition (iD), varying the concentration of inhibitor to a fixed concentration of substrate, L-dopa (D). The non-linear regression adjustment of iD with respect to the initial inhibitor concentration [I]0 allows for the calculation of the inhibitor concentration necessary to inhibit the activity by 50%, at a given substrate concentration (IC₅₀), thus avoiding making interpolations between different values of iD. The analytical expression of the IC₅₀, for the different types of inhibition, are related to the apparent inhibition constant (KIapp). Therefore, this parameter can be used: (a) To classify a series of inhibitors of an enzyme by their power. Determining these values at a fixed substrate concentration, the lower IC₅₀, the more potent the inhibitor. (b) Checking an inhibitor for which the type and the inhibition constant have been determined (using the usual methods), must confirm the IC₅₀ value according to the corresponding analytical expression. (c) The type and strength of an inhibitor can be analysed from the study of the variation in iD and IC₅₀ with substrate concentration. The dependence of IC₅₀ on the substrate concentration allows us to distinguish between non-competitive inhibition (iD does not depend on [D]0) and the rest. In the case of competitive inhibition, this dependence of iD on [D]0 leads to an ambiguity between competitive inhibition and type 1 mixed inhibition. This is solved by adjusting the data to the possible equations; in the case of a competitive inhibitor, the calculation of KI1app is carried out from the IC₅₀ expression. The same occurs with uncompetitive inhibition and type 2 mixed inhibition. The representation of iD vs. n, with n=[D]0/KmD, allows us to distinguish between them. A hyperbolic iD vs. n representation that passes through the origin of coordinates is a characteristic of uncompetitive inhibition; the calculation of KI2app is immediate from the IC₅₀ value. In the case of mixed inhibitors, the values of the apparent inhibition constant of meta-tyrosinase (Em) and oxy-tyrosinase (Eox), KI1app and the apparent inhibition constant of metatyrosinase/Dopa complexes (EmD) and oxytyrosinase/Dopa (EoxD), KI2app are obtained from the dependence of iD vs. n, and the results obtained must comply with the IC₅₀ value.     

Millimeter-Wave Spectroscopy of Methylfuran Isomers: Local vs. Global Treatments of the Internal Rotation

Methylfurans are methylated aromatic heterocyclic volatile organic compounds and primary or secondary pollutants in the atmosphere due to their capability to form secondary organic aerosols in presence of atmospheric oxidants. There is therefore a significant interest to monitor these molecules in the gas phase. High resolution spectroscopic studies of methylated furan compounds are generally limited to pure rotational spectroscopy in the vibrational ground state. This lack of results might be explained by the difficulties arisen from the internal rotation of the methyl group inducing non-trivial patterns in the rotational spectra. In this study, we discuss the benefits to assign the mm-wave rotational-torsional spectra of methylfuran with the global approach of the BELGI-Cs code compared to local approaches such as XIAM and ERHAM. The global approach reproduces the observed rotational lines of 2-methylfuran and 3-methylfuran in the mm-wave region at the experimental accuracy for the ground vt=0 and the first torsional vt=1 states with a unique set of molecular parameters. In addition, the V3 and V6 parameters describing the internal rotation potential barrier may be determined with a high degree of accuracy with the global approach. Finally, a discussion with other heterocyclic compounds enables the study of the influence of the electronic environment on the hindered rotation of the methyl group.     

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.     

Exogenous Glutathione Protects IPEC-J2 Cells against Oxidative Stress through a Mitochondrial Mechanism

The accumulation of reactive oxygen species (ROS) triggers oxidative stress in cells by oxidizing and modifying various cellular components, preventing them from performing their inherent functions, ultimately leading to apoptosis and autophagy. Glutathione (GSH) is a ubiquitous intracellular peptide with multiple functions. In this study, a hydrogen peroxide (H2O2)-induced oxidative damage model in IPEC-J2 cells was used to investigate the cellular protection mechanism of exogenous GSH against oxidative stress. The results showed that GSH supplement improved the cell viability reduced by H2O2-induced oxidative damage model in IPEC-J2 cells in a dose-dependent manner. Moreover, supplement with GSH also attenuated the H2O2-induced MMP loss, and effectively decreased the H2O2-induced mitochondrial dysfunction by increasing the content of mtDNA and upregulating the expression TFAM. Exogenous GSH treatment significantly decreased the ROS and MDA levels, improved SOD activity in H2O2-treated cells and reduced H2O2-induced early apoptosis in IPEC-J2 cells. This study showed that exogenous GSH can protect IPEC-J2 cells against apoptosis induced by oxidative stress through mitochondrial mechanisms.     

Ca+ Ions Solvated in Helium Clusters

We present a combined experimental and theoretical investigation on Ca+ ions in helium droplets, HeNCa+. The clusters have been formed in the laboratory by means of electron-impact ionization of Ca-doped helium nanodroplets. Energies and structures of such complexes have been computed using various approaches such as path integral Monte Carlo, diffusion Monte Carlo and basin-hopping methods. The potential energy functions employed in these calculations consist of analytical expressions following an improved Lennard-Jones formula whose parameters are fine-tuned by exploiting ab initio estimations. Ion yields of HeNCa+ -obtained via high-resolution mass spectrometry- generally decrease with N with a more pronounced drop between N=17 and N=25, the computed quantum HeNCa+ evaporation energies resembling this behavior. The analysis of the energies and structures reveals that covering Ca+ with 17 He atoms leads to a cluster with one of the smallest energies per atom. As new atoms are added, they continue to fill the first shell at the expense of reducing its stability, until N=25, which corresponds to the maximum number of atoms in that shell. Behavior of the evaporation energies and radial densities suggests liquid-like cluster structures.     

Influence of Liquid Crystallinity and Mechanical Deformation on the Molecular Relaxations of an Auxetic Liquid Crystal Elastomer

Liquid Crystal Elastomers (LCEs) combine the anisotropic ordering of liquid crystals with the elastic properties of elastomers, providing unique physical properties, such as stimuli responsiveness and a recently discovered molecular auxetic response. Here, we determine how the molecular relaxation dynamics in an acrylate LCE are affected by its phase using broadband dielectric relaxation spectroscopy, calorimetry and rheology. Our LCE is an excellent model system since it exhibits a molecular auxetic response in its nematic state, and chemically identical nematic or isotropic samples can be prepared by cross-linking. We find that the glass transition temperatures (Tg) and dynamic fragilities are similar in both phases, and the T-dependence of the α relaxation shows a crossover at the same T* for both phases. However, for T>T*, the behavior becomes Arrhenius for the nematic LCE, but only more Arrhenius-like for the isotropic sample. We provide evidence that the latter behavior is related to the existence of pre-transitional nematic fluctuations in the isotropic LCE, which are locked in by polymerization. The role of applied strain on the relaxation dynamics and mechanical response of the LCE is investigated; this is particularly important since the molecular auxetic response is linked to a mechanical Fréedericksz transition that is not fully understood. We demonstrate that the complex Young’s modulus and the α relaxation time remain relatively unchanged for small deformations, whereas for strains for which the auxetic response is achieved, significant increases are observed. We suggest that the observed molecular auxetic response is coupled to the strain-induced out-of-plane rotation of the mesogen units, in turn driven by the increasing constraints on polymer configurations, as reflected in increasing elastic moduli and α relaxation times; this is consistent with our recent results showing that the auxetic response coincides with the emergence of biaxial order.     

TiCl4 Dissolved in Ionic Liquid Mixtures from Аb Initio Molecular Dynamics Simulations

To gain a deeper understanding of the TiCl4 solvation effects in multi-component ionic liquids, we performed ab initio molecular dynamics simulations of 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]−, chloride [Cl]− both with and without water and titanium tetrachloride TiCl4. Complex interactions between cations and anions are observed in all investigated systems. By further addition of water and TiCl4 this complex interaction network is extended. Observations of the radial distribution functions and number integrals show that water and TiCl4 not only compete with each other to interact mainly with [Cl]−, which strongly influences the cation-[BF4]− interaction, but also interact with each other, which leads to the fact that in certain systems the cation-anion interaction is enhanced. Further investigations of the Voronoi polyhedra analysis have demonstrated that water has a greater impact on the nanosegregated system than TiCl4 which is also due to the fact of the shear amount of water relative to all other components and its higher mobility compared to TiCl4. Overall, the polar network of the IL mixture collapses by including water and TiCl4. In the case of [Cl]− chloride enters the water continuum, while [BF4]− remains largely unaffected, which deeply affects the interaction of the ionic liquid (IL) network.     

Photochemistry of Thymine in Protic Polar Nanomeric Droplets Using Electrostatic Embeding TD-DFT/MM

Thymine photochemistry is important for understanding DNA photodamage. In the gas phase, thymine undergoes a fast non-radiative decay from S 2 to S 1 . In the S 1 state, it gets trapped for several picoseconds until returning to the ground-state S 0 . Here, we explore the electrostatic effects of nanomeric droplets of methanol and water on the excited states of thymine. For this purpose, we develop and implement an electrostatic embedding TD-DFT/MM method based on a QM/MM coupling defined through electrostatic potential fitting charges. We show that both in methanol and water, the mechanism is similar to the gas phase. The solvent molecules participate in defining the branching plane of S 0 /S 1 intersection and have a negligible effect on the S 1 /S 2 intersection. Despite the wrong topology of the ground/excited state intersections, electrostatic embedding TD-DFT/MM allows for a fast exploration of the potential energy surfaces and a qualitative picture of the photophysics of thymine in solvent droplets.     

Optical Dynamic Nuclear Polarization of 13C Spins in Diamond at a Low Field with Multi-Tone Microwave Irradiation

Majority of dynamic nuclear polarization (DNP) experiments have been requiring helium cryogenics and strong magnetic fields for a high degree of nuclear polarization. In this work, we instead demonstrate an optical hyperpolarization of naturally abundant 13C nuclei in a diamond crystal at a low magnetic field and the room temperature. It exploits continuous laser irradiation for polarizing electronic spins of nitrogen vacancy centers and microwave irradiation for transferring the electronic polarization to 13C nuclear spins. We have studied the dependence of 13C polarization on laser and microwave powers. For the first time, a triplet structure corresponding to the 14N hyperfine splitting has been observed in the 13C polarization spectrum. By simultaneously exciting three microwave frequencies at the peaks of the triplet, we have achieved 13C bulk polarization of 0.113 %, leading to an enhancement of 90,000 over the thermal polarization at 17.6 mT. We believe that the multi-tone irradiation can be extended to further enhance the 13C polarization at a low magnetic field.     

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.     

A Platform Approach to Protein Encapsulates with Controllable Surface Chemistry

The encapsulation of proteins into core-shell structures is a widely utilised strategy for controlling protein stability, delivery and release. Despite the recognised utility of these microstructures, however, core-shell fabrication routes are often too costly or poorly scalable to allow for industrial translation. Furthermore, many scalable routes rely upon emulsion-techniques implicating denaturing or environmentally harmful organic solvents. Herein, we investigate core-shell protein encapsulation through single-feed, aqueous spray drying: a cheap, industrially ubiquitous particle-formation technology in the absence of organic solvents. We show that an excipient’s preference for the surface of the spray dried particle is well-predicted by its hydrodynamic diameter (Dh) under relevant feed buffer conditions (pH and ionic strength) and that the predictive power of Dh is improved when measured at the spray dryer outlet temperature compared to room temperature (R2 = 0.64 vs. 0.59). Lastly, we leverage these findings to propose an adaptable design framework for fabricating core-shell protein encapsulates by single-feed aqueous spray drying.