Diffusion of curcumin in PLGA-based carriers for drug delivery: a molecular dynamics study.

CONTEXT: The rapid growth and diversification of drug delivery systems have been significantly supported by advancements in micro- and nano-technologies, alongside the adoption of biodegradable polymeric materials like poly(lactic-co-glycolic acid) (PLGA) as microcarriers. These developments aim to reduce toxicity and enhance target specificity in drug delivery. The use of in silico methods, particularly molecular dynamics (MD) simulations, has emerged as a pivotal tool for predicting the dynamics of species within these systems. This approach aids in investigating drug delivery mechanisms, thereby reducing the costs associated with design and prototyping. In this study, we focus on elucidating the diffusion mechanisms in curcumin-loaded PLGA particles, which are critical for optimizing drug release and efficacy in therapeutic applications. METHODS: We utilized MD to explore the diffusion behavior of curcumin in PLGA drug delivery systems. The simulations, executed with GROMACS, modeled curcumin molecules in a representative volume element of PLGA chains and water, referencing molecular structures from the Protein Data Bank and employing the CHARMM force field. We generated PLGA chains of varying lengths using the Polymer Modeler tool and arranged them in a bulk-like environment with Packmol. The simulation protocol included steps for energy minimization, T and p equilibration, and calculation of the isotropic diffusion coefficient from the mean square displacement. The Taguchi method was applied to assess the effects of hydration level, PLGA chain length, and density on diffusion. RESULTS: Our results provide insight into the influence of PLGA chain length, hydration level, and polymer density on the diffusion coefficient of curcumin, offering a mechanistic understanding for the design of efficient drug delivery systems. The sensitivity analysis obtained through the Taguchi method identified hydration level and PLGA density as the most significant input parameters affecting curcumin diffusion, while the effect of PLGA chain length was negligible within the simulated range. We provided a regression equation capable to accurately fit MD results. The regression equation suggests that increases in hydration level and PLGA density result in a decrease in the diffusion coefficient.

Enhancing ReaxFF for molecular dynamics simulations of lithium-ion batteries: an interactive reparameterization protocol.

Lithium-ion batteries (LIBs) have become an essential technology for the green economy transition, as they are widely used in portable electronics, electric vehicles, and renewable energy systems. The solid-electrolyte interphase (SEI) is a key component for the correct operation, performance, and safety of LIBs. The SEI arises from the initial thermal metastability of the anode-electrolyte interface, and the resulting electrolyte reduction products stabilize the interface by forming an electrochemical buffer window. This article aims to make a first-but important-step towards enhancing the parametrization of a widely-used reactive force field (ReaxFF) for accurate molecular dynamics (MD) simulations of SEI components in LIBs. To this end, we focus on Lithium Fluoride (LiF), an inorganic salt of great interest due to its beneficial properties in the passivation layer. The protocol relies heavily on various Python libraries designed to work with atomistic simulations allowing robust automation of all the reparameterization steps. The proposed set of configurations, and the resulting dataset, allow the new ReaxFF to recover the solid nature of the inorganic salt and improve the mass transport properties prediction from MD simulation. The optimized ReaxFF surpasses the previously available force field by accurately adjusting the diffusivity of lithium in the solid lattice, resulting in a two-order-of-magnitude improvement in its prediction at room temperature. However, our comprehensive investigation of the simulation shows the strong sensitivity of the ReaxFF to the training set, making its ability to interpolate the potential energy surface challenging. Consequently, the current formulation of ReaxFF can be effectively employed to model specific and well-defined phenomena by utilizing the proposed interactive reparameterization protocol to construct the dataset. Overall, this work represents a significant initial step towards refining ReaxFF for precise reactive MD simulations, shedding light on the challenges and limitations of ReaxFF force field parametrization. The demonstrated limitations emphasize the potential for developing more versatile and advanced force fields to upscale ab initio simulation through our interactive reparameterization protocol, enabling more accurate and comprehensive MD simulations in the future.

Method for predicting the wettability of micro-structured surfaces by continuum phase-field modelling.

Numerical prediction of material properties is attracting the attention of the scientific community and industry because of its usefulness in the design process. In the fields of fluid dynamics and microfluidics, several simulation methods have been proposed and adopted to evaluate the properties of surfaces and material interfaces, thanks to the increasing computational power available. However, despite the efforts made, a general and standardized methodology for implementing such methods is still lacking, thus requiring a trial-and-error approach for each new problem, making them difficult to implement and creating a bottleneck at the initial stage of surface design. Here, we report a validated protocol to evaluate the wettability of micro-structured surfaces with a phase-field model. Summarizing:•Simulating physical phenomena with multi-phase flows and moving contact lines can be a challenging task, due to the coupling among disparate length scales.•Using the Cahn-Hilliard diffuse-interface model, moving contact lines can be extensively investigated, although difficulties may arise when implementing numerical simulations, e.g., model parameter calibration, selection of boundary conditions, post-processing of fluid dynamics/equilibrium.•A method for employing this model and evaluating the physical consistency of the results is proposed here, considering the wettability of micro-structured surfaces as a case study.

Atomistic to Mesoscopic Modelling of Thermophysical Properties of Graphene-Reinforced Epoxy Nanocomposites.

This research addresses the need for a multiscale model for the determination of the thermophysical properties of nanofiller-enhanced thermoset polymer composites. Specifically, we analyzed the thermophysical properties of an epoxy resin containing bisphenol-A diglyceryl ether (DGEBA) as an epoxy monomer and dicyandiamide (DICY) and diethylene triamine (DETA) as cross-linking agents. The cross-linking process occurs at the atomistic scale through the formation of bonds among the reactive particles within the epoxy and hardener molecules. To derive the interatomic coarse-grained potential for the mesoscopic model and match the density of the material studied through atomic simulations, we employed the iterative Boltzmann inversion method. The newly developed coarse-grained molecular dynamics model effectively reproduces various thermophysical properties of the DGEBA-DICY-DETA resin system. Furthermore, we simulated nanocomposites made of the considered epoxy additivated with graphene nanofillers at the mesoscopic level and verified them against continuum approaches. Our results demonstrate that a moderate amount of nanofillers (up to 2 wt.%) increases the elastic modulus and thermal conductivity of the epoxy resin while decreasing the Poisson's ratio. For the first time, we present a coarse-grained model of DGEBA-DICY-DETA/graphene materials, which can facilitate the design and development of composites with tunable thermophysical properties for a potentially wide range of applications, e.g., automotive, aerospace, biomedical, or energy ones.

The impact of physicochemical features of carbon electrodes on the capacitive performance of supercapacitors: a machine learning approach.

Hybrid electric vehicles and portable electronic systems use supercapacitors for energy storage owing to their fast charging/discharging rates, long life cycle, and low maintenance. Specific capacitance is regarded as one of the most important performance-related characteristics of a supercapacitor's electrode. In the current study, Machine Learning (ML) algorithms were used to determine the impact of various physicochemical properties of carbon-based materials on the capacitive performance of electric double-layer capacitors. Published experimental datasets from 147 references (4899 data entries) were extracted and then used to train and test the ML models, to determine the relative importance of electrode material features on specific capacitance. These features include current density, pore volume, pore size, presence of defects, potential window, specific surface area, oxygen, and nitrogen content of the carbon-based electrode material. Additionally, categorical variables as the testing method, electrolyte, and carbon structure of the electrodes are considered as well. Among five applied regression models, an extreme gradient boosting model was found to best correlate those features with the capacitive performance, highlighting that the specific surface area, the presence of nitrogen doping, and the potential window are the most significant descriptors for the specific capacitance. These findings are summarized in a modular and open-source application for estimating the capacitance of supercapacitors given, as only inputs, the features of their carbon-based electrodes, the electrolyte and testing method. In perspective, this work introduces a new wide dataset of carbon electrodes for supercapacitors extracted from the experimental literature, also giving an instance of how electrochemical technology can benefit from ML models.

Mesoscopic Modeling and Experimental Validation of Thermal and Mechanical Properties of Polypropylene Nanocomposites Reinforced By Graphene-Based Fillers.

The development of nanocomposites relies on structure-property relations, which necessitate multiscale modeling approaches. This study presents a modeling framework that exploits mesoscopic models to predict the thermal and mechanical properties of nanocomposites starting from their molecular structure. In detail, mesoscopic models of polypropylene (PP)- and graphene-based nanofillers (graphene (Gr), graphene oxide (GO), and reduced graphene oxide (rGO)) are considered. The newly developed mesoscopic model for the PP/Gr nanocomposite provides mechanistic information on the thermal and mechanical properties at the filler-matrix interface, which can then be exploited to enhance the prediction accuracy of traditional continuum simulations by calibrating the thermal and mechanical properties of the filler-matrix interface. Once validated through a dedicated experimental campaign, this multiscale model demonstrates that with the modest addition of nanofillers (up to 2 wt %), the Young's modulus and thermal conductivity show up to 35 and 25% enhancement, respectively, whereas the Poisson's ratio slightly decreases. Among the different combinations tested, the PP/Gr nanocomposite shows the best mechanical properties, whereas PP/rGO demonstrates the best thermal conductivity. This validated mesoscopic model can contribute to the development of smart materials with enhanced mechanical and thermal properties based on polypropylene, especially for mechanical, energy storage, and sensing applications.

Ultrafast nano generation of acoustic waves in water via a single carbon nanotube.

Generation of ultra high frequency acoustic waves in water is key to nano resolution sensing, acoustic imaging and theranostics. In this context water immersed carbon nanotubes (CNTs) may act as an ideal optoacoustic source, due to their nanometric radial dimensions, peculiar thermal properties and broad band optical absorption. The generation mechanism of acoustic waves in water, upon excitation of both a single-wall (SW) and a multi-wall (MW) CNT with laser pulses of temporal width ranging from 5 ns down to ps, is theoretically investigated via a multiscale approach. We show that, depending on the combination of CNT size and laser pulse duration, the CNT can act as a thermophone or a mechanophone. As a thermophone, the CNT acts as a nanoheater for the surrounding water, which, upon thermal expansion, launches the pressure wave. As a mechanophone, the CNT acts as a nanopiston, its thermal expansion directly triggering the pressure wave in water. Activation of the mechanophone effect is sought to trigger few nanometers wavelength sound waves in water, matching the CNT acoustic frequencies. This is at variance with respect to the commonly addressed case of water-immersed single metallic nano-objects excited with ns laser pulses, where only the thermophone effect significantly contributes. The present findings might be of impact in fields ranging from nanoscale non-destructive testing to water dynamics at the meso to nanoscale.

Anisotropic Electrostatic Interactions in Coarse-Grained Water Models to Enhance the Accuracy and Speed-Up Factor of Mesoscopic Simulations.

Water models with realistic physical-chemical properties are essential to study a variety of biomedical processes or engineering technologies involving molecules or nanomaterials. Atomistic models of water are constrained by the feasible computational capacity, but calibrated coarse-grained (CG) ones can go beyond these limits. Here, we compare three popular atomistic water models with their corresponding CG model built using finite-size particles such as ellipsoids. Differently from previous approaches, short-range interactions are accounted for with the generalized Gay-Berne potential, while electrostatic and long-range interactions are computed from virtual charges inside the ellipsoids. Such an approach leads to a quantitative agreement between the original atomistic models and their CG counterparts. Results show that a timestep of up to 10 fs can be achieved to integrate the equations of motion without significant degradation of the physical observables extracted from the computed trajectories, thus unlocking a significant acceleration of water-based mesoscopic simulations at a given accuracy.

Machine learning and materials modelling interpretation of in vivo toxicological response to TiO2 nanoparticles library (UV and non-UV exposure).

Assessing the risks of nanomaterials/nanoparticles (NMs/NPs) under various environmental conditions requires a more systematic approach, including the comparison of effects across many NMs with identified different but related characters/descriptors. Hence, there is an urgent need to provide coherent (eco)toxicological datasets containing comprehensive toxicity information relating to a diverse spectra of NPs characters. These datasets are test benches for developing holistic methodologies with broader applicability. In the present study we assessed the effects of a custom design Fe-doped TiO2 NPs library, using the soil invertebrate Enchytraeus crypticus (Oligochaeta), via a 5-day pulse via aqueous exposure followed by a 21-days recovery period in soil (survival, reproduction assessment). Obviously, when testing TiO2, realistic conditions should include UV exposure. The 11 Fe-TiO2 library contains NPs of size range between 5-27 nm with varying %Fe (enabling the photoactivation of TiO2 at energy wavelengths in the visible-light range). The NPs were each described by 122 descriptors, being a mixture of measured and atomistic model descriptors. The data were explored using single and univariate statistical methods, combined with machine learning and multiscale modelling techniques. An iterative pruning process was adopted for identifying automatically the most significant descriptors. TiO2 NPs toxicity decreased when combined with UV. Notably, the short-term water exposure induced lasting biological responses even after longer-term recovery in clean exposure. The correspondence with Fe-content correlated with the band-gap hence the reduction of UV oxidative stress. The inclusion of both measured and modelled materials data benefitted the explanation of the results, when combined with machine learning.

Effect of water nanoconfinement on the dynamic properties of paramagnetic colloidal complexes.

The anomalous behavior of confined water at the nanoscale has remarkable implications in a number of nanotechnological applications. In this work, we analyze the effect of water self-diffusion on the dynamic properties of a solvated gadolinium-based paramagnetic complex, typically used for contrast enhancement in magnetic resonance imaging. In particular, we examine the effect of silica-based nanostructures on water behavior in the proximity of the paramagnetic complex via atomistic simulations, and interpret the resulting tumbling dynamics in the light of the local solvent modification based on the Lipari-Szabo formalism and of the fractional Stokes-Einstein relation. It is found that the local water confinement induces an increased "stiffness" on the outer sphere of the paramagnetic complex, which eventually reduces its tumbling properties. These model predictions are found to explain well the relaxivity enhancement observed experimentally by confining paramagnetic complexes into porous nanoconstructs, and thus offer mechanistic guidelines to design improved contrast agents for imaging applications.

Multistage and passive cooling process driven by salinity difference.

Space cooling in buildings is anticipated to rise because of an increasing thermal comfort demand worldwide, and this calls for cost-effective and sustainable cooling technologies. We present a proof-of-concept multistage device, where a net cooling capacity and a temperature difference are demonstrated as long as two water solutions at disparate salinity are maintained. Each stage is made of two hydrophilic layers separated by a hydrophobic membrane. An imbalance in water activity in the two layers naturally causes a non-isothermal vapor flux across the membrane without requiring any mechanical ancillaries. One prototype of the device developed a specific cooling capacity of up to 170 W m-2 at a vanishing temperature difference, considering a 3.1 mol/kg calcium chloride solution. To provide perspective, if successfully up-scaled, this concept may help satisfy at least partially the cooling needs in hot, humid regions with naturally available salinity gradients.

Mechanistic correlation between water infiltration and framework hydrophilicity in MFI zeolites.

Hydrophobic zeolites are nanoporous materials that are attracting an increasing interest, especially for catalysis, desalination, energy storage and biomedical applications. Nevertheless, a more profound understanding and control of water infiltration in their nanopores is still desirable to rationally design zeolite-based materials with tailored properties. In this work, both atomistic simulations and previous experimental data are employed to investigate water infiltration in hydrophobic MFI zeolites with different concentration of hydrophilic defects. Results show that limited concentrations of defects (e.g. 1%) induce a change in the shape of infiltration isotherms (from type-V to type-I), which denotes a sharp passage from typical hydrophobic to hydrophilic behavior. A correlation parametrized on both energy and geometric characteristics of the zeolite (infiltration model) is then adopted to interpolate the infiltration isotherms data by means of a limited number of physically-meaningful parameters. Finally, the infiltration model is combined with the water-zeolite interaction energy computed by simulations to correlate the water intrusion mechanism with the atomistic details of the zeolite crystal, such as defects concentration, distribution and hydrophilicity. The suggested methodology may allow a faster (more than one order of magnitude) and more systematic preliminary computational screening of innovative zeolite-based materials for energy storage, desalination and biomedical purposes.

Nano-metering of Solvated Biomolecules Or Nanoparticles from Water Self-Diffusivity in Bio-inspired Nanopores.

Taking inspiration from the structure of diatom algae frustules and motivated by the need for new detecting strategies for emerging nanopollutants in water, we analyze the potential of nanoporous silica tablets as metering devices for the concentration of biomolecules or nanoparticles in water. The concept relies on the different diffusion behavior that water molecules exhibit in bulk and nanoconfined conditions, e.g., in nanopores. In this latter situation, the self-diffusion coefficient of water reduces according to the geometry and surface properties of the pore and to the concentration of suspended biomolecules or nanoparticles in the pore, as extensively demonstrated in a previous study. Thus, for a given pore-liquid system, the self-diffusivity of water in nanopores filled with biomolecules or nanoparticles provides an indirect measure of their concentration. Using molecular dynamics and previous results from the literature, we demonstrate the correlation between the self-diffusion coefficient of water in silica nanopores and the concentration of proteins or nanoparticles contained therein. Finally, we estimate the time required for the nanoparticles to fill the nanopores, in order to assess the practical feasibility of the overall nano-metering protocol. Results show that the proposed approach may represent an alternative method for assessing the concentration of some classes of nanopollutants or biomolecules in water.

From GROMACS to LAMMPS: GRO2LAM : A converter for molecular dynamics software.

Atomistic simulations have progressively attracted attention in the study of physical-chemical properties of innovative nanomaterials. GROMACS and LAMMPS are currently the most widespread open-source software for molecular dynamics simulations thanks to their good flexibility, numerous functionalities and responsive community support. Nevertheless, the very different formats adopted for input and output files are limiting the possibility to transfer GROMACS simulations to LAMMPS. In this article, we present GRO2LAM, a modular and open-source Python 2.7 code for rapidly translating input files and parameters from GROMACS to LAMMPS format. The robustness of the tool has been assessed by comparing the simulation results obtained by GROMACS and LAMMPS, after the format conversion by GRO2LAM. Specifically, three nanoscale configurations of interest in both engineering and biomedical fields are studied, namely a carbon nanotube, an iron oxide nanoparticle, and a protein immersed in water. In perspective, GRO2LAM may be the first step to achieve a full interoperability between molecular dynamics software. This would allow to easily exploit their complementary potentialities and post-processing functionalities. Moreover, GRO2LAM could facilitate the cross-check of simulation results, guaranteeing the reproducibility of molecular dynamics models and testing their robustness. Graphical Abstract GRO2LAM, a modular and open-source Python code for rapidly translating input files and parameters from GROMACS to LAMMPS format.

Coffee-based colloids for direct solar absorption.

Despite their promising thermo-physical properties for direct solar absorption, carbon-based nanocolloids present some drawbacks, among which the unpleasant property of being potentially cytotoxic and harmful to the environment. In this work, a sustainable, stable and inexpensive colloid based on coffee is synthesized and its photo-thermal properties investigated. The proposed colloid consists of distilled water, Arabica coffee, glycerol and copper sulphate, which provide enhanced properties along with biocompatibility. The photo-thermal performance of the proposed fluid for direct solar absorption is analysed for different dilutions and compared with that of a traditional flat-plate collector. Tailor-made collectors, opportunely designed and realized via 3D-printing technique, were used for the experimental tests. The results obtained in field conditions, in good agreement with two different proposed models, show similar performance of the volumetric absorption using the proposed coffee-based colloids as compared to the classical systems based on a highly-absorbing surface. These results may encourage further investigations on simple, biocompatible and inexpensive colloids for direct solar absorption.

Mesoscopic Moment Equations for Heat Conduction: Characteristic Features and Slow-Fast Mode Decomposition.

In this work, we derive different systems of mesoscopic moment equations for the heat-conduction problem and analyze the basic features that they must hold. We discuss two- and three-equation systems, showing that the resulting mesoscopic equation from two-equation systems is of the telegraphist's type and complies with the Cattaneo equation in the Extended Irreversible Thermodynamics Framework. The solution of the proposed systems is analyzed, and it is shown that it accounts for two modes: a slow diffusive mode, and a fast advective mode. This latter additional mode makes them suitable for heat transfer phenomena on fast time-scales, such as high-frequency pulses and heat transfer in small-scale devices. We finally show that, if proper initial conditions are provided, the advective mode disappears, and the solution of the system tends asymptotically to the transient solution of the classical parabolic heat-conduction equation.

Efficient steam generation by inexpensive narrow gap evaporation device for solar applications.

Technologies for solar steam generation with high performance can help solving critical societal issues such as water desalination or sterilization, especially in developing countries. Very recently, we have witnessed a rapidly growing interest in the scientific community proposing sunlight absorbers for direct conversion of liquid water into steam. While those solutions can possibly be of interest from the perspective of the involved novel materials, in this study we intend to demonstrate that efficient steam generation by solar source is mainly due to a combination of efficient solar absorption, capillary water feeding and narrow gap evaporation process, which can also be achieved through common materials. To this end, we report both numerical and experimental evidence that advanced nano-structured materials are not strictly necessary for performing sunlight driven water-to-vapor conversion at high efficiency (i.e. ≥85%) and relatively low optical concentration (≈10 suns). Coherently with the principles of frugal innovation, those results unveil that solar steam generation for desalination or sterilization purposes may be efficiently obtained by a clever selection and assembly of widespread and inexpensive materials.

Thermal transport across nanoparticle-fluid interfaces: the interplay of interfacial curvature and nanoparticle-fluid interactions.

We investigate the general dependence of the thermal transport across nanoparticle-fluid interfaces using molecular dynamics computations. We show that the thermal conductance depends strongly both on the wetting characteristics of the nanoparticle-fluid interface and on the nanoparticle size. Strong nanoparticle-fluid interactions, leading to full wetting states in the host fluid, result in high thermal conductances and efficient interfacial transport of heat. Weak interactions result in partial drying or full drying states, and low thermal conductances. The variation of the thermal conductance with particle size is found to depend on the fluid-nanoparticle interactions. Strong interactions coupled with large interfacial curvatures lead to optimum interfacial heat transport. This complex dependence can be modelled using an equation that includes the interfacial curvature as a parameter. In this way, we rationalise the existing experimental and computer simulation results and show that the thermal transport across nanoscale interfaces is determined by the correlations of both interfacial curvature and nanoparticle-fluid interactions.

Thermal transport phenomena in nanoparticle suspensions.

Nanoparticle suspensions in liquids have received great attention, as they may offer an approach to enhance thermophysical properties of base fluids. A good variety of applications in engineering and biomedicine has been investigated with the aim of exploiting the above potential. However, the multiscale nature of nanosuspensions raises several issues in defining a comprehensive modelling framework, incorporating relevant molecular details and much larger scale phenomena, such as particle aggregation and their dynamics. The objectives of the present topical review is to report and discuss the main heat and mass transport phenomena ruling macroscopic behaviour of nanosuspensions, arising from molecular details. Relevant experimental results are included and properly put in the context of recent observations and theoretical studies, which solved long-standing debates about thermophysical properties enhancement. Major transport phenomena are discussed and in-depth analysis is carried out for highlighting the role of geometrical (nanoparticle shape, size, aggregation, concentration), chemical (pH, surfactants, functionalization) and physical parameters (temperature, density). We finally overview several computational techniques available at different scales with the aim of drawing the attention on the need for truly multiscale predictive models. This may help the development of next-generation nanoparticle suspensions and their rational use in thermal applications.

Interplay between hydrophilicity and surface barriers on water transport in zeolite membranes.

A comprehensive understanding of molecular transport within nanoporous materials remains elusive in a broad variety of engineering and biomedical applications. Here, experiments and atomistic simulations are synergically used to elucidate the non-trivial interplay between nanopore hydrophilicity and surface barriers on the overall water transport through zeolite crystals. At these nanometre-length scales, these results highlight the dominating effect of surface imperfections with reduced permeability on the overall water transport. A simple diffusion resistance model is shown to be sufficient to capture the effects of both intracrystalline and surface diffusion resistances, thus properly linking simulation to experimental evidence. This work suggests that future experimental work should focus on eliminating/overcoming these surface imperfections, which promise an order of magnitude improvement in permeability.