Experimental Activation Energy Research Articles

Machine learning meets mechanistic modelling for accurate prediction of experimental activation energies.

Accurate prediction of chemical reactions in solution is challenging for current state-of-the-art approaches based on transition state modelling with density functional theory. Models based on machine learning have emerged as a promising alternative to address these problems, but these models currently lack the precision to give crucial information on the magnitude of barrier heights, influence of solvents and catalysts and extent of regio- and chemoselectivity. Here, we construct hybrid models which combine the traditional transition state modelling and machine learning to accurately predict reaction barriers. We train a Gaussian Process Regression model to reproduce high-quality experimental kinetic data for the nucleophilic aromatic substitution reaction and use it to predict barriers with a mean absolute error of 0.77 kcal mol−1 for an external test set. The model was further validated on regio- and chemoselectivity prediction on patent reaction data and achieved a competitive top-1 accuracy of 86%, despite not being trained explicitly for this task. Importantly, the model gives error bars for its predictions that can be used for risk assessment by the end user. Hybrid models emerge as the preferred alternative for accurate reaction prediction in the very common low-data situation where only 100–150 rate constants are available for a reaction class. With recent advances in deep learning for quickly predicting barriers and transition state geometries from density functional theory, we envision that hybrid models will soon become a standard alternative to complement current machine learning approaches based on ground-state physical organic descriptors or structural information such as molecular graphs or fingerprints.

Mechanism of synergistic effects and kinetics analysis in catalytic co-pyrolysis of water hyacinth and HDPE

Co-pyrolysis of water hyacinth (Eichhornia crassipes) and HDPE is a commercially effective method to mitigate the pressure of environmental pollution and the problem of invasive alien species. In this paper, the thermal behavior and kinetic evaluation of the two raw materials with or without catalysts have been investigated. The mechanism of catalytic pyrolysis was also investigated by analyzing the free radicals generated during catalytic pyrolysis. The results showed that the heating rates had a negligible effect on the pyrolysis process. The addition of HDPE could reduce the formation of char and be conducive to the complete reaction with the higher value of the comprehensive pyrolysis index (CPI). The sample containing 75% water hyacinth exhibited the strongest synergistic effect in co-pyrolysis, with the lowest value of the difference between experimental and theoretical weight loss (ΔW) and activation energy among the co-pyrolysis samples. However, the sample with 25% water hyacinth showed an inhibiting effect due to the positive value of ΔW in the main decomposition temperature. Whereas the sample containing HZSM-5 had the lowest ΔW and activation energy compared to the blends without catalysts indicating that the presence of HZSM-5 was conducive to the synergistic effect in co-pyrolysis. The promotion was then demonstrated to be due to the abundance of free radicals generated in the catalytic process. This paper aims to provide some suggestions on the solution for the alleviation of environmental and energy crisis.

First-principles study of intrinsic point defects and Xe impurities in uranium monocarbide

Based on density functional theory (DFT) calculations, we perform an extensive investigation of intrinsic point defects and Xe impurities in uranium monocarbide (UC). The DFT calculations involve both the conventional generalized gradient approximation (GGA) and the GGA+U approach with the Hubbard parametric term (U), using up to 5×5×5 supercells. GGA calculations for the formation energy of intrinsic defects demonstrate the significant effect of using larger supercells than in previous studies. Results confirm that the ⟨111⟩ and ⟨100⟩ dumbbell interstitials are the most stable interstitial configurations for U and C, respectively. The interstitial mechanisms are favored for self-diffusion of both uranium and carbon and diffusion of Xe under equilibrium conditions. Calculations also reveal that the Xe substitutional defect at the C lattice site tends to adopt an off-site configuration, which can be interpreted as a Xe interstitial–C vacancy complex. We also utilize GGA+U to assess the impact of effective U parameter (Ueff) on the results. Moreover, we introduce a method to estimate the carbon chemical potential by fitting the phase diagram composition data and propose a selection Ueff=1.25 eV based on the experimental Xe diffusion activation energy. With this approach, GGA+U calculations reproduce the available experimental data for the formation energy of the carbon Frenkel pair and can explain the stepwise recovery of intrinsic properties and burst Xe release behavior in UC observed in annealing experiments.

Adsorption and Reaction of Trimethyl and Triethyl Phosphite on Fe3O4 by Density Functional Theory

First-principle density functional theory (DFT) calculations are used to calculate the most stable structures and heats of adsorption of trimethyl and triethyl phosphites on an Fe3O4 surface in order to understand their behavior as lubricant additives. Previous surface science studies of phosphate and phosphite esters on iron oxide surface have shown that they react by sequentially desorbing the corresponding alcohol and aldehyde in equimolar amounts. This implies that the reactions are limited by the rate of alkoxide removal, followed by a rapid reaction to form the alcohol and aldehyde. It is found that the trialkyl phosphites, dialkyl phosphite, and monoalkyl phosphite all adsorb with the phosphorus atoms bound to surface Fe3+ ions, with the alkoxy groups close to parallel to the surface. The heat of adsorption increases as the alkoxide groups are removed. The postulate that the rate-limiting step in the decomposition of the phosphate esters is the sequential removal of the alkoxide group is tested by plotting the experimental activation energy obtained from temperature-programed desorption experiments versus the corresponding heats of adsorption calculated by DFT. A linear dependence shows that the reaction obeys the so-called Evans–Polanyi relation, thereby confirming the above postulate. The slope of the plot is close to unity and thus implies that the structure of the transition state of the reaction resembles the product.

Mechanisms of Solid–Gas Reactions: Reduction of Air Pollutants on Carbons

Ozone is a strong oxidizer and sulfur dioxide is a precursor to acid rain, both are air pollutants that can damage the respiratory tissues of animals and plants making them hazardous to the environment. They are isoelectronic valence O = X = O (X = S, O) molecules and their mechanism of reduction on carbons is similar. The solid–gas kinetics were studied in a flow system with a tubular reactor under differential and steady state conditions. Initial rates and product distribution were used to determine the stoichiometry of the reaction. The reduction of XO2 on carbons proceeds through a common primary mechanism with oxidized and reduced intermediates. The reactivity of the intermediates that were inserted on carbons (graphite, activated carbon, graphene oxide) modified by SO2 is selective with respect to aminolysis and thiolysis. A theoretical study of the chemisorption of SO2 on dehydrogenated pyrene as graphite active site model showed that at 900 °C the chemisorption occurs mainly on the diradical zigzag edge. Computational quantum chemistry calculations were carried out for the reduction of SO2 on graphite to produce elemental sulfur and CO2 using tetradehydrogenated-benzo[α]anthracene (TBA) as model. The mechanisms of the decarboxylation and sulfur transport steps were postulated. Ozonation of graphite showed that the 1,2,3-trioxolane decomposes to an oxyrene, eliminating O2. Both reactions, the SO2 and O3 with graphite, have the same experimental free energy of activation for the decarboxylation reaction. The results show that for SO2 the desulfurization has a much lower energetic demand than the decarboxylation route raising the important possibility of using the reaction of reduction of SO2 on carbons to reduce the acid rain, producing elemental sulfur as the main product, without increasing the greenhouse effect due to the formation of CO2.

Oxidation of Al-Co Alloys at High Temperatures.

In this work, the high temperature oxidation behavior of Al71Co29 and Al76Co24 alloys (concentration in at.%) is presented. The alloys were prepared by controlled arc-melting of Co and Al granules in high purity argon. The as-solidified alloys were found to consist of several different phases, including structurally complex m-Al13Co4 and Z-Al3Co phases. The high temperature oxidation behavior of the alloys was studied by simultaneous thermal analysis in flowing synthetic air at 773–1173 K. A protective Al2O3 scale was formed on the sample surface. A parabolic rate law was observed. The rate constants of the alloys have been found between 1.63 × 10−14 and 8.83 × 10−12 g cm−4 s−1. The experimental activation energies of oxidation are 90 and 123 kJ mol−1 for the Al71Co29 and Al76Co24 alloys, respectively. The oxidation mechanism of the Al-Co alloys is discussed and implications towards practical applications of these alloys at high temperatures are provided.

Activation energy barriers for Na migration in Na12A zeolite: The main contribution to ionic current via doubly occupied NaII site?

Abstract The activation barriers for Na jumps between the cationic sites NaI, NaII, and NaIII are calculated using VASP in the presence of anions (Cl−, Br−) or H2O as well as for the empty Na12A zeolite model. The inter-cage cationic NaIII→NaII’→NaIII′ path, where NaIII′ is located in the neighbor α–cage, was modeled via an intermediate state (NaII′) corresponding to two NaII cations in one 8R window. This NaII′ point can also serve as an intermediate for the second type of intra-cage Na drift, i.e., to a NaIII site in the same cage. Reasonable agreement with experimental activation energy was achieved for NaII’→NaIII′, NaII→NaIII, and NaIII→NaIII jumps. The lower activation energy is obtained for the NaII’→NaIII' jump than for NaII→NaIII, which indicates the possible role of the doubly occupied 8R window as an intermediate state in the Na transfer with high frequency and lower activation energy related to ionic current (inter-cage transfer) in Na12A.

The study of decomposition of 1-ethyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide by using Termogravimetry: Dissecting vaporization and decomposition of ILs

The question of thermal stability of ionic liquids (ILs) is often raised in the scope of industrial application of ILs. The most widely used approach for this purpose is Thermo Gravimetric Analysis (TGA) with the determination of decomposition “Tonset” incorrectly called Decomposition Temperature. In the present study with the help of TGA study of 1-ethyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide we have shown the Tonset or specified mass-loss rate cannot be applied as the innitial point of decomposition. We have proposed a new approach allowing dissecting the mass loss caused by vaporization of the sample and its thermal decomposition. The further analysis of the decomposition mass loss by using “isoconversional kinetics” allows to determine kinetic parameters of decomposition such as pre-exponential factor A and activation energy EA. The experimental activation energy EA = (250 ± 15) kJ·mol−1 agreed with the gas phase B3LYP calculated value by Kroon EA = 250 kJ·mol−1.

From Carbodiimides to Carbon Dioxide: Quantification of the Electrophilic Reactivities of Heteroallenes.

Kinetics of the reactions of isocyanates, isothiocyanates, carbodiimides, carbon disulfide, and carbon dioxide with carbanions or enamines (reference nucleophiles) have been measured photometrically in acetonitrile or DMSO solution at 20 °C. The resulting second-order rate constants and the previously published reactivity parameters N and sN of the reference nucleophiles were substituted into the correlation logk2(20°C) = sN(N + E) to determine the electrophilicity parameters of the heteroallenes: TsNCO (E = -7.69) ≫ PhNCO (E = -15.38) > CS2 (E = -17.70) ≈ PhNCS (E = -18.15) > PhNCNPh (E = -20.14) ≫ CyNCNCy (E ≈ -30). An approximate value could be derived for CO2 (-16 < E < - 11). Quantum chemical calculations were performed at the IEFPCM(DMSO)/B3LYP-D3/6-311+G(d,p) level of theory and compared with experimental Gibbs activation energies. The distortion-interaction model was used to rationalize the different reactivities of O- and S-substituted heteroallenes. Eventually it is demonstrated that the electrophilicity parameters determined in this work can be used as ordering principle for literature-known reactions of heteroallenes.

Stress‐life prediction of 25°C polypropylene materials based on calibration of Zhurkov fatigue life model

AbstractIn this study, to evaluate the chemical and mechanical properties of polypropylene (PP), activation‐energy and tensile tests were performed at room temperature (25°C) on pure PP and PP reinforced with glass fibre (GF). To improve the prediction accuracy of the fatigue life, three models based on the calibration of the Zhurkov model were proposed: a regression model, modified strain‐rate model and lethargy coefficient‐based model. Based on the experimental data analysis and statistical assessment results, we proposed a modified strain‐rate model that satisfies the dependency of the physical parameters and is congruent with the predicted fatigue life data. The experimental data and modified strain‐rate model were compared with the direct cyclic analysis results. The tendency of the frequency factor as a correction parameter in the modified strain‐rate model corresponded to the experimental activation energy and the increasing GF content.

Infinite dilution in doped ceria and high activation energies

Oxygen ion diffusion determines the performance of materials in energy conversion, energy storage and catalysis. For nominally pure cerium oxide, experiments measure high activation enthalpies while calculations predict low activation enthalpies. Moreover, for doped oxides, e.g. doped ceria, experiments show a high activation enthalpy for both pure ceria and for high dopant fractions, leading to a minimum in activation enthalpy for small dopant fractions. While for high dopant fractions the increase in activation enthalpy is correlated with the association of oxygen vacancies and dopant ions, which are both created by doping, the minimum in activation enthalpy is assumed in the literature to be related to the maximum in ionic conductivity at similar dopant fractions. In this study, density functional theory (DFT) calculations and Kinetic Monte Carlo (KMC) simulations are combined in order to calculate the ionic conductivity and activation enthalpy in doped oxides. We show that the experimental ionic conductivity and activation energy in nominally pure cerium oxide is dominated by impurities. We resolve the discrepancy between activation enthalpies of nominally pure oxides in experiments as opposed to calculations. This will lead to a more comprehensive understanding of the oxygen ion conductivity and its underlying atomistic mechanisms. Moreover, such a focus will be of great benefit to the future development of sustainable and efficient materials.

Master-eğriler kinetik yönteminin izotermal olmayan selüloz pirolizine uygulanması ve piroliz işleminin termodinamik analizi

Kinetic modeling of thermochemical conversion methods such as pyrolysis is one of the most challenging issues for bio-refineries. It is known that cellulose together with hemicellulose and lignin mainly affect the characteristics of biomass pyrolysis. However, there is still limited knowledge about the thermal behaviors of biopolymers that go into complex multi-phase pyrolysis reactions in the literature. Therefore, cellulose pyrolysis kinetics and thermodynamics were investigated in this study. Kinetic parameters of the pyrolysis process were calculated by a combined method of master-plots and Friedman method. Active pyrolysis of cellulose is found to occur between 263 and 455 °C. Applied Friedman method was perfectly fitted with the experimental data and activation energy of the thermochemical conversion process was found between 150.8 and 190.2 kJ/mol while the mean activation energy was calculated as 164.3 kJ/mol. The comparison of kinetic models used of solid-state thermal decomposition processes indicated that the cellulose pyrolysis mechanism is a diffusion-controlled (D3) degradation process at lower conversions (0&amp;lt;α&amp;lt;0.5) and the process can be explained by reaction-based mechanisms at higher conversion degrees.

Experimental and theoretical study of the thermal decomposition of ethyl acetate during fast pyrolysis

The thermal decomposition of ethyl acetate was experimentally studied in a newly designed fast pyrolysis set-up. The results were compared to theoretical calculations and literature values in order to proof the experimental concept. The reaction was carried out in a free fall tubular reactor with a residence time of 0.15 s. The identification and quantification of products stream composition was performed online using a GC-TCD/FID. First, an overview of the reaction rate at feed volumes of 0.25, 0.50 and 0.75 mL was obtained at reaction temperature between 400 to 600 °C in intervals of 50 °C. As a result mass transfer limitation for feeds larger than 0.5 mL were identified. For the second approach, a constant feed volume of 0.25 mL and temperatures between 420 to 550 °C were investigated. Using the experimental results, a global kinetic model is proposed for the thermal decomposition of ethyl acetate into ethylene and acetic acid through a first order unimolecular reaction. Also, theoretical calculations were performed at ωB97XD/6-311++G(d,p) level. A concerted mechanism through a six-membered transition state was identified in the reaction path. The theoretical and experimental activation energy values lie within the literature values between 193 and 213 kJ/mol.

Polymeric Electrolyte Comprising a NAFION Membrane Plasticized by Dimethyl Sulfoxide and the Transport Specifics of Alkali Metal Ions in It: Quantum-Chemical Simulation

The first theoretical study of alkali-metal ion transport in a polymeric inorganic electrolyte based on a dimethyl sulfoxide-plasticized Nafion membrane is reported. The structure and intermolecular interactions in XNafion · nDMSO (X = Li, Na, K, Rb, and Cs; n = 8 and 12) ionomers are simulated by the DFT method with hybrid density functionals B3LYP, wB97XD, and PBE taking into account periodic boundary conditions and the projector augmented wave (PAW) method in the VASP and GAUSSIAN program packages. According to the calculations, the barriers rise from 0.2 to 0.4 eV in the series Li–K but lower to 0.3–0.2 eV as the radius increases further in the series Rb, Cs. These results are quantitatively consistent with experimental conductivity activation energy data: 0.26 (Li+), 0.37–0.38 (Na+, K+), 0.27 (Rb+), and 0.20 (Cs+) eV. The conclusion is drawn about the structure and conductivity of the electrolyte depending on the nature of the cations.

Selective hydrogen production at Pt(111) investigated by Quantum Monte Carlo methods for metal catalysis

AbstractThis rapid communication gives the salient points and results of the theoretical investigation of a chemical reaction for efficient selective hydrogen production. The clean fuel produced is a sustainable energy source. Accurate methods based on quantum theory are used because the changing electronic structure is a probe that monitors reactions. The reaction between water and carbon monoxide is used industrially with metal catalysts, usually platinum. There is a considerable economic and environmental challenge underpinning this fundamental investigation where bond dissociation plays an essential role. A bond dissociation process is often the limiting step of reaction rates for industrial catalysis. Most mainstream quantum approaches fail to a greater or lesser degree in the description of this process. The present work advocates a promising alternative: the initial analysis of statistical data generated by the Quantum Monte Carlo (QMC) method demonstrated very stringent statistical accuracy for essential information on hydrogen production via the water‐gas shift reaction with platinum catalyst. The transition state structure is obtained from QMC force constants and illustrated here. It corresponds to water OH‐stretch concerted with Pt‐H bond formation, whilst the OH oxygen atom begins to interact with the CO carbon. The present QMC evaluation of the corresponding activation barrier is low: 17.0 ± 0.2 kcal/mol. It is close to the experimental apparent activation energy of 17.05 kcal/mol. This method is applicable to a wide range of similar systems.

High Temperature Dielectric Response and AC Conductivity Mechanism of (Nd, Ni) codoped BiFeO3

BiFeO3(BFO)and Bi0.9Nd0.1Fe0.9Ni0.1O3 (BNFNO) nanocrystalline samples have been synthesized through the sol-gel auto combustion process. Frequency and temperature dependent dielectric measurements have been performed in order to apprehend the influence of (Nd,Ni) co-doping on the conduction mechanism of pristine BFO.Thereal part of the dielectric constant is found to increase with the Nd and Ni doping, which is accredited to higher space charge polarization induced by the dopant ions.The dielectric constant measured as a function of frequency exhibits a dispersivebehaviour within low frequency region and attains a constant value at higher frequencies. A thermally activated relaxation peak is found at around 500K in the BNFNO samplewhich shifts towards lower temperatures at higher frequencies. The AC conductivity follows the Jonscher’s universal dynamic law, which has been utilized to evaluate the value of frequency exponent ‘s’. The decrease in the value of the ‘s’ parameter with the increase in temperature establishes that the conduction mechanism is governed by the correlated barrier hopping (CBH) model. The ac conductivity data has been utilized to evaluate the minimum hopping distance (Rmin) and the binding energies (Wm) of the electrons in BNFNO. Arrhenius equation has been used to fit the experimental data and activation energy has been evaluated at different temperatures. The Activation energy is found to increase with the increase in the temperatures, owing to the creation of more oxygen vacancies instigated by (Nd, Ni) co-doping.

Magnetic chains of Fe3 clusters in the {Fe3YO2} butterfly molecular compound.

The "butterfly" molecule [Fe3Y(μ3-O)2(CCl3COO)8(H2O)(THF)3] (in brief {Fe3YO2}) includes three Fe3+ ions which build a robust Fe3 cluster with a strong intracluster antiferromagnetic exchange and a total spin S = 5/2. It represents the starting magnetic system to study further interactions with magnetic rare earths when Y is replaced with lanthanides. We present heat capacity and equilibrium susceptibility measurements below 2 K, which show that each cluster has a sizeable magnetic anisotropy pointing to the existence of intercluster interactions. However, no phase transition to a long-range magnetically ordered phase is observed down to 20 mK. The intercluster interaction is analysed in the framework of the one-dimensional Blume-Capel model with an antiferromagnetic chain interaction constant J/kB = -40(2) mK between Fe3 cluster spins, and a uniaxial anisotropy with parameter D/kB = -0.56(3) K. This is associated to single chains of Fe3 clusters oriented along the shortest intercluster distances displayed by the crystal structure of {Fe3YO2}. Ac susceptibility measurements reveal that the magnetic relaxation is dominated by a quantum tunnelling process below 0.2 K, and by thermally activated processes above this temperature. The experimental activation energy of this single chain magnet, Ea/kB = 3.4(6) K, can be accounted for by the combination of contributions arising from single-molecule magnetic anisotropy and spin-spin correlations along the chains.

Enhancing a de novo enzyme activity by computationally-focused ultra-low-throughput screening.

Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers (∼2 × 104 M-1 s-1 and ∼102 s-1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within ∼2 kcal mol-1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stability-guidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design.

Hot air drying kinetics of novel functional carrot snack: Impregnated using polyphenolic rich osmotic solution with ultrasound pretreatment

AbstractIn this study, osmotic impregnation technique combined with ultrasonic pretreatment was applied to carrot slices followed by hot air drying in order to production of functional and healthy snack enriched with polyphenolic compounds derived from roselle extract (at least 400 ppm). The purposes of osmotic impregnation‐dehydration process were the maximum amount of dehydration and phenolic compounds diffusion along with the minimum sugar gain. Results revealed that increasing osmotic solution concentration and process time over optimal points, resulted in reduction of phenolic compounds saturation inside the samples. Also, it was found that prolonged ultrasound pretreatment had a damaging effect on fruit tissue microchannels and led to a decrease in the mass transfer rate. The optimal conditions of this stage were osmotic solution concentration of 60% and the process time of 60 min without ultrasonic pretreatment. In the second stage, the results of hot air drying showed that the highest amount in the total phenolic content was achieved at 70°C. Furthermore, among the several models evaluated, Jena &amp; Das model with a determination coefficient of 0.99 and minimum error indexes, was the best fitted to drying experimental data and activation energy for effective diffusivity coefficient was calculated as 35.987 kJ/mol.Practical applicationsDue to elevation in consumer's awareness in recent decades about the importance of food composition for healthier life and diseases prevention, demand for functional foods has been increasing. Food enrichment with bioactive compounds like polyphenol compounds from different plant sources is an applicable method for this purpose. In this study, carrot, which is a valuable vegetable itself, was further extended to a more functional food by enrichment with rosella plant extract and then hot air dried to produce a snack which has more nutritional quality than conventional ones. Also, osmotic dehydration results in lower moisture content of carrot slices which helps reduced hot air drying time and lower energy consumption. Overall, results from this research is usable to develop a functional snack based on carrot and rosella extract.

Morphology Evolution during Lithium-Based Vertically Aligned Nanocomposite Growth.

Ceramic-based nanocomposites are a rapidly evolving research area as they are currently being used in a wide range of applications. Epitaxial vertically aligned nanocomposites (VANs) offer promising advantages over conventional planar multilayers as key functionalities are tailored by the strong coupling at their vertical interfaces. However, limited knowledge exists of which material systems are compatible in composite films and which types of structures are optimal for a given functionality. No lithium-based VANs have yet been explored for energy storage, while 3D solid-state batteries offer great promise for enhanced energy and power densities. Although solid-on-solid kinetic Monte Carlo simulation (KMCS) models of VAN growth have previously been developed, phase separation was forced into the systems by limiting hopping directions and/or tuning the activation energies for hopping. Here, we study the influence of the temperature and deposition rate on the morphology evolution of lithium-based VANs, consisting of a promising LiMn2O4 cathode and a Li0.5La0.5TiO3 electrolyte, by applying a KMCS model with activation energies for hopping obtained experimentally and with minimum restrictions for hopping directions. Although the model considers only the kinetic processes away from thermodynamic equilibrium, which would determine the final shape of the pillars within the matrix, the trends in pillar size and distribution within the simulated VANs are in good agreement with experiments. This provides an elegant tool to predict the growth of VAN materials as the experimental activation energies and higher degrees of freedom for hopping result in a more realistic and low computational cost model to obtain accurate simulations of VAN materials.