Defect Dynamics in MAPbI3 Polycrystalline Films: The Trapping Effect of Grain Boundaries

Metal halide perovskites, chemical compounds of ABX₃ stoichiometry (A=CH₃NH₃⁺, Cs⁺, …; B=Pb²⁺, Sn²⁺; X=I⁻, Br⁻), have attracted great interest as they emerged as one of the most promising class of materials for low-cost photovoltaics with over 25 % certified power conversion efficiency. An important open question for further improving their efficiency and stability is the formation and dynamics of point defects, for example, iodide vacancies and interstitials. In particular, recently it has been shown that defects strongly interact with grain boundaries, which, for example, prevents a quick restoration of initial conditions of film when kept in the dark after illumination [Phung et al., Adv. Energy Mater. 2020 , 1903735]. It has also been shown that iodide defects may accumulate at grain boundaries, where they induce carriers’ recombination [Park et al., ACS Energy Lett. 2019 , 4, 1321–1327]. In this article, we make use of molecular dynamics and ab initio simulations to follow the evolution and compute the energetics of a iodide vacancy, , and an iodide Interstitial, , interacting with Σ5/(102) grain boundaries of different termination, MAI and PbI₂. We show that the polarization charge of Σ5/(102) grain boundary associated to a prescribed termination drives the dynamics of charged defects, and . The long-range interaction of grain boundaries with charged species might induce the accumulation of point defects present in crystallites or formed under operation conditions. Moreover, the selective attraction of specific defects by a grain boundary may help splitting Frenkel pairs formed in solar cells under illumination, thus preventing the quick annihilation of defects and enhancing the effect of light in inducing degradation processes.     

The Trifluoromethyl Anion: Evidence for [K(crypt‐222)]+

[K(crypt‐222)]⁺ ( 1 ) and [K(crypt‐222)]⁺ ( 3 ) are isostructural, displaying nearly identical unit cell parameters. The two structures are similar to the extent that the previously reported [K(crypt‐222)]⁺ model can be refined against the new data for [K(crypt‐222)]⁺ , with extra electron density being observed from the fourth fluorine atom of the . In agreement with experimental observations, theoretical calculations suggest that deprotonated [K(crypt‐222)]⁺ is highly unstable even at as low as 195 K. The previously considered 1:1 CHF ₃ clathrate of deprotonated [K(crypt‐222)]⁺ (crystallographically indistinguishable from 1 ) is ruled out on the basis of all available data.     

Correlating Boron–Hydrogen Stretching Frequencies with Boron–Hydrogen Bond Lengths in Closoboranes: An Approach Using DFT Calculations

Using harmonic and anharmonic DFT calculations, we have established a general correlation between B–H stretching frequencies and B–H bond lengths valid for the closoboranes (x = 6 – 12), substituted closoboranes B₁₂H_{12 –} (with X = F, Cl, Br and n = 1 – 3 and 9 – 12) and the carboranes and , suggesting that this correlation is also applicable to other similar species. It appears that the average B–H stretching frequency observed around 2500 cm⁻¹ shift by about −100 cm⁻¹ if the average B–H bond length increases by 1 pm. In contrast to , the B–H bond in closoboranes is practically covalent and the correlation evidenced between its stretching frequency and its length proves to be similar to the one observed for the C–H bond.     

Investigating Possible Dipole-Bound States of Cyanopolyynes: the Case for the C5N− Anion Detected in Interstellar Space

We present results of quantum structure calculations aimed at demonstrating the possible existence of dipole-bound states (DBS) for the anion , a species already detected in the Interstellar medium (ISM). The positive demonstration of DBS existence using ab initio studies is an important step toward elucidating possible pathways for the formation of the more tightly bound valence bound states (VBS) in environments where free electrons from starlight ionization processes are known to be available to interact with the radical partner of the title molecule. Our current calculations show that such excited DBS states can exist in , in agreement with what we had previously found for the smallercyanopolyyne in the series: the anion. This system has a very weakly bound anion with binding energies of about 3 and 9 cm⁻¹ for the and DBS, respectively.     

Kinetics of two redox reactions in 2-butoxyethanol and water near its critical point of solution

The rate constants of two redox reactions and in the critical solution of 2-butoxyethanol and water have been measured by using the UV spectrophotometry at the initial reaction stage. It was found that the rate constants at various temperatures for two reactions were well described by the Arrhenius equation in the noncritical region. The critical slowing down effect was detected in the critical region. The critical slowing down exponents were determined to be 0.044 ± 0.004 and 0.046 ± 0.005 for reactions and , respectively. The values of the critical slowing down exponents showed that only dynamic critical slowing down effect, and no thermodynamic singularity could be observed for the two reactions.     

Simulations of a weakly conducting droplet under the influence of an alternating electric field

We investigate the electrohydrodynamics of an initially spherical droplet under the influence of an external alternating electric field by conducting axisymmetric numerical simulations using a charge-conservative volume-of-fluid based finite volume flow solver. The mean amplitude of shape oscillations of a droplet subjected to an alternating electric field for leaky dielectric fluids is similar to the steady-state deformation under an equivalent root mean squared direct electric field for all possible electrical conductivity ratio and permittivity ratio of the droplet to the surrounding fluid. In contrast, our simulations for weakly conducting media show that this equivalence between alternating and direct electric fields does not hold for . Moreover, for a range of parameters, the deformation obtained using the alternating and direct electric fields is qualitatively different, that is, for low and high , the droplet becomes prolate under alternating electric field but deforms to an oblate shape in the case of the equivalent direct electric field. A parametric study is conducted by varying the time period of the applied alternating electric field, the permittivity and the electrical conductivity ratios. It is observed that while increasing has a negligible effect on the deformation dynamics of the droplet for , it enhances the deformation of the droplet when for both alternating and direct electric fields. We believe that our results may be of immense consequence in explaining the morphological evolution of droplets in a plethora of scenarios ranging from nature to biology.     

On the normalized Laplacian of Möbius phenylene chain and its applications

Let denote a molecular graph of linear [n] phenylene with n hexagons and n squares, and let the Möbius phenylene chain be the graph obtained from the by identifying the opposite lateral edges in reversed way. Utilizing the decomposition theorem of the normalized Laplacian characteristic polynomial, we study the normalized Laplacian spectrum of , which consists of the eigenvalues of two symmetric matrices ℒ _R and ℒ _Q of order 3n. By investigating the relationship between the roots and coefficients of the characteristic polynomials of the two matrices above, we obtain an explicit closed-form formula of the multiplicative degree-Kirchhoff index as well as the number of spanning trees of . Furthermore, we determine the limited value for the quotient of the multiplicative degree-Kirchhoff index and the Gutman index of .     

The lift force on a charged sphere that translates and rotates in an electrolyte

The distortion of the charge cloud around a uniformly charged, dielectric, rigid sphere that translates and rotates in an unbounded binary, symmetric electrolyte at zero Reynolds number is examined. The zeta potential of the particle ζ is assumed small relative to the thermal voltage scale. It is assumed that the equilibrium structure of the cloud is slightly distorted, which requires that the Péclet numbers characterizing distortion due to particle translation, , and rotation, , are small compared to unity. Here, a is radius of the particle; D is the ionic diffusion coefficient; and , where U and Ω are the rectilinear and angular velocities of the particle, respectively. Perturbation expansions for small and are employed to calculate the nonequilibrium structure of the cloud, whence the force and torque on the particle are determined. In particular, we predict that the sphere experiences a force orthogonal to its directions of translation and rotation. This “lift” force arises from the nonlinear distortion of the cloud under the combined actions of particle translation and rotation. The lift force is given by . Here, ε is the permittivity of the electrolyte; is the Debye length; and is a negative function that decreases in magnitude with increasing . The lift force implies that an unconstrained particle would follow a curved path; an electrokinetic analog of the inertial Magnus effect. Finally, the implication of the lift force on cross‐streamline migration of an electrophoretic particle in shear flow is discussed.     

Size Evolution of Photoabsorption Spectra of Small Clusters: A Computational Study

Photoabsorption spectra of clusters, N=5–9, have been calculated using a diatomics-in-molecules like electronic structure model and a path-integral Monte Carlo sampling method. A qualitative change in the calculated spectra has been observed at N=9, which has been interpreted in terms of a structural transformation in the clusters consisting in a transition from trimer-like ionic cores observed for N≤7 to dimer-like ionic cores prevailing in through an intermediate state (comparable abundances of both types of ionic cores) observed in . The calculated spectra have been thoroughly compared with an earlier calculation on , , and reported from our group and data available for the same cluster sizes from an experiment.     

Cover Feature: How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings? (ChemPhysChem 9/2023)

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, , and to the zz component of the off-nucleus shielding tensor, , around benzene (C₆H₆) and cyclobutadiene (C₄H₄) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between and around the aromatic C₆H₆ and the antiaromatic C₄H₄, the diamagnetic and paramagnetic contributions to these quantities, and , and and , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C₆H₆ and C₄H₄ are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.     

The Threshold Photoelectron Spectrum of Fulvenone: A Reactive Ketene Derivative in Lignin Valorization

Unveiling reaction mechanisms by isomer-selective detection of reactive intermediates requires advanced spectroscopic knowledge. We study the photoionization of fulvenone (c-C₅H₄=C=O), a reactive ketene species relevant in catalytic pyrolysis of lignin, which was generated by pyrolysis of 2-methoxy acetophenone. The high-resolution threshold photoelectron spectrum (TPES) with vacuum ultraviolet synchrotron radiation revealed well-resolved vibrational transitions, assigned to ring deformation modes of the cyclopentadiene moiety. The adiabatic ionization energy was determined to be 8.25±0.01 eV and is assigned to the ²A₂← ¹A₁ transition. A broad and featureless band arising at 9 eV is associated with the ²B₁← ¹A₁ excitation. A conical intersection is responsible for the ultrafast relaxation of the fulvenone cation from the into the state resulting in a featureless and lifetime broadened band. These insights will increase the detection capabilities for fulvenone and thereby help to elucidate reaction mechanisms in lignin catalytic pyrolysis.     

On normalized Laplacians,multiplicative degree-Kirchhoff indices,and spanning trees of the linear [n]phenylenes and their dicyclobutadieno derivatives

Let be the molecular graph of the linear [n] phenylene with n hexagons and n − 1 squares, and let be the graph obtained by attaching four-membered rings to the terminal hexagons of . In this article, the normalized Laplacian spectrum of consisting of the eigenvalues of two symmetric tridiagonal matrices of order 3n is determined. An explicit closed-form formula of the multiplicative degree-Kirchhoff index (respectively the number of spanning trees) of is derived. Similarly, explicit closed-form formulas of the multiplicative degree-Kirchhoff index and the number of spanning trees of are obtained. It is interesting to see that the multiplicative degree-Kirchhoff index of (respectively ) is approximately to one half of its Gutman index.     

An Artificial Single Molecular Channel Showing High Chloride Transport Selectivity and pH-Responsive Conductance

Inspired by the unique structure and function of the natural chloride channel (ClC) selectivity filter, we present herein the design of a ClC-type single channel molecule. This channel displays high ion transport activity with half-maximal effective concentration, EC₅₀, of 0.10 μM, or 0.075 mol % (channel molecule to lipid ratio), as determined by fluorescent analysis using lucigenin-encapsulated vesicles. Planar bilayer lipid membrane conductance measurements indicated an excellent Cl⁻/K⁺ selectivity with a permeability ratio P /P up to 12.31, which is comparable with the chloride selectivity of natural ClC proteins. Moreover, high anion/anion selectivity (P /P =66.21) and pH-dependent conductance and ion selectivity of the channel molecule were revealed. The ClC-like transport behavior is contributed by the cooperation of hydrogen bonding and anion–π interactions in the central macrocyclic skeleton, and by the existence of pH-responsive terminal phenylalanine residues.     

Potential of AMnO3 (A=Ca,Sr, Ba,La) as Active Layer in Inorganic Perovskite Solar Cells

Inorganic perovskite CaMnO was proposed as a substitution for the TiO anatase in electron transport layers of solar cells containing the hybrid perovskite CH NH PbI based on increased mobility of electrons and better optical matching. Due to a suitable band gap concerning the absorption of sunlight, we investigate the potential of CaMnO and similar manganite perovskites, where Ca is replaced by either Sr, Ba or La, as an absorber layer in inorganic perovskite solar cells. In this study, we have used optical measurements on the synthesized AMnO (A=Ca, Sr, Ba, La) samples to aid density functional theory calculations (DFT) in order to accurately simulate the electronic and optical properties of AMnO compounds and gauge their potential for the role of absorber layer. Both experimental measurements and theoretical calculations show suitable band gap of 1.1-1.5 eV, depending on the compound, and absorption coefficients of the order of cm in the visible part of the spectrum.     

NMR Self-Diffusion and Transverse Relaxation Time in Bitumen: The Effect of Aging

In this investigation the dynamics of two types of bitumens with different penetration grade were tested by using dynamic shear rheometry (DSR) and Nuclear Magnetic Resonance (NMR) at unaged conditions, and upon both short- and long-term artificial aging. The gel-sol transition temperature was found to increase with increasing the time of aging treatment. Arrhenius parameters of the viscosity were found, unexpectedly, to be correlated with those of simple liquids, suggesting that the two kinds of systems, although chemically and physically quite different, share the same basic process at the molecular level. The molecular dynamics has been then investigated by NMR Pulsed Field Gradient Stimulated-Echo (PFGSE) and relaxometry (Carr-Purcell-Meiboom-Gill, CPMG, spin-echo pulse sequence) to capture the effect of aging upon dynamics variables such as self-diffusion coefficients D and transverse relaxation times T₂. The translational diffusion at T> of the light molecular components of both types of bitumens was characterized by broad distributions of D which were found independent of the experimental time scale up to 0.2 s. Similarly, T₂ data could be described as a continuous unimodal distributions of relaxation times determined both at T< and T> .     

Real-Space Interpretation of Interatomic Charge Transfer and Electron Exchange Effects by Combining Static and Kinetic Potentials and Associated Vector Fields**

Intricate behaviour of one-electron potentials from the Euler equation for electron density and corresponding gradient force fields in crystals was studied. Channels of locally enhanced kinetic potential and corresponding saddle Lagrange points were found between chemically bonded atoms. Superposition of electrostatic and kinetic potentials and electron density allowed partitioning any molecules and crystals into atomic - and potential-based -basins; -basins explicitly account for the electron exchange effect, which is missed for -ones. Phenomena of interatomic charge transfer and related electron exchange were explained in terms of space gaps between zero-flux surfaces of - and -basins. The gap between - and -basins represents the charge transfer, while the gap between - and -basins is a real-space manifestation of sharing the transferred electrons caused by the static exchange and kinetic effects as a response against the electron transfer. The regularity describing relative positions of -, -, and - basin boundaries between interacting atoms was proposed. The position of -boundary between - and -ones within an electron occupier atom determines the extent of transferred electron sharing. The stronger an H⋅⋅⋅O hydrogen bond is, the deeper hydrogen atom's -basin penetrates oxygen atom's -basin, while for covalent bonds a -boundary closely approaches a -one indicating almost complete sharing of the transferred electrons. In the case of ionic bonds, the same region corresponds to electron pairing within the -basin of an electron occupier atom.     

On the Mechanical Properties of Popgraphene-Based Nanotubes: a Reactive Molecular Dynamics Study

Carbon-based tubular materials have sparked a great interest in future electronics and optoelectronics device applications. In this work, we computationally studied the mechanical properties of nanotubes generated from popgraphene (PopNTs). Popgraphene is a 2D carbon allotrope composed of 5-8-5 rings. We carried out fully atomistic reactive (ReaxFF) molecular dynamics for PopNTs of different chiralities ( and ) and/or diameters and at different temperatures (from 300 up to 1200 K). Results showed that the tubes are thermally stable (at least up to 1200 K). All tubes presented stress/strain curves with a quasi-linear behavior followed by an abrupt drop of stress values. Interestingly, armchair-like PopNTs ( ) can stand a higher strain load before fracturing when contrasted to the zigzag-like ones ( ). Moreover, it was obtained that Young's modulus (Y_Mod) (750–900 GPa) and ultimate strength (σ_US) (120–150 GPa) values are similar to the ones reported for conventional armchair and zigzag carbon nanotubes. Y_Mod values obtained for PopNTs are not significantly temperature-dependent. While the σ_US values for the showed a quasi-linear dependence with the temperature, the exhibited no clear trends.     

Fluorine Substitution of TCNQ Alters the Redox-Driven Catalytic Pathway for the Ferricyanide-Thiosulfate Reaction

Mechanistic variation in catalysis through substituent-based redox tuning is well established. Fluorination of TCNQ (TCNQ=tetracyanoquinodimethane) provides ~850 mV variation in the redox potentials of the and (n=0, 2, 4) processes. With , catalysis of the kinetically very slow ferrocyanide-thiosulfate redox reaction in aqueous solution occurs via a mechanism in which the catalyst is reduced to when reacting with which is oxidised to . Subsequently, reacts with to form and reform the catalyst, in another thermodynamically favoured process. An analogous mechanism applies with as a catalyst. In contrast, since the reaction of with is thermodynamically unfavourable, an alternative mechanism is required to explain the catalytic activity observed in this non-fluorinated system. Here, upon addition of , reduction of to occurs with concomitant oxidation of to , which then acts as the catalyst for oxidation. Thermodynamic data explain the observed differences in the catalytic mechanisms. (n=0, 4) also act as catalysts for the ferricyanide-thiosulfate reaction in aqueous solution. The present study shows that homogeneous pathways are available following addition of these dissolved materials. Previously, these (n=0, 4) coordination polymers have been regarded as insoluble in water and proposed as heterogeneous catalysts for the ferricyanide-thiosulfate reaction. Details and mechanistic differences were established using UV-visible spectrophotometry and cyclic voltammetry.     

Insights into the Kramers’ flux‐over‐population rate for chemical reactions in liquid phases through the matrix transport equation

Kramers’ equation models a chemical reaction as a Brownian particle diffusing over a potential barrier under the influence of medium viscosity. In the case of high viscosity, the equation reduces to a simpler Smoluchowski equation. In this report, we have contrived an equivalent matrix‐transport equation that relates the ordered pair (activity, flux) of the output (activated complex) to that of the input (reactant). With an initial condition of the Dirac delta type placed at the location of the reactant, and a reflecting boundary condition set on the reactant state, and an absorbing boundary condition on the activated complex state, we are able to prove the equality relation between the mean first passage time, , for the diffusion and the inverse of the rate constant, k^?1, for the reaction counterpart. We have also derived , where λ_i is the ith eigenvalue of the Smoluchowski differential operator stipulated with the above‐mentioned boundary conditions. We have also deduced that, in the long time limit, the number of particles remaining inside the diffusion domain decays exponentially with a relaxation time just the same as the concentration of the reactant does for a first‐order reaction system.