RESUMEN
We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.
RESUMEN
The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.
RESUMEN
One of the most widely used and computationally efficient models that accounts for London dispersion interactions within density functional theory (DFT) is the D3 dispersion correction model. In this work, we demonstrate that this model can induce the appearance of unphysical minima on the potential energy surface (PES) when the coordination number of atoms changes. Optimizing to these artifactual structures can lead to significant errors in determining the interaction energy between two molecules and in estimating the thermodynamic properties of the system. In several specific examples, such as Kuratowski-type H2-NiKur and H2-PdKur clusters, these local minima exhibited extremely high PES curvature, resulting in incorrect estimations of harmonic frequencies and significant overestimations of zero-point energy and enthalpy values. Although such erroneous behavior of the D3 model is relatively rare, it can occur across a wide range of chemical species, including molecules like the [Li(C6H6)]+ complex and the dispiro(acridan)-substituted pyracene (DSAP) molecule. Our analysis reveals that the root of the problem lies in the definition of the AB atomic-pair dependent C6AB coefficients in the D3 model. To address this issue, we propose a reparameterization of the D3 model by introducing a modified C6AB functional form that now depends on the specific pair of considered atoms. This new model, termed D3-Smooth (or D3S for short), is designed to smooth out the PES associated with the dispersion correction. By doing so, we demonstrate that D3S eliminates unphysical local minima while maintaining the quite satisfactory accuracy of the parent D3 method in interaction energy benchmark sets. For example, the RMS difference between using the D3(BJ) correction to B3LYP and the D3S(BJ) correction across the large MGCDB84 data set of nearly 5000 data points is only 0.12 kJ/mol. Similar results are obtained for every other D3-corrected functional tested. Consistent with this result, no significant improvement could be obtained for the B3LYP-D3S(0) correction by reoptimizing the damping function.
RESUMEN
X-ray Transient Absorption Spectroscopy (XTAS) is a powerful probe for ultrafast molecular dynamics. The evolution of XTAS signal is controlled by the shapes of potential energy surfaces of the associated core-excited states, which are difficult to directly measure. Here, we study the vibrational dynamics of Raman activated CCl4 with XTAS targeting the C 1s and Cl 2p electrons. The totally symmetric stretching mode leads to concerted elongation or contraction in bond lengths, which in turn induce an experimentally measurable red or blue shift in the X-ray absorption energies associated with inner-shell electron excitations to the valence antibonding levels. The ratios between slopes of different core-excited potential energy surfaces (CEPESs) thereby extracted agree very well with Restricted Open-Shell Kohn-Sham calculations. The other, asymmetric, modes do not measurably contribute to the XTAS signal. The results highlight the ability of XTAS to reveal coherent nuclear dynamics involving < 0.01 Å atomic displacements and also provide direct measurement of forces on CEPESs.
RESUMEN
The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g., triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using ab initio electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. In this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of 1013 s-1 for each dye molecule.
RESUMEN
In order to alleviate the computational burden associated with superlinear compute scalings with molecular size in electron correlation methods, researchers have developed local correlation methods that wisely treat relatively small contributions as zeros but still yield accurate energy approximation. Such local correlation techniques can also be combined with parallel computing resources to obtain further efficiency and scalability. This work focuses on the distributed memory parallel implementation of a local correlation method for second order MoÌ·ller-Plesset (MP2) theory. This method also only has a single threshold to control the dropping of terms and accuracy of different computing kernels in the algorithm. The process partitioning strategy and distributed parallel implementation with the message passing interface (MPI) are discussed. In particular, the algorithm relies on a fixed sparsity pattern matrix multiplication and a corresponding distributed conjugate gradient solver, which exhibits almost linear scaling in both strong and weak scaling analyses. Numerical experiments on a range of molecules, including linear chains and molecules with 2 and 3-dimensional characters, are reported. For example, with only 32 MPI ranks, this MP2 implementation can calculate the correlation energy of vancomycin in def2-TZVP basis within 0.003% accuracy (10-6.5 threshold) in half an hour, where the same problem is unfeasible to solve with sequential or pure shared memory implementations.
RESUMEN
In periodic systems, the Hartree-Fock (HF) exchange energy exhibits the slowest convergence of all HF energy components as the system size approaches the thermodynamic limit. We demonstrate that the recently proposed staggered mesh method for Fock exchange energy [Xing, Li, and Lin, Math. Comp., 2024], which is specifically designed to sidestep certain singularities in exchange energy evaluation, can expedite the finite-size convergence rate for the exact exchange energy across a range of insulators and semiconductors when compared to the regular and truncated Coulomb methods. This remains true even for two computationally cheaper versions of this new method, which we call non-SCF and split-SCF staggered mesh. Additionally, a sequence of numerical tests on simple solids showcases the staggered mesh method's ability to improve convergence toward the thermodynamic limit for band gaps, bulk moduli, equilibrium lattice dimensions, energies, and phonon force constants.
RESUMEN
We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy.
RESUMEN
Porous solids can accommodate and release molecular hydrogen readily, making them attractive for minimizing the energy requirements for hydrogen storage relative to physical storage systems. However, H2 adsorption enthalpies in such materials are generally weak (-3 to -7 kJ/mol), lowering capacities at ambient temperature. Metal-organic frameworks with well-defined structures and synthetic modularity could allow for tuning adsorbent-H2 interactions for ambient-temperature storage. Recently, Cu2.2Zn2.8Cl1.8(btdd)3 (H2btdd = bis(1H-1,2,3-triazolo-[4,5-b],[4',5'-i])dibenzo[1,4]dioxin; CuI-MFU-4l) was reported to show a large H2 adsorption enthalpy of -32 kJ/mol owing to π-backbonding from CuI to H2, exceeding the optimal binding strength for ambient-temperature storage (-15 to -25 kJ/mol). Toward realizing optimal H2 binding, we sought to modulate the π-backbonding interactions by tuning the pyramidal geometry of the trigonal CuI sites. A series of isostructural frameworks, Cu2.7M2.3X1.3(btdd)3 (M = Mn, Cd; X = Cl, I; CuIM-MFU-4l), was synthesized through postsynthetic modification of the corresponding materials M5X4(btdd)3 (M = Mn, Cd; X = CH3CO2, I). This strategy adjusts the H2 adsorption enthalpy at the CuI sites according to the ionic radius of the central metal ion of the pentanuclear cluster node, leading to -33 kJ/mol for M = ZnII (0.74 Å), -27 kJ/mol for M = MnII (0.83 Å), and -23 kJ/mol for M = CdII (0.95 Å). Thus, CuICd-MFU-4l provides a second, more stable example of optimal H2 binding energy for ambient-temperature storage among reported metal-organic frameworks. Structural, computational, and spectroscopic studies indicate that a larger central metal planarizes trigonal CuI sites, weakening the π-backbonding to H2.
RESUMEN
Excitation of iron pentacarbonyl [Fe(CO)5], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO)5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental-theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M2,3-edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO)5. We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO)4 along the lowest singlet and triplet surfaces to form Fe(CO)3. Calculations of transient spectra in both singlet and triplet states as well as spin-orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO)4. Thus, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO)5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.
RESUMEN
We investigate the early stages of cesium lead bromide perovskite formation through absorption spectroscopy of stopped-flow reactions, high-throughput mapping, and direct synthesis and titration of potential precursor species. Calorimetric and spectroscopic measurements of lead bromide complex titrations combined with theoretical calculations suggest that bromide complexes with higher coordination numbers than previously considered for nonpolar systems can better explain observed behaviors. Synthesis mapping of binary lead halides reveals multiple lead bromide species with absorption peaks higher than 300 nm, including a previously observed species with a peak at 313 nm and two species with peaks at 345 and 370 nm that also appear as reaction intermediates during the formation of lead bromide perovskites. Based on theoretical calculations of excitonic energies that match within 50 meV, we give a preliminary assignment of these species as two-dimensional magic-sized clusters with side lengths of 2, 3, and 4 unit cells. Kinetic measurements of the conversion of benzoyl bromide precursor are connected to stopped-flow measurements of product formation and demonstrate that the formation of complexes and magic-sized clusters (i.e., nucleation) is controlled by precursor decomposition, whereas the growth rate of 2D and 3D perovskites is significantly slower.
RESUMEN
X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core-electron binding energy (CEBE) analysis. We extend the development of real-space pseudopotential methods to calculating 1s, 2s, and 2p3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn-Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.
RESUMEN
In this article, we introduce the occupied-virtual orbitals for chemical valence (OVOCV). The OVOCVs can replace or complement the closely related idea of the natural orbitals for chemical valence (NOCV). The input is a difference density matrix connecting any initial single determinant to any final determinant, at a given molecular geometry, and a given one-particle basis. This arises in problems such as orbital rearrangement or charge transfer (CT) in energy decomposition analysis (EDA). The OVOCVs block-diagonalize the density difference operator into 2 × 2 blocks, which are spanned by one level that is filled in the initial state (the occupied OVOCV) and one that is empty (the virtual OVOCV). By contrast, the NOCVs fully diagonalize the density difference matrix and therefore are orbitals with mixed occupied-virtual character. Use of the OVOCVs makes it much easier to identify the donor and acceptor orbitals. We also introduce two different types of EDA methods with the OVOCVs and, most importantly, a charge decomposition analysis method that fixes the unreasonably large CT amount obtained directly from NOCV analysis. The square of the CT amount associated with each NOCV pair emerges as the appropriate value from the OVOCV analysis. When connecting the same initial and final states, this value is identical to the CT amount obtained from the independent absolutely localized molecular orbital (ALMO) complementary occupied-virtual orbital pair (COVP) analysis. The total, summed over all pairs, is also exactly the same as the independently suggested excitation number, as proved herein. Several examples are presented to compare NOCVs and OVOCVs: stretched H2+, a strong halogen bond between tetramethylthiourea and iodine, coordination of ethene in Zeise's salt, and binding in the Cp3La···C≡NCy complex.
RESUMEN
BACKGROUND: Hyaluronan (HA) is an extracellular glycosaminoglycan polysaccharide with widespread roles throughout development and in healthy and neoplastic tissues. In pluripotent stem cell culture it can support both stem cell renewal and differentiation. However, responses to HA in culture are influenced by interaction with a range of cognate factors and receptors including components of blood serum supplements, which alter results. These may contribute to variation in cell batch production yield and phenotype as well as heighten the risks of adventitious pathogen transmission in the course of cell processing for therapeutic applications. MAIN: Here we characterise differentiation of a human embryo/pluripotent stem cell derived Mesenchymal Stromal Cell (hESC/PSC-MSC)-like cell population by culture on a planar surface coated with HA in serum-free media qualified for cell production for therapy. Resulting cells met minimum criteria of the International Society for Cellular Therapy for identification as MSC by expression of. CD90, CD73, CD105, and lack of expression for CD34, CD45, CD14 and HLA-II. They were positive for other MSC associated markers (i.e.CD166, CD56, CD44, HLA 1-A) whilst negative for others (e.g. CD271, CD71, CD146). In vitro co-culture assessment of MSC associated functionality confirmed support of growth of hematopoietic progenitors and inhibition of mitogen activated proliferation of lymphocytes from umbilical cord and adult peripheral blood mononuclear cells, respectively. Co-culture with immortalized THP-1 monocyte derived macrophages (Mɸ) concurrently stimulated with lipopolysaccharide as a pro-inflammatory stimulus, resulted in a dose dependent increase in pro-inflammatory IL6 but negligible effect on TNFα. To further investigate these functionalities, a bulk cell RNA sequence comparison with adult human bone marrow derived MSC and hESC substantiated a distinctive genetic signature more proximate to the former. CONCLUSION: Cultivation of human pluripotent stem cells on a planar substrate of HA in serum-free culture media systems is sufficient to yield a distinctive developmental mesenchymal stromal cell lineage with potential to modify the function of haematopoietic lineages in therapeutic applications.
Asunto(s)
Diferenciación Celular , Ácido Hialurónico , Células Madre Mesenquimatosas , Células Madre Pluripotentes , Humanos , Células Madre Mesenquimatosas/metabolismo , Células Madre Mesenquimatosas/citología , Ácido Hialurónico/farmacología , Ácido Hialurónico/metabolismo , Células Madre Pluripotentes/metabolismo , Células Madre Pluripotentes/citología , Medio de Cultivo Libre de Suero/farmacología , Linaje de la Célula , Células Cultivadas , Técnicas de Cultivo de Célula/métodos , Técnicas de CocultivoRESUMEN
The electrified aqueous/metal interface is critical in controlling the performance of energy conversion and storage devices, but an atomistic understanding of even basic interfacial electrochemical reactions challenges both experiment and computation. We report a combined simulation and experimental study of (reversible) ion-transfer reactions involved in anodic Ag corrosion/deposition, a model system for interfacial electrochemical processes generating or consuming ions. With the explicit modeling of the electrode potential and a hybrid implicit-explicit solvation model, the density functional theory calculations produce free energy curves predicting thermodynamics, kinetics, partial charge profiles, and reaction trajectories. The calculated (equilibrium) free energy barriers (0.2 eV), and their asymmetries, agree with experimental activation energies (0.4 eV) and transfer coefficients, which were extracted from temperature-dependent voltage-step experiments on Au-supported, Ag-nanocluster substrates. The use of Ag nanoclusters eliminates the convolution of the kinetics of Ag+(aq.) generation and transfer with those of nucleation or etch-pit formation. The results indicate that the barrier is controlled by the bias-dependent competition between partial solvation of the incipient ion, metal-metal bonding, and electrostatic stabilization by image charge, with the latter two factors weakened by stronger positive biases. We also report simulations of the bias-dependence of defect generation relevant to nucleating corrosion by removing an atom from a perfect Ag(100) surface, which is predicted to occur via a vacancy-adatom intermediate. Together, these experiments and calculations provide the first validated, accurate, molecular model of the central steps that govern the rates of important dissolution/deposition reactions broadly relevant across the energy sciences.
RESUMEN
The chemical bond is the cornerstone of chemistry, providing a conceptual framework to understand and predict the behavior of molecules in complex systems. However, the fundamental origin of chemical bonding remains controversial and has been responsible for fierce debate over the past century. Here, we present a unified theory of bonding, using a separation of electron delocalization effects from orbital relaxation to identify three mechanisms [node-induced confinement (typically associated with Pauli repulsion, though more general), orbital contraction, and polarization] that each modulate kinetic energy during bond formation. Through analysis of a series of archetypal bonds, we show that an exquisite balance of energy-lowering delocalizing and localizing effects are dictated simply by atomic electron configurations, nodal structure, and electronegativities. The utility of this unified bonding theory is demonstrated by its application to explain observed trends in bond strengths throughout the periodic table, including main group and transition metal elements.
RESUMEN
Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification of molecules in gas, solution, and solid-state phases. However, accurate prediction of chemical shifts using the gold-standard coupled cluster with singles, doubles, and perturbative triple excitations [CCSD(T)] method with a complete basis set (CBS) can be prohibitively expensive. By contrast, machine learning (ML) methods offer inexpensive alternatives for chemical shift predictions but are hampered by generalization to molecules outside the original training set. Here, we propose several new ideas in ML of the chemical shift prediction for H, C, N, and O that first introduce a novel feature representation, based on the atomic chemical shielding tensors within a molecular environment using an inexpensive quantum mechanics (QM) method, and train it to predict NMR chemical shieldings of a high-level composite theory that approaches the accuracy of CCSD(T)/CBS. In addition, we train the ML model through a new progressive active learning workflow that reduces the total number of expensive high-level composite calculations required while allowing the model to continuously improve on unseen data. Furthermore, the algorithm provides an error estimation, signaling potential unreliability in predictions if the error is large. Finally, we introduce a novel approach to keep the rotational invariance of the features using tensor environment vectors (TEVs) that yields a ML model with the highest accuracy compared to a similar model using data augmentation. We illustrate the predictive capacity of the resulting inexpensive shift machine learning (iShiftML) models across several benchmarks, including unseen molecules in the NS372 data set, gas-phase experimental chemical shifts for small organic molecules, and much larger and more complex natural products in which we can accurately differentiate between subtle diastereomers based on chemical shift assignments.
RESUMEN
A detailed chemical understanding of H2 interactions with binding sites in the nanoporous crystalline structure of metal-organic frameworks (MOFs) can lay a sound basis for the design of new sorbent materials. Computational quantum chemical calculations can aid in this quest. To set the stage, we review general thermodynamic considerations that control the usable storage capacity of a sorbent. We then discuss cluster modeling of H2 ligation at MOF binding sites using state-of-the-art density functional theory (DFT) calculations, and how the binding can be understood using energy decomposition analysis (EDA). Employing these tools, we illustrate the connections between the character of the MOF binding site and the associated adsorption thermodynamics using four experimentally characterized MOFs, highlighting the role of open metal sites (OMSs) in accessing binding strengths relevant to room temperature storage. The sorbents are MOF-5, with no open metal sites, Ni2(m-dobdc), containing Lewis acidic Ni(II) sites, Cu(I)-MFU-4l, containing π basic Cu(I) sites and V2Cl2.8(btdd), also containing π-basic V(II) sites. We next explore the potential for binding multiple H2 molecules at a single metal site, with thermodynamics useful for storage at ambient temperature; a materials design goal which has not yet been experimentally demonstrated. Computations on Ca2+ or Mg2+ bound to catecholate or Ca2+ bound to porphyrin show the potential for binding up to 4 H2; there is precedent for the inclusion of both catecholate and porphyrin motifs in MOFs. Turning to transition metals, we discuss the prediction that two H2 molecules can bind at V(II)-MFU-4l, a material that has been synthesized with solvent coordinated to the V(II) site. Additional calculations demonstrate binding three equivalents of hydrogen per OMS in Sc(I) or Ti(I)-exchanged MFU-4l. Overall, the results suggest promising prospects for experimentally realizing higher capacity hydrogen storage MOFs, if nontrivial synthetic and desolvation challenges can be overcome. Coupled with the unbounded chemical diversity of MOFs, there is ample scope for additional exploration and discovery.
RESUMEN
The catalytic transformation of C-H to C-N bonds offers rapid access to fine chemicals and high-performance materials, but achieving high selectivity from undirected aminations of unactivated C(sp3)-H bonds remains an outstanding challenge. We report the origins of the reactivity and selectivity of a Cu-catalyzed C-H amidation of simple alkanes. Using a combination of experimental and computational mechanistic studies and energy decomposition techniques, we uncover a switch in mechanism from inner-sphere to outer-sphere coupling between alkyl radicals and the active Cu(II) catalyst with increasing substitution of the alkyl radical. The combination of computational predictions and detailed experimental validation shows that simultaneous minimization of both Cu-C covalency and alkyl radical size increases the rate of reductive elimination and that both strongly electron-donating and electron-withdrawing substituents on the catalyst accelerate the selectivity-determining C-N bond formation process as a result of a change in mechanism. These findings offer design principles for the development of improved catalyst scaffolds for radical C-H functionalization reactions.
RESUMEN
Background: Augmented renal clearance (ARC) is a phenomenon observed in critically ill patients, leading to supraphysiologic drug clearance and concern for suboptimal antibiotic concentrations. The purpose of this study was to compare the clinical outcomes of our institutional protocolized antibiotic dosing regimen in critically ill patients with bacteremia and ARC compared with critically ill patients without ARC. Patients and Methods: We performed a retrospective study comparing the efficacy of an institutional protocolized antibiotic dosing regimen in critically ill patients with bacteremia and ARC compared with critically ill patients without ARC. The primary end point was in-hospital mortality. Secondary outcomes were intensive care unit (ICU) mortality, days requiring mechanical ventilation, ICU length of stay (LOS), hospital LOS, development of drug resistance to index antibiotic agent, and documented clearance of blood cultures within 72 hours. Results: There were 75 patients included in this study. Twenty percent of patients in the ARC group died in the hospital versus 31% in the non-ARC group (p = 0.26). The results for the ARC group versus the non-ARC group for the secondary outcomes of ICU mortality (20% vs. 26%; p = 0.56), ICU LOS (14.7 days vs. 7 days; p = 0.07), hospital LOS (28.3 days vs. 21.6 days; p = 0.03), days requiring mechanical ventilation (14 days vs. 12 days; p = 0.49), duration of antibiotic therapy (7.5 days vs. 9.0 days; p = 0.39), documented clearance of blood cultures within 72 hours (41% vs. 33%; p = 0.56), and the development of drug resistance to the index antibiotic agent (0% vs. 0%; p > 0.99) were also calculated. Conclusions: Among critically ill patients with bacteremia and ARC, there was no difference in in-hospital mortality compared with critically ill patients without ARC. There was a difference in hospital LOS, with a shorter duration of stay for the non-ARC group. There was no development of multi-drug-resistant organisms in either group.