Neutral Molecules Research Articles

Surface hopping molecular dynamics simulation of ultrafast methyl iodide photodissociation mapped by Coulomb explosion imaging.

We present a highly efficient approach to directly and reliably simulate photodissociation followed by Coulomb explosion of methyl iodide. In order to achieve statistical reliability, more than 40 000 trajectories are calculated on accurate potential energy surfaces of both the neutral molecule and the doubly charged cation. Non-adiabatic effects during photodissociation are treated using a Landau-Zener surface hopping algorithm. The simulation is performed analogous to a recent pump-probe experiment using coincident ion momentum imaging [Ziaee et al., Phys. Chem. Chem. Phys., 2023, 25, 9999-10010]. At large pump-probe delays, the simulated delay-dependent kinetic energy release signals show overall good agreement with the experiment, with two major dissociation channels leading to I(2P3/2) and I*(2P1/2) products. At short pump-probe delays, the simulated kinetic energy release differs significantly from the values obtained by a purely Coulombic approximation or a one-dimensional description of the dicationic potential energy surfaces, and shows a clear bifurcation near 12 fs, owing to non-adiabatic transitions through a conical intersection. The proposed approach is particularly suitable and efficient in simulating processes that highly rely on statistics or for identifying rare reaction channels.

Triggered integer charge transfer: energy-level alignment at an organic-2D material interface.

Weakly interacting systems such as organic molecules on monolayers of hexagonal boron nitride (h-BN) offer the possibility of single integer charge transfer leading to the formation of organic ions. Such open-shell systems exhibit unique optical and electronic properties which differ from their neutral counterparts. In this study, we used a joint experimental and theoretical approach to investigate the charge transfer of perylene-3,4,9,10-tetracarboxylic dianhydride (PTCDA) molecules on h-BN/Ni(111) by using differential reflectance spectroscopy (DRS), scanning tunneling spectroscopy (STS), and photoelectron orbital tomography (POT) measurements in combination with density functional theory (DFT) calculations. Our results show that the PTCDA monolayer consists of highly ordered organic radical anions and neutral molecules. In addition, the occurrence of the integer charge transfer is discussed based on the energy-level alignment. Since the integer charge transfer is not limited to PTCDA, we propose that the h-BN covered Ni(111) surface is a promising substrate for studying the optical and electronic properties of highly ordered organic anions.

Impact of the neutral molecule trapping on beam lifetime and beam profile

Recent investigations have shown that neutral molecules may be trapped by the electromagnetic field of a charged particle beam. This effect changes the residual gas distribution according to a characteristic temperature of the system consisting of beam and residual gas. Here, we investigate the effect of the modified gas density profile on the beam lifetime.

Probing the heterogeneity of molecular level organization of ionic liquids: a comparative study using neutral Nile red and cationic Nile blue sulfate as fluorescent probes for butyrolactam-based protic ionic liquids.

Ionic liquids (ILs) are liquid salts composed of cations and anions, known for their significant local heterogeneity at the molecular level. To understand the microheterogeneity with regard to their local polarity and local viscosity, we have used two structurally similar but chemically distinguishable fluorescent probes: Nile red (NR), a neutral molecule, and Nile blue sulfate (NBS), a charged molecule. A comparative study of the response of the two probes to the molecular level heterogeneity of ILs is expected to provide a better clarity of understanding regarding the charged polar domain and the uncharged hydrophobic domain of ILs. Towards this, we synthesized two butyrolactam-based protic ionic liquids (PILs), i.e., BTF and BTD, with the same ionic headgroup ([BT]+) and different alkyl tails ([RCOO]-), where {R = H, C11H23}. BTF has no significant hydrophobic domain, whereas BTD has a larger hydrophobic domain. Temperature-dependent fluorescence parameters such as fluorescence intensity, lifetime, and anisotropy were measured for both NR and NBS molecules. The use of a pair of structurally similar but ionically different probes enables differential estimation of parameters like the microviscosity of a domain using the fluorescence anisotropy parameter (r). The absorption and emission spectra of both probe molecules are observed to be blue shifted upon going from BTF to BTD. NR showed a significant blue shift in absorption and emission band maxima. Conversely, NBS exhibited a small wavelength shift, possibly influenced by the preferred location of their charged head group domain. Temperature-dependent rotational relaxation time (θ) of NR in BTD is smaller than that of NBS by 60-70%, indicating that stronger charge-charge interactions exist between the polar domain of BTD and NBS. Moreover, it is observed that the local viscosity of the BTF IL around both probes is similar, whereas there is a considerable difference for the BTD IL. These results are an indication that NBS being charged prefers to locate itself in the charged head group region of the IL, whereas NR being neutral tends to reside both in the hydrophobic domain and in the head group but is predominantly located in the hydrophobic domain.

Bidentate boron Lewis acids: synthesis by tin boron exchange reaction and host-guest complex formation.

Four bidentate boron Lewis acids based on the 1,8-diethynylanthracene backbone have been synthesized by a tin-boron exchange reaction with various chloroboranes, yielding the products in good to excellent yields. Complexation experiments of the host compounds with pyridine, pyrimidine and TMEDA demonstrated striking differences in terms of formation and solubility of the supramolecular adducts. The host-guest complexes were investigated by multinuclear NMR spectroscopy and structurally characterized by X-ray diffraction experiments, illustrating the adaptation of the host system upon adduct formation with different neutral guest molecules.

Heterogeneous Li coordination in solvent-in-salt electrolytes enables high Li transference numbers.

The transport properties and the underlying coordination structure of a ternary electrolyte consisting of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), 1,2-dimethoxyethane (DME), and 1,3-dioxolane (DOL) is studied over a wide concentration range, up to that of a Solvent-in-Salt (SiS) electrolyte. Among other advantages for next-generation battery applications, SiS electrolytes offer a high lithium transference number (tLi) of 0.73. We analyze the transport mechanism by electrophoretic NMR (eNMR), providing the mobilities (μi) of all species. Intriguingly, in the SiS region, the mobility of the neutral species DME exceeds the cation mobility (μDME > μLi), suggesting a heterogeneous transport mechanism, where the Li+ mobility is averaged over different species. Based on Raman spectroscopy, NMR spectroscopy and MD simulations, we derive a model for a concentration-dependent Li+ coordination environment with a heterogeneous Li+ coordination in the SiS region, where the 1st coordination shell either consists of TFSI- and DOL only, or of DME, TFSI-, and DOL. Lithium ions partially coordinated by DME migrate faster in an electric field, in contrast to lithium ions solely coordinated by anions and DOL molecules, explaining the peculiarity of the rapidly migrating neutral DME molecules. Further, DME is identified as an exclusively bidentate ligand, while TFSI- and DOL act as bridging ligands coordinating different Li+ ions. Thus, Li+ coordination heterogeneity is the basis for Li+ transport heterogeneity and for achieving very high Li+ transference numbers. In addition, an effective dynamic decoupling of Li+ and anions occurs with an Onsager coefficient σ+- ≈ 0. These results provide a deeper understanding of the very efficient lithium-ion transport in SiS electrolytes, with the potential to bring further improvements for battery applications.

Features of optoelectronic processes in a molecular junction based on a fluorophore with optically active Frontier π-orbitals: electrofluorochromism in a ZnPc-based junction.

A model of the optoelectronic process in a molecular junction has been developed, in which electron transfer occurs through transmission channels associated with the filling of the π and π* orbitals of the fluorophore with transferred electrons. The contribution of each channel to the formation of current and electroluminescence (EL) is determined by the probability of the realization of those electronic states of the molecule that, at a given bias voltage, are involved in electron transfer. It is shown that in the vicinity of critical bias voltage, stepwise changes in current and EL occur, and the height of each step is controlled by kinetic processes associated with both electron transfer and intramolecular transitions. Using the obtained analytical expressions for the relative intensities of the emission lines X0 and X+ and comparing theoretical results with experimental data on STM-induced EL in a ZnPc-based junction, we showed that the method for analyzing the behavior of current and EL near critical voltages can serve as an effective tool for understanding the physical mechanisms responsible for optoelectronic processes at the single-molecule level. The method also made it possible to obtain real values of the energy of the Frontier orbitals of the ZnPc molecule embedded between the electrodes, as well as the energies of those electronic states of the neutral and charged molecules that participate in the optoelectronic process, including electrofluorochromism.

Participation of transition metal atoms in noncovalent bonds.

The existence of halogen, chalcogen, pnicogen, and tetrel bonds as variants of noncovalent σ and π-hole bonds is now widely accepted, and many of their properties have been elucidated. The ability of the d-block transition metals to potentially act as Lewis acids in a similar capacity is examined systematically by DFT calculations. Metals examined span the entire range of the d-block from Group 3 to 12, and are selected from several rows of the periodic table. These atoms are placed in a variety of neutral MXn molecules, with X = Cl and O, and paired with a NH3 nucleophile. The resulting M⋯N bonds tend to be stronger than their p-block analogues, many of them with a substantial degree of covalency. The way in which the properties of these bonds is affected by the row and column of the periodic table from which the M atom is drawn, and the number and nature of ligands, is elucidated.

Trapping an unprecedented octacoordinated iron(II) complex with neutral bis-tetrazolylpyridyl ligands and solvent molecules.

Iron(II) can show a very rich coordination chemistry with concomitant modulation of its properties as promising functional materials. Metalation of the neutral tridentate nitrogen-donor mer-coordinating ligand 2,6-bis(2-(methyl)-2H-tetrazol-5-yl)pyridine (Me2btp) with Fe(ClO4)2·6H2O through accurate solvent polarity control enables the selective crystallization of [FeHS/LS(Me2btp)2](ClO4)2·MeCN·2.75H2O (2HS/LS·MeCN·2.75H2O) as red rods, where half of the iron(II) centres resides in the low spin (LS, S = 0) state and the other half is in the high spin (HS, S = 2) state. The red rods spontaneously convert into yellow crystals once removed from the mother liquor and exposed to air due to solvent rearrangement within the crystal packing; these new crystals can be assigned to [FeHS(Me2btp)2](ClO4)2·solvent (2HS·solvent) where all the iron(II) centres are now blocked in the HS state, as confirmed by magnetic measurements. The polarity of the crystallization solvent, together with the maintenance of the crystals within the mother liquor, are pivotal for the reactivity and interconversion of different species. Indeed, upon long standing in solution, 2HS/LS·MeCN·2.75H2O converts to another form of red crystals belonging to [FeLS(Me2btp)2][FeHS(Me2btp)(MeCN)2(H2O)](ClO4)4·MeCN (2LS·3HS·MeCN), as confirmed by single crystal X-ray diffraction data. In this co-crystal, the iron(II) in 2 resides in the LS state at all temperatures while the iron(II) in 3 is blocked in the HS state. Well-formed yellow crystals could be also isolated among the red crystals of 2HS/LS·MeCN·2.75H2O, and they could be identified as the unprecedented octacoordinated species [Fe(Me2btp)2(MeCN)(H2O)](ClO4)2·H2O (1·H2O) by single-crystal X-ray diffraction. These yellow crystals are stable in the air, but slowly convert into 2LS·3HS·MeCN if kept in the mother liquor for about one week. 1·H2O can be considered the trapped intermediate in the solid state during the conversion of [FeHS(Me2btp)2]2+ into [FeHS(Me2btp)(MeCN)2(H2O)]2+ in solution, where the two tridentate ligands in the starting species can unfold to accommodate coordinated MeCN and H2O molecules, as confirmed by theoretical calculations, and eventually one of the two Me2btp is completely replaced by the solvent.

Time-resolved photoelectron diffraction imaging of methanol photodissociation involving molecular hydrogen ejection.

Imaging ultrafast atomic and molecular hydrogen motion with femtosecond time resolution is a challenge for ultrafast spectroscopy due to the low mass and small scattering cross section of the moving neutral hydrogen atoms and molecules. Here, we propose time- and momentum-resolved photoelectron diffraction (TMR-PED) as a way to overcome limitations of existing methodologies and illustrate its performance using a prototype molecular dissociation process involving the sequential ejection of a neutral hydrogen molecule and a proton from the methanol dication. By combining state-of-the-art molecular dynamics and electron-scattering methods, we show that TMR-PED allows for direct imaging of hydrogen atoms in action. More specifically, the fingerprint of hydrogen dynamics reflects the time evolution of polarization-averaged molecular-frame photoelectron angular distributions (PA-MFPADs) as would be recorded in X-ray pump/X-ray probe experiments with few-femtosecond resolution. We present the results of two precursor experiments that support the feasibility of this approach.

Deciphering the doublet luminescence mechanism in neutral organic radicals: spin-exchange coupling, reversed-quartet mechanism, excited-state dynamics.

Neutral organic radical molecules have recently attracted considerable attention as promising luminescent and quantum-information materials. However, the presence of a radical often shortens their excited-state lifetime and results in fluorescence quenching due to enhanced intersystem crossing (EISC). Recently, an experimental report introduced an efficient luminescent radical molecule, tris(2,4,6-trichlorophenyl)methyl-carbazole-anthracene (TTM-1Cz-An). In this study, we systematically performed quantum theoretical calculations combined with the path integral approach to quantitatively calculate the excited-state dynamics processes and spectral characteristics. Our theoretical findings suggest that the sing-doublet D1 state, originating from the anthracene excited singlet state, is quickly converted to the doublet (trip-doublet) state via EISC, facilitated by a significant nonequivalence exchange interaction, with ΔJ ST = 0.174 cm-1. The formation of the quartet state (Q1, trip-quartet) was predominantly dependent on the exchange coupling 3/2J TR = 0.086 cm-1 between the triplet spin electrons of anthracene and the TTM-1Cz radical. Direct spin-orbit coupling ISC to the Q1 state was minimal due to the nearly identical spatial wavefunctions of the and Q1 levels. The effective occurrence of reverse intersystem crossing (RISC) from the Q1 to D1 state is a critical step in controlling the luminescence of TTM-1Cz-An. The calculated RISC rate k RISC, including the Herzberg-Teller effect, was 3.64 × 105 s-1 at 298 K, significantly exceeding the phosphorescence and nonradiative rates of the Q1 state, thus enabling the D1 repopulation. Subsequently, a strong electronic coupling of 37.4 meV was observed between the D1 and D2 states, along with a dense manifold of doublet states near the D1 state energy, resulting in a larger reverse internal conversion rate k RIC of 9.26 × 1010 s-1. Distributed to the D2 state, the obtained emission rate of k f = 2.98-3.18 × 107 s-1 was in quite good agreement with the experimental value of 1.28 × 107 s-1, and its temperature effect was not remarkable. Our study not only provides strong support for the experimental findings but also offers valuable insights for the molecular design of high-efficiency radical emitters.

Recognition of Bis-thiourea Tweezers to Neutral Molecules in Non-Polar Solvent

Recognition of Bis-thiourea Tweezers to Neutral Molecules in Non-Polar Solvent

A ratiometric fluorescent probe revealing the abnormality of acetylated tau by visualizing polarity in Alzheimer's disease.

Abnormal neuronal polarity leads to early deficits in Alzheimer's disease (AD) by affecting the function of axons. Precise and rapid evaluation of polarity changes is very important for the early prevention and diagnosis of AD. However, due to the limitations of existing detection methods, the mechanism related to how neuronal polarity changes in AD is unclear. Herein, we reported a ratiometric fluorescent probe characterized by neutral molecule to disclose the polarity changes in nerve cells and the brain of APP/PS1 mice. Cy7-K showed a sensitive and selective ratiometric fluorescence response to polarity. Remarkably, unlike conventional intramolecular charge transfer fluorescent probes, the fluorescence quantum yield of Cy7-K in highly polar solvents is higher than that in low polar solvents due to the transition of neutral quinones to aromatic zwitterions. Using the ratiometric fluorescence imaging, we found that beta-amyloid protein (Aβ) inhibits the expression of histone deacetylase 6, thereby increasing the amount of acetylated Tau protein (AC-Tau) and ultimately enhancing cell polarity. There was a high correlation between polarity and AC-Tau. Furthermore, Cy7-K penetrated the blood-brain barrier to image the polarity of different brain regions and confirmed that APP/PS1 mice had higher polarity than Wild-type mice. The probe Cy7-K will be a promising tool for assessing the progression of AD development by monitoring polarity.

Neutral rhodol-based dyes expressing localization in mitochondria.

Neutral rhodol-based red emitters are shown to efficiently localize in mitochondria, as demonstrated by confocal microscopy and co-localization studies. A simple model is proposed to explain the localization mechanism of neutral molecules. The model takes into account the strong coupling between the molecular dipole moment and the electric field of the inner mitochondrial membrane.

Four-membered heterocyclic molecules featuring boron and heavy group 14 elements that exhibit both σ-aromatic and π-aromatic properties: a new synthetic target.

In this study, we present a series of theoretically designed B2G14G14' molecules, featuring four-membered-ring heterocycles containing boron and heavy group 14 elements (G14 and G14' = Si, Ge, Sn, and Pb). Through the use of density functional theory (DFT), natural bond orbital (NBO) analysis, quantum theory of atoms in molecules (QTAIM), and electron localization function (ELF), our studies demonstrate a strong π single bond between the bridgehead G14 and G14' atoms, with minimal participation from a very weak G14-G14' σ bond. Additionally, the nucleus independent chemical shift (NICS), anisotropy of current-induced density (ACID), and adaptive natural density partitioning (AdNDP) analyses definitively establish the presence of both σ-aromaticity and π-aromaticity in these inorganic four-membered heterocyclic neutral molecules.

Waugh Type Polyoxometalate Enneamolybdomangate(IV) as an Acid Catalyst for Hydrolysis and as an Oxidant for Aspirin: A Kinetic Study

The oxidations of aspirin and salicylic acid by enneamolybdomanganate(IV) were studied in perchloric acid medium. The aspirin undergo hydrolysis prior to the oxidation to generate salicylic acid, which was found to be catalyzed by both HClO4 as well as the oxidant. The rate of hydrolysis of aspirin catalyzed by the oxidant was found to be lesser than that of HClO4. A mechanism based on the kinetic results is proposed for oxidation of both aspirin and salicylic acid. Initially, aspirin undergo hydrolysis to give salicylic acid, which is catalyzed by both the oxidant, which also acts as an acid catalyst and the HClO4 used. The salicylic acid generated in prior hydrolysis reacts with an oxidant undergoing oxidative decarboxylation by enneamolybdomanganate(IV) in a rate determining step generating phenol. Further oxidation of phenol by another oxidant molecule gives the final product p-benzoquinone in a fast step. The mechanism of oxidation of salicylic acid was also studied under identical conditions as that of aspirin and found to be similar except its prior hydrolysis. A plausible mechanism was proposed and the rate laws for both the substrates were derived. In order to understand the intervention of free radicals the reaction was studied in presence of added acrylonitrile. The formation of copious precipitate due to polymerization of added acrylonitrile was not observed and the values of kobs also remain constant in presence of added acrylonitrile. The formation of MnIII was also studied in presence of added tetrasodium pyrophosphate, which did not affect the values of kobs. Therefore, the formation of free radical and MnIII was excluded while proposing mechanisms for both the substrates. The effect of added products Mn2+ and molybdate also did not have any effects indicating no prior equilibria involving these two products. The effect of ionic strength and solvent polarity are in favour of the reaction between two neutral molecules. The activation parameters were determined and support an outer sphere reaction with formation of a weak precursor complex between the reactants was proposed for both aspirin and salicylic acid.

Power and Sensitivity Management of Carbon Nanotube Transistor Glucose Biosensors.

Continuous glucose monitoring (CGM), which is significant for the daily management of diabetes, requires a low-power-consumption sensor system that can track low nanomolar levels of glucose in physiological fluids, such as sweat and tears. However, traditional electrochemical methods are limited to analytes in micromolar to millimolar ranges and entail high power consumption. Carbon nanotube (CNT) film field-effect transistors (FETs) are promising for constructing extremely sensitive biosensors, but their wide applications in CGM are limited by the strong screening effect of physiological fluids and the zero charge of glucose molecules. In this study, we demonstrate a glucose aptamer-modified CNT FET biosensor to realize a highly sensitive CGM system with sub-nW power consumption by applying a suitable gate voltage. A positive gate voltage can enlarge the effective Debye screening length at the double layer to reduce the local ion population nearby and then improve the sensitivity of the FET-based biosensors by 5 times. We construct CNT FET sensors for CGM with a limit of detection of 0.5 fM, a record dynamic range up to 109, and a power consumption down to ∼100 pW. The proposed field-modulated sensing performance scheme is applicable to other aptamer-based FET biosensors for detecting neutral or less charged molecules and opens opportunities to develop facilely modulated, highly sensitive, low-power, and noninvasive CGM systems.

Improving the Accuracy of Physics-Based Hydration-Free Energy Predictions by Machine Learning the Remaining Error Relative to the Experiment.

The accuracy of computational models of water is key to atomistic simulations of biomolecules. We propose a computationally efficient way to improve the accuracy of the prediction of hydration-free energies (HFEs) of small molecules: the remaining errors of the physics-based models relative to the experiment are predicted and mitigated by machine learning (ML) as a postprocessing step. Specifically, the trained graph convolutional neural network attempts to identify the "blind spots" in the physics-based model predictions, where the complex physics of aqueous solvation is poorly accounted for, and partially corrects for them. The strategy is explored for five classical solvent models representing various accuracy/speed trade-offs, from the fast analytical generalized Born (GB) to the popular TIP3P explicit solvent model; experimental HFEs of small neutral molecules from the FreeSolv set are used for the training and testing. For all of the models, the ML correction reduces the resulting root-mean-square error relative to the experiment for HFEs of small molecules, without significant overfitting and with negligible computational overhead. For example, on the test set, the relative accuracy improvement is 47% for the fast analytical GB, making it, after the ML correction, almost as accurate as uncorrected TIP3P. For the TIP3P model, the accuracy improvement is about 39%, bringing the ML-corrected model's accuracy below the 1 kcal/mol threshold. In general, the relative benefit of the ML corrections is smaller for more accurate physics-based models, reaching the lower limit of about 20% relative accuracy gain compared with that of the physics-based treatment alone. The proposed strategy of using ML to learn the remaining error of physics-based models offers a distinct advantage over training ML alone directly on reference HFEs: it preserves the correct overall trend, even well outside of the training set.

Co-crystals: The neutral molecules in modification of crystal structure of [W(CN)6(bpy)]2− anion salts

Ten new coformers were incorporated into structure of (PPh4)2[W(CN)6(bpy)] salt (bpy = 2,2′-bipyridine), this includes such neutral molecules as: 1,10-phenanthroline, 4-hydroxyacetophenone, 2-amino-5-nitrobenzonitrile, 5-amino-2-pyridinecarbonitrile, 2-nitrobenzonitrile, 4-formylbenzoic acid, 4-pyridinecarbonitrile, 2-pyridinecarbonitrile, 2-hydroxybenzonitrile and 2-naphtol. The composition of cocrystals were confirmed by elemental analysis while the incorporation of coformers was also proved by IR spectra. The single crystal structures were solved and described for salts with 2-naphtol and 4-hydroxyacetophenone. The formation of crystals was caused by formation of strong hydrogen bonds between coformers and cyanido ligands of the anion. The structures were compared and discussed with other similar systems. Additionally, the influence of the coformer on the packing structure of hexacyanido tungsten(IV) complexes was discussed.

Multi-Scale Simulations of Gas-Phase Particles Generated in Plasma Enhanced Atomic Layer Deposition Processes

Recently, Plasma Enhanced Atomic Layer Deposition (PEALD) has been investigated as an advanced semiconductor deposition method to expand the number and variety of materials other than those available in thermal ALD. The high reactivity of the plasma species at the deposition surface during PEALD has the advantage of increasing the degree of freedom of processing conditions and enabling a wider range of material properties. PEALD consists of four steps in the order of feed, 1st purge, reaction, and 2nd purge. In this reaction step, the secondary gas-phase reaction is important because the properties of the thin film vary greatly depending on the particles that arrive on the growing surface. Therefore, the particle distribution in the gas phase before arrival on the surface should be accurately predicted to estimate the film quality. Generally, the DSMC (Direct Simulation Monte Carlo) method is a commonly used method for analyzing particle distribution in the gas phase. The DSMC method requires reaction and scattering models that simulate the behavior between molecules. Variable Hard Sphere and Variable Soft Sphere are often used as the model in DSMC. These models are about molecular transport, such as molecular diffusion, and cannot deal with molecular reactions. The total collision energy model is determined for molecular reactions, but it may not accurately describe molecular reactions in non-equilibrium flow. The purpose of the study is to develop a new reaction-scattering model in non-equilibrium flow in the DSMC method for determining the particle distribution in the gas phase in the PEALD processes. Here, we focus on ammonia plasma used in the reaction step of boron nitride nanomaterial films, which are used in various applications, such as semiconductors and fine ceramics with two-dimensional and three-dimensional structures. This study focused on neutral and ionic molecules, NH3 and NH4 +, which can be the main collisional species in secondary gas-phase reactions.A molecular dynamics simulation with a reactive force field (ReaxFF) was performed. The left figure shows the simulation system of NH3 and NH4 + collisions. The periodic boundary conditions were applied to the simulation box for all directions in the size of 100 Å × 100 Å × 100 Å. Randomly oriented NH3 and NH4 + were placed in a line at an initial distance of 15 Å. NH4 + was given an electric charge of –0.22892 on N and 0.30723 on H. The initial translational energy of NH4 + was set at 5-60 eV, and the difference in the vertical coordinate at 0.1-5.0 Å. The collisions were sampled 2000 times.The right figure shows the calculated collision cross section, which determines whether a collision occurs or not. The collision coefficient bcoll represents the difference in the maximum vertical coordinate at which the absolute value of the scattering angle is greater than 1 degree. The collision cross section decreases as the translational energy increases because the negligible attraction-force act on NH3 and NH4 + at a long distance. When the translational energy exceeds 40 eV, the collisional cross-section no longer decreases and becomes asymptotic. The reason can be considered as follows. As the cross-section approaches 40 Å2, the collision cross section asymptotically approaches 40 Å2, regardless of the translational energy, because the attractive and repulsive forces are balanced. We compared our results with a reference model in grey color. In the reference model, Etr, 0 represents the minimum translational energy, σT,0 is the collision cross section at the minimum, Etr is the translational energy, and ω is a fitting parameter. Here, ω was set at –3.1. The reference model shows a good agreement with our results. However, the collision cross section is asymptotic above 40 eV, which may affect the results when the DSMC simulations are performed. We will apply the reaction scattering model to the DSMC method and analyze the gas-phase reaction of PEALD, including a comparison of the results of existing models. Figure 1