Search

Single-component organic solar cells (SCOSCs) have attracted extensive attention due to their simplified device manufacturing and excellent stability. However, the relationship between morphology and charge carrier mobility in the active layers of SCOSCs is not well understood. In this work, we present a comprehensive investigation on this issue by studying four dyads (fullerenes as acceptor units) used as materials of active layers in small-molecule single-component organic solar cells (SM-SCOSCs), in which dyad 4 created the record of power conversion efficiency (PCE) of SM-SCOSC until now. Utilizing a multiscale theoretical approach, the results identify that the acceptor-acceptor stacking is dominant in amorphous films, significantly improving electron mobility and lowering hole mobility. We also find the importance of achieving a balance between electron and hole mobility to further improve PCE of SM-SCOSC because dyad 4 exhibits a more balanced electron/hole mobility than the other three molecules. These findings indicate the importance of tuning and enhancing donor-donor and acceptor-acceptor stacking simultaneously, offering insights for the design and optimization of future SM-SCOSC., (© 2024 Wiley-VCH GmbH.)

Using solid nanoparticles (NPs) as catalysts is the most effective method to achieve catalytic growth of single-walled carbon nanotubes (SWCNTs) with ultrapure chirality. Until now, SWCNTs with a suitable chirality purity have not been prepared in experiments. That is, the evolution of solid NPs during the catalytic growth of SWCNTs is in contradiction with the original concept of a changeless structure. Hence, in this work, the evolution mechanism of solid cobalt NPs during the nucleation process of SWCNTs is analyzed through molecular dynamics. Similar to the experimental observations, the results show that a drastic structural fluctuation of the NPs occurs during the nucleation of SWCNTs. This structural fluctuation is caused by the fact that the elastic strain energy and surface energy of the NPs can be tuned when a carbon gradient exists between the subsurface and interior of the NP. Furthermore, such a carbon gradient can be reduced by changing the carbon feeding rate. This work not only reveals the evolution mechanism of solid catalysts during the nucleation of SWCNTs but also provides prospects for realizing solid catalysts with a changeless structure by tuning the experimental parameters., (© 2024 Wiley‐VCH GmbH.)

Context: Nanoscrolls are tube-shaped structures formed when a sheet or ribbon of material is rolled into a cylinder, creating a hollow tube with a diameter on the nanoscale, similar to the papyrus. Carbon nanoscrolls have unique properties that make them useful in various applications, such as energy storage, catalysis, and drug delivery. In this study, we employed classical molecular dynamics simulations to investigate the formation and stability of nanoscrolls composed of graphene and hexagonal boron nitride (hBN) nanoribbons. Using a carbon nanotube (CNT) as a template to trigger their collapsing, we found that graphene/graphene, graphene/hBN, and hBN/hBN could form CNT-wrapped nanoscrolls at ultrafast speeds. We also confirmed that these nanoscrolls are thermally stable and discussed the other products formed from the interaction of these complexes and their temperature dependence. Gr/Gr and hBN/Gr nanoscrolls exhibit similar interlayer distances, while hBN/hBN nanoscrolls have wider interlayer distances than the other two composite nanoscrolls. These features suggest that hBN/hBN composite nanoscrolls could more efficiently capture small molecules because of their greater interlayer spacing., Methods: We conducted molecular dynamics simulations using the Forcite package in the Biovia Materials Studio software, which employs the Universal and Dreiding force fields. We considered an NVT ensemble with a fixed time step of 1.0 fs for a duration of 500 ps. The velocity Verlet algorithm was adopted to integrate the equations of motion of the entire system. We employed the Nosé-Hoover-Langevin thermostat to control the system temperature. The simulations were carried out without periodic boundary conditions, so there was no pressure coupling., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Among various porous solids for gas separation and purification, metal-organic frameworks (MOFs) are promising materials that potentially combine high CO 2 uptake and CO 2 /N 2 selectivity. So far, within the hundreds of thousands of MOF structures known today, it remains a challenge to computationally identify the best suited species. First principle-based simulations of CO 2 adsorption in MOFs would provide the necessary accuracy; however, they are impractical due to the high computational cost. Classical force field-based simulations would be computationally feasible; however, they do not provide sufficient accuracy. Thus, the entropy contribution that requires both accurate force fields and sufficiently long computing time for sampling is difficult to obtain in simulations. Here, we report quantum-informed machine-learning force fields (QMLFFs) for atomistic simulations of CO 2 in MOFs. We demonstrate that the method has a much higher computational efficiency (∼1000×) than the first-principle one while maintaining the quantum-level accuracy. As a proof of concept, we show that the QMLFF-based molecular dynamics simulations of CO 2 in Mg-MOF-74 can predict the binding free energy landscape and the diffusion coefficient close to experimental values. The combination of machine learning and atomistic simulation helps achieve more accurate and efficient in silico evaluations of the chemisorption and diffusion of gas molecules in MOFs.

Surface functionalization of metal nanoparticles (NPs), e.g., using peptides and proteins, has recently attracted a considerable attention in the field of design of therapeutics and diagnostics. The possibility of diverse functionalization allows them to selectively interact with proteins, while the metal core ensures solubility, making them tunable therapeutic agents against diseases due to mis-folding or aggregation. On the other hand, their action is limited by possible self-aggregation, which could be, however, prevented based on the full understanding of their phase diagram as a function of the environmental variables (temperature, ionic strength of the solution, concentration) and intrinsic characteristics (size, charge, amount, and type of functional groups). A common modeling strategy to study the phase behavior is to represent the NPs as spheres interacting via effective potentials implicitly accounting for the solvation effects. Their size put the NPs into the class of colloids, albeit with particularly complex interactions including both attractive and repulsive features, and a consequently complex phase diagram. In this work, we review the studies exploring the phases of these systems starting from those with only attractive or repulsive interactions, displaying a simpler disperse-clustered-aggregated transitions. The phase diagram is here interpreted focusing on the universal aspects, i.e., those dependent on the general feature of the potentials, and available data are organized in a parametric phase diagram. We then consider the potentials with competing attractive short range well and average-long-range repulsive tail, better representing the NPs. Through the proper combination of the attractive only and repulsive only potentials, we are able to interpret the appearance of novel phases, characterized by aggregates with different structural characteristics. We identify the essential parameters that stabilize the disperse phase potentially useful to optimize NP therapeutic activity and indicate how to tune the phase behavior by changing environmental conditions or the NP chemical-physical properties., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Bini, Brancolini and Tozzini.)

Artemis is an endonuclease responsible for breaking hairpin DNA strands during immune system adaptation and maturation as well as the processing of potentially toxic DNA lesions. Thus, Artemis may be an important target in the development of anticancer therapy, both for the sensitization of radiotherapy and for immunotherapy. Despite its importance, its structure has been resolved only recently, and important questions concerning the arrangement of its active center, the interaction with the DNA substrate, and the catalytic mechanism remain unanswered. In this contribution, by performing extensive molecular dynamic simulations, both classically and at the hybrid quantum mechanics/molecular mechanics level, we evidenced the stable interaction modes of Artemis with a model DNA strand. We also analyzed the catalytic cycle providing the free energy profile and key transition states for the DNA cleavage reaction.

The specific heat behavior in bulk nanomaterials (NMs) obtained by adding nanoparticles to pure suspending media has attracted a lot interest in recent years. Controversial results about NMs specific heat ( c p ) have been reported in the literature, where nanoparticles (NPs) of different sizes and materials were suspended in solid and liquid salts at different concentrations and temperatures. However, a unified picture explaining the c p variations with respect to the pure SM and the ideal NM (where NP and SM are assumed to not interact). We then use the model to interpret results of our classical molecular dynamics simulations, which we perform in the solid and liquid phases of NMs representative of three different classes, defined according to the atomic interactions at the interface. The analysis reveals nontrivial and competing effects influencing c p behavior in two-component ionic bulk nanomaterials containing NPs. The model, designed to work in the dilute regime, divides the NM in three regions: bulk suspending medium (SM), nanoparticles, and interface regions. It includes the effects of temperature, NP size, and NP concentration (mass fraction), allowing us to calculate c p variations with respect to the pure SM and the ideal NM (where NP and SM are assumed to not interact). We then use the model to interpret results of our classical molecular dynamics simulations, which we perform in the solid and liquid phases of NMs representative of three different classes, defined according to the atomic interactions at the interface. The analysis reveals nontrivial and competing effects influencing c p , such as system-dependent atomic rearrangements at the interface, vibrations of the NP as a whole and c p variations coming from the individual NP and SM specific heats. Our study contributes to the interpretation of past controversial results and helps in designing NMs with improved thermal properties, which is highly relevant for industrial applications in thermal energy storage and renewable energy production.

The chemical and physical properties of lanthanide coordination complexes can significantly change with small variations in their molecular structure. Further, in solution, coordination structures (e.g., lanthanide-ligand complexes) are dynamic. Resolving solution structures, computationally or experimentally, is challenging because structures in solution have limited spatial restrictions and are responsive to chemical or physical changes in their surroundings. To determine structures of lanthanide-ligand complexes in solution, a molecular simulation approach is presented in this chapter, which concurrently considers chemical reactions and molecular dynamics. Lanthanide ion, ligand, solvent, and anion molecules are explicitly included to identify, in atomic resolution, lanthanide coordination structures in solution. The computational protocol described is applicable to determining the molecular structure of lanthanide-ligand complexes, particularly with ligands known to bind lanthanides but whose structures have not been resolved, as well as with ligands not previously known to bind lanthanide ions. The approach in this chapter is also relevant to elucidating lanthanide coordination in more intricate structures, such as in the active site of enzymes., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Nowadays, boron nitride has attracted a great deal of attention due to its physical (chemical) properties, facile synthesis, and experimental characterization, indicating great potential for industrial application. Based on this, we develop here a theoretical study on boron nitride nanoflakes built-up from hexagonal boron nitride nanosheets exhibiting hexagonal, rectangular, and triangular shapes. In order to investigate geometry effects such as those due to the presence of armchair and zigzag edges and distinct shapes, we analyzed their properties from both classical and quantum viewpoints. Using classical molecular dynamics calculations, we show that the nanosheets preserve their structural stability at high temperatures, while DFT calculations demonstrate HOMO-LUMO energy gap variation within the theoretical energy gaps of h-BN in bulk and 2D crystals. Besides that, we have also found that boron nitride nanoflakes structures have spatially symmetrical spin densities.

Transitions between different conformational states are ubiquitous in proteins, being involved in signaling, catalysis, and other fundamental activities in cells. However, modeling those processes is extremely difficult, due to the need of efficiently exploring a vast conformational space in order to seek for the actual transition path for systems whose complexity is already high in the stable states. Here we report a strategy that simplifies this task attacking the complexity on several sides. We first apply a minimalist coarse-grained model to Calmodulin, based on an empirical force field with a partial structural bias, to explore the transition paths between the apo-closed state and the Ca-bound open state of the protein. We then select representative structures along the trajectory based on a structural clustering algorithm and build a cleaned-up trajectory with them. We finally compare this trajectory with that produced by the online tool MinActionPath, by minimizing the action integral using a harmonic network model, and with that obtained by the PROMPT morphing method, based on an optimal mass transportation-type approach including physical constraints. The comparison is performed both on the structural and energetic level, using the coarse-grained and the atomistic force fields upon reconstruction. Our analysis indicates that this method returns trajectories capable of exploring intermediate states with physical meaning, retaining a very low computational cost, which can allow systematic and extensive exploration of the multi-stable proteins transition pathways., (Copyright © 2019 Delfino, Porozov, Stepanov, Tamazian and Tozzini.)

Fluidity of lipid membranes is known to play an important role in the functioning of living organisms. The fluorescent probe Laurdan embedded in a lipid membrane is typically used to assess the fluidity state of lipid bilayers by utilizing the sensitivity of Laurdan emission to the properties of its lipid environment. In particular, Laurdan fluorescence is sensitive to gel vs liquid⁻crystalline phases of lipids, which is demonstrated in different emission of the dye in these two phases. Still, the exact mechanism of the environment effects on Laurdan emission is not understood. Herein, we utilize dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC) lipid bilayers, which at room temperature represent gel and liquid⁻crystalline phases, respectively. We simulate absorption and emission spectra of Laurdan in both DOPC and DPPC bilayers with quantum chemical and classical molecular dynamics methods. We demonstrate that Laurdan is incorporated in heterogeneous fashion in both DOPC and DPPC bilayers, and that its fluorescence depends on the details of this embedding.

Focal adhesion kinase (FAK) is required for signaling in the heart. S910 phosphorylated FAK is known to cause pathological cardiac hypertrophy. The switching of FAK between its inactive (-i), activated (-a) and hyperactive (-h) state is controlled by phosphorylation. FAK consists of three domains, namely: FERM, Kinase, and FAT joined by linkers L1 and L2. The structural basis of FAK phosphorylation and signaling to the downstream pathways is not understood. In this work, we carried out homology modeling and domain assembly of full length human iFAK and aFAK. 100 ns classical molecular dynamic simulations were performed using AMBER14 and effect of S910 phosphorylation on FAK was investigated. The iFAK model superposed on a small angel X-ray scattering (SAXS) derived model with RMSD of 1.18 Å for 590 Cα atoms. aFAK showed S910 phosphorylation site in L2 shielded by FERM. S910 phosphorylation in hFAK led to its exposure accompanied by a large conformational change and exposing the previously buried Grb2 interaction site responsible for causing cardiac hypertrophy. The models of FAK are in agreement with diverse experimental data and observed differences in biological action. Understanding the structure activity relationships of FAK in response to phosphorylation is important for its future therapeutic modulation., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Models of the hydrogenoxalate (bioxalate, charge -1) and oxalate (charge -2) anions were developed for classical molecular dynamics (CMD) simulations and parametrized against ab initio molecular dynamics (AIMD) data from our previous study (Kroutil et al. (2016) J Mol Model 22:210). The interactions of the anions with water were described using charges scaled according to the electronic continuum correction approach with rescaling of nonbonded parameters (ECCR), and those descriptions of anion interactions were found to agree well with relevant AIMD and experimental results. The models with full RESP charges showed excessively strong electrostatic interactions between the solute and water molecules, leading to an overstructured solvation shell around the anions and thus to a diffusion coefficient that was much too low. The effect of charge scaling was more evident for the oxalate dianion than for the hydrogenoxalate anion. Our work provides CMD models for ions of oxalic acid and extends previous studies that showed the importance of ECCR for modeling divalent ions and ions of organic compounds. Graphical abstract The radial distribution function between a water oxygen (Ow) and an oxygen of the oxalate dianion (Ox) significantly improved when scaled charges were applied to the anion.

The origin of the dissolution of methylammonium lead trihalide (MAPI) crystals in liquid water is clarified by finite-temperature molecular dynamics by developing a MYP-based force field (MYP1) for water-MAPI systems. A thermally activated process is found with an energy barrier of 0.36 eV consisting of a layer-by-layer degradation with generation of inorganic PbI 2 films and solvation of MA and I ions. We rationalize the effect of water on MAPI by identifying a transition from a reversible absorption and diffusion in the presence of vapor to the irreversible destruction of the crystal lattice in liquid due to a cooperative action of water molecules. A strong water-MAPI interaction is found with a binding energy of 0.41 eV/H 2 O and wetting energy of 0.23 N/m. The water vapor absorption is energetically favored (0.29 eV/H 2 O), and the infiltrated molecules can migrate within the crystal with a diffusion coefficient D = 1.7 × 10 -8 cm 2 /s and activation energy of 0.28 eV.

The most representative member of the Ras subfamily is its HRas isoform. Ras proteins being GTPases, possess an intrinsic activity to hydrolyze the GTP molecule to GDP. During the transition phases, between active and inactive states, P-loop and switch regions show maximum variations. Various hot-spot Ras mutants (G12V, A59G, Q61L etc) have been reported, that limit the protein's conformation in the permanent active state. In the present study, we aim to explore the structural dynamics of one such crucial mutant of Ras namely A59G which belongs to the conserved Switch II region of the protein. Approximately ∼15μs of Classical Molecular Dynamics (CMD) simulations have been carried out on the mutant and wild-type complexes. Further, a metadynamics simulation of 500ns was also carried out, which suggests an energy barrier of ∼9.56kcal/mol between wild-type and mutant conformation. We demonstrate the role of water molecule in maintaining the required interaction networks in the pre-hydrolysis state, its impact on A59G mutation, distinct orientation of the Gln61 residue in two conformations, disruption of crucial Gly60 and γ phosphate and the change in the Switch II region. The outcome of our study captures the pre-hydrolysis state of the HRas protein. It also establishes the fact that this mutation makes the movement of Switch II region and the conserved DXXGQ motif highly constrained, which is known to be an important requirement for hydrolysis. This suggests that the A59G mutation may decrease the rate of intrinsic hydrolysis as well as GAP-mediated hydrolysis., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Flexibility is involved in a wide range of biological processes, such as protein assembly and binding recognition. EPSP synthase is an enzyme that must undergo a large conformational change to accommodate its ligands into its binding cavity. However, although the structure of EPSP synthase has been determined, its plasticity has not been explored in depth. Therefore, in this work, we extensively examined the influence of the flexibility of Mycobacterium tuberculosis EPSP (MtEPSP) synthase on the function of this protein using classical and replica-exchange metadynamics simulations. We were able to identify five well-populated conformational clusters for MtEPSP synthase: two corresponding to open, one to ajar, and two to closed conformations. We also pinpointed three hydrophobic regions that are responsible for guiding transitions among these states. Taken together, the new findings presented here indicate how the hydrophobic regions modulate the flexibility of MtEPSP synthase, and they highlight the importance of considering these dynamic features in drug design projects employing this enzyme as a target. Graphical abstract The flexibility of EPSP synthase as a function of the pincer angles.

The present work combines experimental and theoretical studies of the collision between keV ion projectiles and clusters of pyrene, one of the simplest polycyclic aromatic hydrocarbons (PAHs). Intracluster growth processes induced by ion collisions lead to the formation of a wide range of new molecules with masses larger than that of the pyrene molecule. The efficiency of these processes is found to strongly depend on the mass and velocity of the incoming projectile. Classical molecular dynamics simulations of the entire collision process-from the ion impact (nuclear scattering) to the formation of new molecular species-reproduce the essential features of the measured molecular growth process and also yield estimates of the related absolute cross sections. More elaborate density functional tight binding calculations yield the same growth products as the classical simulations. The present results could be relevant to understand the physical chemistry of the PAH-rich upper atmosphere of Saturn's moon Titan.

A detailed understanding of the connections of fabrication and processing to structural and thermal properties of low-dimensional nanostructures is essential to design materials and devices for phononics, nanoscale thermal management, and thermoelectric applications. Silicon provides an ideal platform to study the relations between structure and heat transport since its thermal conductivity can be tuned over 2 orders of magnitude by nanostructuring. Combining realistic atomistic modeling and experiments, we unravel the origin of the thermal conductivity reduction in ultrathin suspended silicon membranes, down to a thickness of 4 nm. Heat transport is mostly controlled by surface scattering: rough layers of native oxide at surfaces limit the mean free path of thermal phonons below 100 nm. Removing the oxide layers by chemical processing allows us to tune the thermal conductivity over 1 order of magnitude. Our results guide materials design for future phononic applications, setting the length scale at which nanostructuring affects thermal phonons most effectively.

The generation of bond, angle, and torsion parameters for classical molecular dynamics force fields typically requires fitting parameters such that classical properties such as energies and gradients match precalculated quantum data for structures that scan the value of interest. We present a program, Paramfit, distributed as part of the AmberTools software package that automates and extends this fitting process, allowing for simplified parameter generation for applications ranging from single molecules to entire force fields. Paramfit implements a novel combination of a genetic and simplex algorithm to find the optimal set of parameters that replicate either quantum energy or force data. The program allows for the derivation of multiple parameters simultaneously using significantly fewer quantum calculations than previous methods, and can also fit parameters across multiple molecules with applications to force field development. Paramfit has been applied successfully to systems with a sparse number of structures, and has already proven crucial in the development of the Assisted Model Building with Energy Refinement Lipid14 force field., (© 2014 Wiley Periodicals, Inc.)

The diffusion coefficients of carbon dioxide (CO2) and ethanol (EtOH) in carbonated hydroalcoholic solutions and Champagne wines are evaluated as a function of temperature by classical molecular dynamics (MD) simulations and (13)C NMR spectroscopy measurements. The excellent agreement between theoretical and experimental diffusion coefficients suggest that ethanol is the main molecule, apart from water, responsible for the value of the CO2 diffusion coefficients in typical Champagne wines, a result that could likely be extended to most sparkling wines with alike ethanol concentrations. CO2 and EtOH hydrodynamical radii deduced from viscometry measurements by applying the Stokes-Einstein relationship are found to be mostly constant and in close agreement with MD predictions. The reliability of our approach should be of interest to physical chemists aiming to model transport phenomena in supersaturated aqueous solutions or water/alcohol mixtures.

Background: Free-standing carbon nanomembranes (CNM) with molecular thickness and macroscopic size are fascinating objects both for fundamental reasons and for applications in nanotechnology. Although being made from simple and identical precursors their internal structure is not fully known and hard to simulate due to the large system size that is necessary to draw definite conclusions., Results: We performed large-scale classical molecular dynamics investigations of biphenyl-based carbon nanomembranes. We show that one-dimensional graphene-like stripes constitute a highly symmetric quasi one-dimensional energetically favorable ground state. This state does not cross-link. Instead cross-linked structures are formed from highly excited precursors with a sufficient amount of broken phenyls., Conclusion: The internal structure of the CNM is very likely described by a disordered metastable state which is formed in the energetic initial process of electron irradiation and depends on the process of relaxation into the sheet phase.