Increased Absorption of Thyroxine in a Murine Model of Hypothyroidism Using Water/CO2 Nanobubbles.

Thyroxine (T4) is a drug extensively utilized for the treatment of hypothyroidism. However, the oral absorption of T4 presents certain limitations. This research investigates the efficacy of CO2 nanobubbles in water as a potential oral carrier for T4 administration to C57BL/6 hypothyroid mice. Following 18 h of fasting, the formulation was administered to the mice, demonstrating that the combination of CO2 nanobubbles and T4 enhanced the drug's absorption in blood serum by approximately 40%. To comprehend this observation at a molecular level, we explored the interaction mechanism through which T4 engages with the CO2 nanobubbles, employing molecular simulations, semi-empirical quantum mechanics, and PMF calculations. Our simulations revealed a high affinity of T4 for the water-gas interface, driven by additive interactions between the hydrophobic region of T4 and the gas phase and electrostatic interactions of the polar groups of T4 with water at the water-gas interface. Concurrently, we observed that at the water-gas interface, the cluster of T4 formed in the water region disassembles, contributing to the drug's bioavailability. Furthermore, we examined how the gas within the nanobubbles aids in facilitating the drug's translocation through cell membranes. This research contributes to a deeper understanding of the role of CO2 nanobubbles in drug absorption and subsequent release into the bloodstream. The findings suggest that utilizing CO2 nanobubbles could enhance T4 bioavailability and cell permeability, leading to more efficient transport into cells. Additional research opens the possibility of employing lower concentrations of this class of drugs, thereby potentially reducing the associated side effects due to poor absorption.

Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: Application to liquid water.

Knowledge of the electronic stopping curve for swift ions, Se(v), particularly around the Bragg peak, is important for understanding radiation damage. Experimentally, however, the determination of such a feature for light ions is very challenging, especially in disordered systems such as liquid water and biological tissue. Recent developments in real-time time-dependent density functional theory (rt-TDDFT) have enabled the calculation of Se(v) along nm-sized trajectories. However, it is still a challenge to obtain a meaningful statistically averaged Se(v) that can be compared to observations. In this work, taking advantage of the correlation between the local electronic structure probed by the projectile and the distance from the projectile to the atoms in the target, we devise a trajectory pre-sampling scheme to select, geometrically, a small set of short trajectories to accelerate the convergence of the averaged Se(v) computed via rt-TDDFT. For protons in liquid water, we first calculate the reference probability distribution function (PDF) for the distance from the proton to the closest oxygen atom, ϕR(rp→O), for a trajectory of a length similar to those sampled experimentally. Then, short trajectories are sequentially selected so that the accumulated PDF reproduces ϕR(rp→O) to increasingly high accuracy. Using these pre-sampled trajectories, we demonstrate that the averaged Se(vp) converges in the whole velocity range with less than eight trajectories, while other averaging methods using randomly and uniformly distributed trajectories require approximately ten times the computational effort. This allows us to compare the Se(vp) curve to experimental data and assess widely used empirical tables based on Bragg's rule.

First principles simulation of damage to solvated nucleotides due to shock waves.

We present a first-principles molecular dynamics study of the effect of shock waves (SWs) propagating in a model biological medium. We find that the SW can cause chemical modifications through varied and complex mechanisms, in particular, phosphate-sugar and sugar-base bond breaks. In addition, the SW promotes the dissociation of water molecules, thus enhancing the ionic strength of the medium. Freed protons can hydrolyze base and sugar rings previously opened by the shock. However, many of these events are only temporary, and bonds reform rapidly. Irreversible damage is observed for pressures above 15-20 GPa. These results are important to gain a better understanding of the microscopic damage mechanisms underlying cosmic-ray irradiation in space and ion-beam cancer therapy.

Mechanisms of Iodide⁻Triiodide Exchange Reactions in Ionic Liquids: A Reactive Molecular-Dynamics Exploration.

Efficient charge transport has been observed in iodine-doped, iodide-based room-temperature ionic liquids, yielding high ionic conductivity. To elucidate preferred mechanistic pathways for the iodide ( I - )-to-triiodide ( I 3 - ) exchange reactions, we have performed 10 ns reactive molecular-dynamics calculations in the liquid state for 1-butyl-3-methylimidazolium iodide ([BMIM][I]) at 450 to 750 K. Energy-barrier distributions for the iodine-swapping process were determined as a function of temperature, employing a charge-reassignment scheme drawn in part from electronic-structure calculations. Bond-exchange events were observed with rate-determining energy barriers ranging from ~0.19 to 0.23 ± 0.06 eV at 750 and 450 K, respectively, with an approximately Arrhenius temperature dependence for iodine self-diffusivity and reaction kinetics, although diffusion dominates/limits the bond-exchange events. This charge transfer is not dissimilar in energetics to those in solid-state superionic conductors.

Cement As a Waste Form for Nuclear Fission Products: The Case of (90)Sr and Its Daughters.

One of the main challenges faced by the nuclear industry is the long-term confinement of nuclear waste. Because it is inexpensive and easy to manufacture, cement is the material of choice to store large volumes of radioactive materials, in particular the low-level medium-lived fission products. It is therefore of utmost importance to assess the chemical and structural stability of cement containing radioactive species. Here, we use ab initio calculations based on density functional theory (DFT) to study the effects of (90)Sr insertion and decay in C-S-H (calcium-silicate-hydrate) in order to test the ability of cement to trap and hold this radioactive fission product and to investigate the consequences of its ß-decay on the cement paste structure. We show that (90)Sr is stable when it substitutes the Ca(2+) ions in C-S-H, and so is its daughter nucleus (90)Y after ß-decay. Interestingly, (90)Zr, daughter of (90)Y and final product in the decay sequence, is found to be unstable compared to the bulk phase of the element at zero K but stable when compared to the solvated ion in water. Therefore, cement appears as a suitable waste form for (90)Sr storage.

Protection of DNA against low-energy electrons by amino acids: a first-principles molecular dynamics study.

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

Dynamics of excess electronic charge in aliphatic ionic liquids containing the bis(trifluoromethylsulfonyl)amide anion.

In a recent article (J. Am. Chem. Soc. 2011, 133, 20186) we investigated the initial spatial distribution of dry excess electrons in a series of room-temperature ionic liquids (RTILs). Perhaps unexpectedly, we found that in some alkylammonium-based systems the excess negative charge resided on anions and not on the positive cations. Following on these results, in the current paper we describe the time evolution of an excess electronic charge introduced in alkylammonium- and pyrrolidinium-based ionic liquids coupled with the bis(trifluoromethylsulfonyl)amide ([Tf2N(-)]) anion. We find that on a 50 fs time scale an initially delocalized excess electron localizes on a single [Tf2N(-)] anion which begins a fragmentation process. Low-energy transitions have a very different physical origin on the several femtoseconds time scale when compared to what occurs on the picosecond time scale. At time zero, these are intraband transitions of the excess electron. However after 40 fs when the excess electronic charge localizes on a single anion, these transitions disappear, and the spectrum is dominated by electron-transfer transitions between the fragments of the doubly charged breaking anion.

Excess electron interactions with solvated DNA nucleotides: strand breaks possible at room temperature.

When biological matter is subjected to ionizing radiation, a wealth of secondary low-energy (<20 eV) electrons are produced. These electrons propagate inelastically, losing energy to the medium until they reach energies low enough to localize in regions of high electron affinity. We have recently shown that in fully solvated DNA fragments, nucleobases are particularly attractive for such excess electrons. The next question is what is their longer-term effect on DNA. It has been advocated that they can lead to strand breaks by cleavage of the phosphodiester C(3')-O(3') bond. Here we present a first-principles study of free energy barriers for the cleavage of this bond in fully solvated nucleotides. We have found that except for dAMP, the barriers are on the order of 6 kcal/mol, suggesting that bond cleavage is a regular feature at 300 K. Such low barriers are possible only as a result of solvent and thermal fluctuations. These findings support the notion that low-energy electrons can indeed lead to strand breaks in DNA.

Nonadiabatic forces in ion-solid interactions: the initial stages of radiation damage.

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

Inelastic scattering of electrons in water from first principles: cross sections and inelastic mean free path for use in Monte Carlo track-structure simulations of biological damage.

Modelling the inelastic scattering of electrons in water is fundamental, given their crucial role in biological damage. In Monte Carlo track-structure (MC-TS) codes used to assess biological damage, the energy loss function (ELF), from which cross sections are extracted, is derived from different semi-empirical optical models. Only recently have first ab initio results for the ELF and cross sections in water become available. For benchmarking purpose, in this work, we present ab initio linear-response time-dependent density functional theory calculations of the ELF of liquid water. We calculated the inelastic scattering cross sections, inelastic mean free paths, and electronic stopping power and compared our results with recent calculations and experimental data showing a good agreement. In addition, we provide an in-depth analysis of the contributions of different molecular orbitals, species and orbital angular momenta to the total ELF. Moreover, we present single-differential cross sections computed for each molecular orbital channel, which should prove useful for MC-TS simulations.

Characterization and Performance Enhancement of Cement-Based Thermoelectric Materials.

Thermoelectric materials enable the direct conversion of thermal to electrical energy. One application of this is ambient heat energy harvesting where relatively stable temperature gradients existing between the inside and outside of a building could be utilized to produce electricity. Buildings can thus change from energy consumers to energy generators. This could ultimately help reduce the surface temperatures and energy consumption of buildings, especially in urban areas. In this paper, research work carried out on developing and characterizing a cement-based thermoelectric material is presented. Cement-based samples are doped with different metal oxides (Bi2O3 and Fe2O3) to enhance their thermoelectric properties, which are defined through their Seebeck coefficient, electrical conductivity and thermal conductivity. The study also discusses the positive impact of moisture content on the electrical conductivity.

Dry excess electrons in room-temperature ionic liquids.

In this article, we describe general trends to be expected at short times when an excess electron is generated or injected in different room-temperature ionic liquids (RTILs). Perhaps surprisingly, the excess electron does not localize systematically on the positively charged cations. Rather, the excess charge localization pattern is determined by the cation and anion HOMO/LUMO gaps and, more importantly, by their relative LUMO alignments. As revealed by experiments, the short-time (ps/ns) transient UV spectrum of excess electrons in RTILs is often characterized by two bands, a broad band at low energies (above 1000 nm) and another weaker band at higher energies (around 400 nm). Our calculations show that the dry or presolvated electron spectrum (fs) also has two similar features. The broad band at low energies is due to transitions between electronic states with similar character on ions of the same class but in different locations of the liquid. The lower-intensity band at higher energies is due to transitions in which the electron is promoted to electronic states of different character, in some cases on counterions. Depending on the chemical nature of the RTIL, and especially on the anions, excess electrons can localize on cations or anions. Our findings hint at possible design strategies for controlling electron localization, where electron transfer or transport across species can be facilitated or blocked depending on the alignment of the electronic levels of the individual species.

Excess electron localization in solvated DNA bases.

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Dispersion interactions in room-temperature ionic liquids: results from a non-empirical density functional.

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Román-Pérez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by ≈7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Pópolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields.

Effect of intrinsic defects on the thermal conductivity of PbTe from classical molecular dynamics simulations.

Despite being the archetypal thermoelectric material, still today some of the most exciting advances in the efficiency of these materials are being achieved by tuning the properties of PbTe. Its inherently low lattice thermal conductivity can be lowered to its fundamental limit by designing a structure capable of scattering phonons over a wide range of length scales. Intrinsic defects, such as vacancies or grain boundaries, can and do play the role of these scattering sites. Here we assess the effect of these defects by means of molecular dynamics simulations. For this we purposely parametrize a Buckingham potential that provides an excellent description of the thermal conductivity of this material over a wide temperature range. Our results show that intrinsic point defects and grain boundaries can reduce the lattice conductivity of PbTe down to a quarter of its bulk value. By studying the size dependence we also show that typical defect concentrations and grain sizes realized in experiments normally correspond to the bulk lattice conductivity of pristine PbTe.

Hollow Gold Nanoparticles Produced by Femtosecond Laser Irradiation.

Metallic hollow nanoparticles exhibit interesting optical properties that can be controlled by geometrical parameters. Irradiation with femtosecond laser pulses has emerged recently as a valuable tool for reshaping and size modification of plasmonic metal nanoparticles, thereby enabling the synthesis of nanostructures with unique morphologies. In this Letter, we use classical molecular dynamics simulations to investigate the solid-to-hollow conversion of gold nanoparticles upon femtosecond laser irradiation. Here, we suggest an efficient method for producing hollow nanoparticles under certain specific conditions, namely that the particles should be heated to a maximum temperature between 2500 and 3500 K, followed by a fast quenching to room temperature, with cooling rates lower than 120 ps. Therefore, we describe the experimental conditions for efficiently producing hollow nanoparticles, opening a broad range of possibilities for applications in key areas, such as energy storage and catalysis.

On the role of magnesium in a LiF:Mg,Ti thermoluminescent dosimeter.

LiF doped with Mg and Ti is the most widely used thermoluminescent (TL) dosimeter for a large variety of applications. It has been argued that the Mg dopant is the most important defect in the TL process. Besides the common F-centre defects in LiF, optical absorption measurements have suggested the presence of Mg-related absorption bands at 380 nm (3.26 eV), and 310 nm (4.0 eV) when LiF:Mg is exposed to ionizing radiation, whose origin is not yet well understood. This work presents an investigation of the role of defects induced by Mg interstitials in LiF through electronic structure calculations. The calculations show that Mg interstitials induce a local lattice distortion characterized by the displacement of two opposite fluorine atoms, adjacent to the magnesium, away from their original sites by an average distance of 0.6 Å each, while the closest Li atoms are displaced by 0.1 Å. This defect introduces electronic states in the band-gap that can trap excess electrons produced during irradiation, thus enhancing the efficiency of the detector. Holes, on the other hand, are created and trapped in orbitals of mainly Mg-3s character. Additionally, the results suggest that irradiation can simultaneously remove a Li atom nearby a Mg interstitial; substitute a Li by a Mg atom or create a Li vacancy plus a Mg substitutional, giving rise to defects within the LiF gap that are more stable thermodynamically than the Mg interstitial itself. Interestingly, under irradiation the energy levels obtained for LiF:Mg-Lisub + e - (3.486 eV) and LiF:Mg + e - (4.224 eV) defects are very close to the experimental absorption bands.

Solvation Effects on Dissociative Electron Attachment to Thymine.

Ionizing radiation can excite the cellular medium to produce secondary electrons that can subsequently cause damage to DNA. The damage is believed to occur via dissociative electron attachment (DEA). In DEA, the electron is captured by a molecule in a resonant antibonding state and a transient negative ion is formed. If this ion survives against electron autodetachment, then bonds within the molecule may dissociate as energy is transferred from the electronic degrees of freedom into vibrational modes of the molecule. We present a model for studying the effect that transferring kinetic energy into the vibrational modes of a molecule in this way has on a DNA nucleobase. We show that when the base is in an aqueous environment, dissociation is affected by interactions with the surrounding water molecules. In particular, hydrogen bonding between the nucleobase and the solvent can suppress the dissociative channel.

Insights into mechanochemical reactions at the molecular level: simulated indentations of aspirin and meloxicam crystals.

Although solvent-free mechanochemical synthesis continues to gain ever greater importance, the molecular scale processes that occur during such reactions remain largely uncharacterised. Here, we apply computational modelling to indentations between particles of crystals of aspirin and meloxicam under a variety of conditions to mimic the early stages of their mechanochemical cocrystallisation reaction. The study also extends to the effects of the presence of small amounts of solvent. It is found that, despite the solid crystalline nature of the reactants and the presence of little or no solvent, mixing occurs readily at the molecular level even during relatively low-energy collisions. When indented crystals are subsequently drawn apart, a connective neck formed by a mixture of the reactant molecules is observed, suggesting plastic-like behaviour of the reacting materials. Overall the work reveals some striking new insights including (i) relatively facile mixing of crystals under solvent-free conditions, (ii) no appreciable local temperature increases, (iii) localised amorphisation at the contact region and neck of the reacting crystals, and (iv) small amounts of solvent have relatively little effect during this early stage of the reaction, suggesting that their accelerating effect on the reaction may be exerted at later stages.

Polarization relaxation in an ionic liquid confined between electrified walls.

The response of a room temperature molten salt to an external electric field when it is confined to a nanoslit is studied by molecular dynamics simulations. The fluid is confined between two parallel and oppositely charged walls, emulating two electrified solid-liquid interfaces. Attention is focused on structural, electrostatic, and dynamical properties, which are compared with those of the nonpolarized fluid. It is found that the relaxation of the electrostatic potential, after switching the electric field off, occurs in two stages. A first, subpicosecond process accounts for 80% of the decay and is followed by a second subdiffusive process with a time constant of 8 ps. Diffusion is not involved in the relaxation, which is mostly driven by small anion translations. The relaxation of the polarization in the confined system is discussed in terms of the spectrum of charge density fluctuations in the bulk.