Bromine Incorporation Affects Phase Transformations and Thermal Stability of Lead Halide Perovskites.

Mixed-cation and mixed-halide lead halide perovskites show great potential for their application in photovoltaics. Many of the high-performance compositions are made of cesium, formamidinium, lead, iodine, and bromine. However, incorporating bromine in iodine-rich compositions and its effects on the thermal stability of the perovskite structure has not been thoroughly studied. In this work, we study how replacing iodine with bromine in the state-of-the-art Cs0.17FA0.83PbI3 perovskite composition leads to different dynamics in the phase transformations as a function of temperature. Through a combination of structural characterization, cathodoluminescence mapping, X-ray photoelectron spectroscopy, and first-principles calculations, we reveal that the incorporation of bromine reduces the thermodynamic phase stability of the films and shifts the products of phase transformations. Our results suggest that bromine-driven vacancy formation during high temperature exposure leads to irreversible transformations into PbI2, whereas materials with only iodine go through transformations into hexagonal polytypes, such as the 4H-FAPbI3 phase. This work sheds light on the structural impacts of adding bromine on thermodynamic phase stability and provides new insights into the importance of understanding the complexity of phase transformations and secondary phases in mixed-cation and mixed-halide systems.

2D Metal/Graphene and 2D Metal/Graphene/Metal Systems for Electrocatalytic Conversion of CO2 to Formic Acid.

Efficiently transforming CO2 into renewable energy sources is crucial for decarbonization efforts. Formic acid (HCOOH) holds great promise as a hydrogen storage compound due to its high hydrogen density, non-toxicity, and stability under ambient conditions. However, the electrochemical reduction of CO2 (CO2 RR) on conventional carbon black-supported metal catalysts faces challenges such as low stability through dissolution and agglomeration, as well as suffering from high overpotentials and the necessity to overcome the competitive hydrogen evolution reaction (HER). In this study, we modify the physical/chemical properties of metal surfaces by depositing metal monolayers on graphene (M/G) to create highly active and stable electrocatalysts. Strong covalent bonding between graphene and metal is induced by the hybridization of sp and d orbitals, especially the sharp d z 2 ${{d}_{{z}^{2}}}$ , d y z ${{d}_{yz}}$ , and d x z ${{d}_{xz}}$ orbitals of metals near the Fermi level, playing a decisive role. Moreover, charge polarization on graphene in M/G enables the deposition of another thin metallic film, forming metal/graphene/metal (M/G/M) structures. Finally, evaluating overpotentials required for CO2 reduction to HCOOH, CO, and HER, we find that Pd/G, Pt/G/Ag, and Pt/G/Au exhibit excellent activity and selectivity toward HCOOH production. Our novel 2D hybrid catalyst design methodology may offer insights into enhanced electrochemical reactions through the electronic mixing of metal and other p-block elements.

Thermoresponsive Conductivity of Graphene-Based Fibers.

Smart materials are versatile material systems which exhibit a measurable response to external stimuli. Recently, smart material systems have been developed which incorporate graphene in order to share on its various advantageous properties, such as mechanical strength, electrical conductivity, and thermal conductivity as well as to achieve unique stimuli-dependent responses. Here, a graphene fiber-based smart material that exhibits reversible electrical conductivity switching at a relatively low temperature (60 °C), is reported. Using molecular dynamics (MD) simulation and density functional theory-based non-equilibrium Green's function (DFT-NEGF) approach, it is revealed that this thermo-response behavior is due to the change in configuration of amphiphilic triblock dispersant molecules occurring in the graphene fiber during heating or cooling. These conformational changes alter the total number of graphene-graphene contacts within the composite material system, and thus the electrical conductivity as well. Additionally, this graphene fiber fabrication approach uses a scalable, facile, water-based method, that makes it easy to modify material composition ratios. In all, this work represents an important step forward to enable complete functional tuning of graphene-based smart materials at the nanoscale while increasing commercialization viability.

Nonorthogonal Cascade Catalysis in Multicompartment Micelles.

Multicompartment micelles (MCMs) containing acid and base sites in discrete domains are prepared from poly(norbornene)-based amphiphilic bottlebrush copolymers in aqueous media. The acid and base sites are localized in different compartments of the micelle, enabling the nonorthogonal reaction sequence: deacetalization - Knoevenagel condensation - Michael addition of acetals to 2-amino chromene derivatives. Computational simulations using dissipative particle dynamics (DPD) elucidated the bottlebrush composition required to effectively site-isolate the nonorthogonal catalysts. This contribution presents MCMs as a new class of nanostructures for one-pot multistep nonorthogonal cascade catalysis, laying the groundwork for the isolation of three or more incompatible catalysts to synthesize value-added compounds in a single reaction vessel, in water.

The selective heating effect of microwave irradiation on a binary mixture of water and polyethylene oxide: a molecular dynamics simulation approach.

In this study, we investigate the molecular mechanisms of a microwave-driven selective heating process by performing molecular dynamics simulations for three different systems including pure water, pure polyethylene oxide (PEO), and water-PEO mixed systems in the presence of a microwave with two different intensities of electric field such as 0.001 V Å-1 and 0.01 V Å-1 at a frequency of 100 GHz. First, from performing molecular dynamics simulations of CO and CO2 in the presence of the microwave, it is confirmed that the molecular dipole moment is responsible for the rotational motion induced by the oscillating electric field. Second, by analyzing the MD simulations of the pure water system, we discover that the dipole moment of water exhibits a time lag with respect to the microwave. During the heating process, however, the temperature, kinetic, and potential energies increase synchronously with the oscillating electric field of the microwave, showing that the heating of the water system is caused by the molecular reaction of water to the microwave. Comparing the water-PEO mixed system to the pure water and pure PEO systems, the water-PEO mixed system has a higher heating rate than the pure PEO system but a lower heating rate than the pure water system. Therefore, we conclude that heating the water-PEO mixed system is driven by water molecules selectively activated by microwave irradiation. We also calculate the diffusion coefficients of water molecules and PEO chains by describing their mean square displacements, demonstrating that the diffusion coefficients are increased in the presence of microwaves for both water and PEO in pure and mixed systems. Lastly, during the microwave heating process, the structures of the water-PEO mixed system are altered as a function of the intensity of electric field, which is mainly driven by the response of water molecules.

Accelerating Solvent Selection for Type II Porous Liquids.

Type II porous liquids, comprising intrinsically porous molecules dissolved in a liquid solvent, potentially combine the adsorption properties of porous adsorbents with the handling advantages of liquids. Previously, discovery of appropriate solvents to make porous liquids had been limited to direct experimental tests. We demonstrate an efficient screening approach for this task that uses COSMO-RS calculations, predictions of solvent pKa values from a machine-learning model, and several other features and apply this approach to select solvents from a library of more than 11,000 compounds. This method is shown to give qualitative agreement with experimental observations for two molecular cages, CC13 and TG-TFB-CHEDA, identifying solvents with higher solubility for these molecules than had previously been known. Ultimately, the algorithm streamlines the downselection of suitable solvents for porous organic cages to enable more rapid discovery of Type II porous liquids.

Carbon Quantum Dot Modified Reduced Graphene Oxide Framework for Improved Alkali Metal Ion Storage Performance.

Organic materials with redox-active oxygen functional groups are of great interest as electrode materials for alkali-ion storage due to their earth-abundant constituents, structural tunability, and enhanced energy storage properties. Herein, a hybrid carbon framework consisting of reduced graphene oxide and oxygen functionalized carbon quantum dots (CQDs) is developed via the one-pot solvothermal reduction method, and a systematic study is undertaken to investigate its redox mechanism and electrochemical properties with Li-, Na-, and K-ions. Due to the incorporation of CQDs, the hybrid cathode delivers consistent improvements in charge storage performance for the alkali-ions and impressive reversible capacity (257 mAh g-1 at 50 mA g-1 ), rate capability (111 mAh g-1 at 1 A g-1 ), and cycling stability (79% retention after 10 000 cycles) with Li-ion. Furthermore, density functional theory calculations uncover the CQD structure-electrochemical reactivity trends for different alkali-ion. The results provide important insights into adopting CQD species for optimal alkali-ion storage.

Graphene Nanoribbon Hybridization of Zeolitic Imidazolate Framework Membranes for Intrinsic Molecular Separation.

Zeolitic imidazolate frameworks (ZIFs) are promising for gas separation membrane, but their molecular cut-off differs from that expected from its intrinsic aperture structure because of their flexibility. Herein, we introduced graphene nanoribbons (GNRs) to rigidify the ZIF framework. Because the sp2 edge of the GNRs induces strong anchoring effects, the modified layer can be rigidified. Particularly, when the GNRs were embedded and distributed in the ZIF-8 layer, an intrinsic aperture size of 3.4 Šwas observed, resulting in high H2 /CO2 separation (H2 permeance: 5.2×10-6  mol/m2 Pa s, ideal selectivity: 142). The performance surpasses the upper bound of polycrystalline MOF membrane performance. In addition, the membrane can be applied to blue H2 production, as demonstrated with a simulated steam reformed gas containing H2 /CO2 /CH4 . The separation performance was retained in the presence of water. The fundamentals of the molecular transport through the rigid ZIF-8 framework were revealed using molecular dynamics simulations.

DFT-Machine Learning Approach for Accurate Prediction of pKa.

In this study, we propose a novel method of pKa prediction in a diverse set of acids, which combines density functional theory (DFT) method with machine learning (ML) methods. First, the DFT method with B3LYP/6-31++G**/SM8 is used to predict pKa, yielding a mean absolute error of 1.85 pKa units. Subsequently, such pKa values predicted from the DFT method are employed as one of 10 molecular descriptors for developing ML models trained on experimental data. Kernel Ridge Regression (KRR), Gaussian Process Regression, and Artificial Neural Network are optimized using three Pipelines: Pipeline 1 involving only hyperparameter optimization (HPO), Pipeline 2 involving HPO followed by a relative contribution analysis (RCA) and recursive feature elimination (RFE), and Pipeline 3 involving HPO followed by RCA and RFE on an expanded set of composite features. Finally, it is demonstrated that KRR with Pipeline 3 yields optimal pKa prediction at an MAE of 0.60 log units. This algorithm was then utilized to predict the pKa of 37 novel acids. The two most important features were determined to be the number of hydrogen atoms in the molecule and the degree of oxidation of the acid. The predicted pKa values were documented for future reference.

Collagen intrafibrillar mineralization as a result of the balance between osmotic equilibrium and electroneutrality.

Mineralization of fibrillar collagen with biomimetic process-directing agents has enabled scientists to gain insight into the potential mechanisms involved in intrafibrillar mineralization. Here, by using polycation- and polyanion-directed intrafibrillar mineralization, we challenge the popular paradigm that electrostatic attraction is solely responsible for polyelectrolyte-directed intrafibrillar mineralization. As there is no difference when a polycationic or a polyanionic electrolyte is used to direct collagen mineralization, we argue that additional types of long-range non-electrostatic interaction are responsible for intrafibrillar mineralization. Molecular dynamics simulations of collagen structures in the presence of extrafibrillar polyelectrolytes show that the outward movement of ions and intrafibrillar water through the collagen surface occurs irrespective of the charges of polyelectrolytes, resulting in the experimentally verifiable contraction of the collagen structures. The need to balance electroneutrality and osmotic equilibrium simultaneously to establish Gibbs-Donnan equilibrium in a polyelectrolyte-directed mineralization system establishes a new model for collagen intrafibrillar mineralization that supplements existing collagen mineralization mechanisms.

Electrochemical Properties of Boron-Doped Fullerene Derivatives for Lithium-Ion Battery Applications.

The high electron affinity of fullerene C60 coupled with the rich chemistry of carbon makes it a promising material for cathode applications in lithium-ion batteries. Since boron has one electron less than carbon, the presence of boron on C60 cages is expected to generate electron deficiency in C60 , and thereby to enhance its electron affinity. By using density functional theory (DFT), we studied the redox potentials and electronic properties of C60 and C59 B. We have found that doping C60 with one boron atom results in a substantial increase in redox potential from 2.462 V to 3.709 V, which was attributed to the formation of an open shell system. We also investigated the redox and electronic properties of C59 B functionalized with various redox-active oxygen containing functional groups (OCFGs). For the combination of functionalization with OCFGs and boron doping, it is found that the enhancement of redox potential is reduced, which is mainly attributed to the open shell structure being changed to a closed-shell one. Nevertheless, the redox potentials are still higher than that of pristine C60 . From the observation that the lowest unoccupied molecular orbital of closed-shell OCFG- functionalized C59 B is correlated well with the redox potential, it was confirmed that the spin state is crucial to be considered to understand the relationship between electronic structure and redox properties.

Density Functional Theory - Machine Learning Approach to Analyze the Bandgap of Elemental Halide Perovskites and Ruddlesden-Popper Phases.

In this study, we have developed a protocol for exploring the vast chemical space of possible perovskites and screening promising candidates. Furthermore, we examined the factors that affect the band gap energies of perovskites. The Goldschmidt tolerance factor and octahedral factor, which range from 0.98 to 1 and from 0.45 to 0.7, respectively, are used to filter only highly cubic perovskites that are stable at room temperature. After removing rare or radioactively unstable elements, quantum mechanical density functional theory calculations are performed on the remaining perovskites to assess whether their electronic properties such as band structure are suitable for solar cell applications. Similar calculations are performed on the Ruddlesden-Popper phase. Furthermore, machine learning was utilized to assess the significance of input parameters affecting the band gap of the perovskites.

Molecular Modeling Approach to Determine the Flory-Huggins Interaction Parameter for Polymer Miscibility Analysis.

In this work, we present a thorough procedure for estimating the Flory-Huggins χ-parameter for use in atomistic and mesoscale molecular simulations in computational materials science. In particular, we propose improvements upon traditional Flory-Huggins theory by implementing a Connolly volume normalization (CVN). We apply this technique to several test systems, including a blend of poly (epichlorohydrin) and poly (methyl acrylate), a blend of polyethylene glycol and poly (methyl methacrylate), a blend of polystyrene and deuterated polystyrene, and three molecular-weight variants (monomer, dimer, and trimer) of a triblock copolymer for use in multicompartment micelle applications. Our results demonstrate that the newly developed procedure offers high accuracy and efficiency in predicting the Flory-Huggins χ-parameter for miscibility analysis compared to traditional experimental and computational methods. There are still several factors that cause the magnitude of the χ-parameter to vary between simulations performed on molecular species with the same identity but different degrees of polymerization; although we discuss possible explanations for these factors, this is nonetheless a primary focus for further exploration into this new methodology.

Molecular Dynamics Simulations of Aldol Condensation Catalyzed by Alkylamine-Functionalized Crystalline Silica Surfaces.

Molecular dynamics simulations are performed to investigate the cooperatively catalyzed aldol condensation between acetone and 4-nitrobenzaldehyde on alkylamine (or alkylenamine)-grafted silica surfaces, focusing on the mechanism of the catalytic activation of the acetone and 4-nitrobenzaldehyde by the acidic surface silanols followed by the nucleophilic attack of the basic amine functional group toward the activated reactant. From the analysis of the correlations between the catalytically active acid-base sites and reactants, it is concluded that the catalytic cooperativity of the acid-base pair can be affected by two factors: (1) the competition between the silanol and the amine (or enamine) to form a hydrogen bond with a reactant and (2) the flexibility of the alkylamine (or alkylenamine) backbone. Increasing the flexibility of the alkylamine facilitates the nucleophilic attack of the amine on the reactants. From the molecular dynamics simulations, it is found that C3 propylamine and C4 butylamine linkers exhibit the highest probability of reaction, which is consistent with the experimental observation that the activity of the aldol reaction on mesoporous silica depends on the length of alkylamine grafted on the silica surface. This simulation work serves as a pioneering study demonstrating how the molecular simulation approach can be successfully employed to investigate the cooperative catalytic activity of such bifunctional acid-base catalysts.

First-Principles Density Functional Theory Modeling of Li Binding: Thermodynamics and Redox Properties of Quinone Derivatives for Lithium-Ion Batteries.

The Li-binding thermodynamics and redox potentials of seven different quinone derivatives are investigated to determine their suitability as positive electrode materials for lithium-ion batteries. First, using density functional theory (DFT) calculations on the interactions between the quinone derivatives and Li atoms, we find that the Li atoms primarily bind with the carbonyl groups in the test molecules. Next, we observed that the redox properties of the quinone derivatives can be tuned in the desired direction by systematically modifying their chemical structures using electron-withdrawing functional groups. Further, DFT-based investigations of the redox potentials of the Li-bound quinone derivatives provide insights regarding the changes induced in their redox properties during the discharging process. The redox potential decreases as the number of bound Li atoms is increased. However, we found that the functionalization of the quinone derivatives with carboxylic acids can improve their redox potential as well as their charge capacity. Through this study, we also determined that the cathodic activity of quinone derivatives during the discharging process relies strongly on the solvation effect as well as on the number of carbonyl groups available for further Li binding.

Dissipative particle dynamics simulation study of poly(2-oxazoline)-based multicompartment micelle nanoreactor.

We investigate multicompartment micelles consisting of poly(2-oxazoline)-based triblock copolymers for nanoreactor applications, using the DPD simulation method to characterize the internal structure of the micelles and the distribution of reactant. The DPD simulation parameters are determined from the Flory-Huggins interaction parameter (χFH). From the snapshots of the micellar structures and radial distribution function of polymer blocks, it is clearly presented that the micelle is multicompartmental. In addition, by implementing the DPD simulations in the presence of reactants, it is found that Reac-C4 and Reac-OPh are associate well with the hydrophilic shell of the micelle, whereas the other two reactants, Reac-Ph and Reac-Cl, are not incorporated into the micelle. From our DPD simulations, we confirm that the miscibility (solubility) of reactant with the micelle has a strong correlation with the rate of hydrolysis kinetic resolution. Utilizing accurate methods evaluating accurate χFH parameters for molecular interactions in micelle system, this DPD simulation can have a great potential to predict the structures of micelles consisting of designed multiblock copolymers for useful reactions.

Thermodynamic and redox properties of graphene oxides for lithium-ion battery applications: a first principles density functional theory modeling approach.

Understanding the thermodynamic stability and redox properties of oxygen functional groups on graphene is critical to systematically design stable graphene-based positive electrode materials with high potential for lithium-ion battery applications. In this work, we study the thermodynamic and redox properties of graphene functionalized with carbonyl and hydroxyl groups, and the evolution of these properties with the number, types and distribution of functional groups by employing the density functional theory method. It is found that the redox potential of the functionalized graphene is sensitive to the types, number, and distribution of oxygen functional groups. First, the carbonyl group induces higher redox potential than the hydroxyl group. Second, more carbonyl groups would result in higher redox potential. Lastly, the locally concentrated distribution of the carbonyl group is more beneficial to have higher redox potential compared to the uniformly dispersed distribution. In contrast, the distribution of the hydroxyl group does not affect the redox potential significantly. Thermodynamic investigation demonstrates that the incorporation of carbonyl groups at the edge of graphene is a promising strategy for designing thermodynamically stable positive electrode materials with high redox potentials.

A first-principles study of lithium adsorption on a graphene-fullerene nanohybrid system.

The mechanism of Li adsorption on a graphene-fullerene (graphene-C60 ) hybrid system has been investigated using density functional theory (DFT). The adsorption energy for Li atoms on the graphene-C60 hybrid system (-2.285 eV) is found to be higher than that on bare graphene (-1.375 eV), indicating that the Li adsorption on the former system is more stable than on the latter. This is attributed to the high affinity of Li atoms to C60 and the charge redistribution that occurs after graphene is mixed with C60 . The electronic properties of the graphene-C60 system such as band structure, density of states, and charge distribution have been characterized as a function of the number of Li atoms adsorbed in comparison to those of the pure graphene and C60 . Li adsorption is found to preferentially occur on the C60 side due to the high adsorption energy of Li on C60 , which imparts a metallic character to the C60 in the graphene-C60 hybrid system.

Interaction of C60 with water: first-principles modeling and environmental implications.

The nature of fullerene-water interactions has been the subject of much research and debate. Specifically, the presence of a stabilizing, negative surface potential on colloidal aggregates of C60 in water is unexpected, given the neutral nature of pure carbon, and is not well understood. Previous simulation efforts have focused on the C60-water interaction using molecular dynamics simulations that lacked the ability to account for charge transfer and distribution interactions. In this study, first-principles density functional theory was used to analyze the fundamental electronic interactions to elucidate the polarization and charge transfer between water and C60. Simulations show that charge is inductively transferred to the C60 from water molecules, with subsequent polarization of the C60 molecule. In a case with two neighboring C60 molecules, the charge polarization induces a charge onto the second C60. Simulation suggests that this charge transfer and polarization may contribute at least partly to the observed negative surface potential of fullerene aggregates and, combined with hydrogen bonding network formation around C60, provides a fundamental driving force for aggregate formation in water.

Functionalized fullerenes in water: a closer look.

The excellent photophysical properties of C60 fullerenes have spurred much research on their application to aqueous systems for biological and environmental applications. Spontaneous aggregation of C60 in water and the consequent diminution of photoactivity present a significant challenge to aqueous applications. The mechanisms driving the reduction of photoactivity in fullerene aggregates and the effects of functionalization on these processes, however, are not well understood. Here, we take a closer look at the molecular phenomena of functionalized fullerene interactions in water utilizing simulation and experimental tools. Molecular dynamic simulations were performed to investigate time-evolved molecular interactions in systems containing fullerenes with water, oxygen, and/or neighboring fullerene molecules, complimented by physical and chemical characterizations of the fullerenes pre- and postaggregation. Aggregates with widely different photoactivities exhibit similar fullerene-water interactions as well as surface and aggregation characteristics. Photoactive fullerene aggregates had weaker fullerene-fullerene and fullerene-O2 interactions, suggesting the importance of molecular interactions in the sensitization route.