Enhancement of atmospheric nucleation precursors on formic sulfuric anhydride induced nucleation: Theoretical mechanism

As an intermediate formed by H2SO4 (SA), formic sulfate anhydride (FSA) has been hypothesized to play a role in the nucleation of atmospheric aerosols. It is the first time that the clusters (SA)x(A)y(W)n and (FSA)x(A)y(W)n (x = 1-2; y = 1-2; n = 0-4) were systematically studied in theory on the structures, thermodynamics, intermolecular interactions, humidity dependence, atmospheric dependence and optical properties. FSA is predicted to be more stronger to promote the clustering with ammonia (A) than SA, suggesting that substituent group enhances nucleation capability of FSA. Whereas, the substituent group does not influence the humidity sensitivity of hydrated clusters. The clusters trend to form small hydrated clusters (nwater≦3). The study on atmospheric dependence indicates that the stability of the clusters depends more on temperature other than pressure. Moreover, FSA shows a stronger ability on reducing atmospheric visibility than A, SA and water molecules. This finding aims to draw attention to FSA about atmospheric nucleation.

Reactions with Criegee intermediates are the dominant gas-phase sink for formyl fluoride in the atmosphere

Atmospheric oxidation processes are of central importance in atmospheric climate models. It is often considered that volatile organic molecules are mainly removed by hydroxyl radical; however, the kinetics of some reactions of hydroxyl radical with volatile organic molecules are slow. Here we report rate constants for rapid reactions of formyl fluoride with Criegee intermediates. These rate constants are calculated by dual-level multistructural canonical variational transition state theory with small-curvature tunneling (DL-MS-CVT/SCT). The treatment contains beyond-CCSD(T) electronic structure calculations for transition state theory, and it employs validated density functional input for multistructural canonical variational transition state theory with small-curvature tunneling and for variable-reaction-coordinate variational transition state theory. We find that the M11-L density functional has higher accuracy than CCSD(T)/CBS for the HC(O)F + CH2OO and HC(O)F + anti-CH3CHOO reactions. We find significant negative temperature dependence in the ratios of the rate constants for HC(O)F + CH2OO/anti-CH3CHOO to the rate constant for HC(O)F + OH. We also find that different Criegee intermediates have different rate-determining-steps in their reactions with formyl fluoride, and we find that the dominant gas-phase removal mechanism for HC(O)F in the atmosphere is the reaction with CH2OO and/or anti-CH3CHOO Criegee intermediates.

Quantitative kinetics reveal that reactions of HO2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation

Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.

Unimolecular Reactions of E-Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH₂OH)CHOO], for its reactions with H₂O and (H₂O)₂, and for the reaction of the E-(CH₂OH)CHOO…H₂O complex with H₂O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH₂OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH₂OH)CHOO…H₂O + H₂O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H₂O)₂. In addition, we find that the atmospheric lifetime of E-(CH₂OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10^-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH₂OH)CHOO enhances its reactivity.

Synergistic Effects of SO2 and NH3 Coexistence on SOA Formation from Gasoline Evaporative Emissions

Vehicular evaporative emissions make an increasing contribution to anthropogenic sources of volatile organic compounds (VOCs), thus contributing to secondary organic aerosol (SOA) formation. However, few studies have been conducted on SOA formation from vehicle evaporative VOCs under complex pollution conditions with the coexistence of NO_x, SO₂, and NH₃. In this study, the synergistic effects of SO₂ and NH₃ on SOA formation from gasoline evaporative VOCs with NO_x were examined using a 30 m³ smog chamber with the aid of a series of mass spectrometers. Compared with the systems involving SO₂ or NH₃ alone, SO₂ and NH₃ coexistence had a greater promotion effect on SOA formation, which was larger than the cumulative effect of the two promotions alone. Meanwhile, contrasting effects of SO₂ on the oxidation state (OSc) of SOA in the presence or absence of NH₃ were observed, and SO₂ could further increase the OSc with the coexistence of NH₃. The latter was attributed to the synergistic effects of SO₂ and NH₃ coexistence on SOA formation, wherein N-S-O adducts can be formed from the reaction of SO₂ with N-heterocycles generated in the presence of NH₃. Our study contributes to the understanding of SOA formation from vehicle evaporative VOCs under highly complex pollution conditions and its atmospheric implications.

Formation mechanism of typical aromatic sulfuric anhydrides and their potential role in atmospheric nucleation process

Sulfuric anhydrides, generated from the cycloaddition reaction of SO₃ with carboxylic acids, have been revealed to be potential participants in the nucleation process of new particle formation (NPF). Hence the reaction mechanisms of typical aromatic acids (benzoic acid (BA), phenylacetic acid (PAA), phthalic acid (PA), isophthalic acid (mPA), and terephthalic acid (PTA)) with SO₃ to generate the corresponding aromatic sulfuric anhydrides were investigated by density functional theory calculations at the level of M06-2X/6-311++G(3df,3pd). As a result, these reactions were found to be feasible in the gas phase with barriers of 0.34, 0.30, 0.18, 0.08 and 0.12 kcal/mol to generate corresponding aromatic sulfuric anhydrides, respectively. The thermodynamic stabilities of clusters containing aromatic sulfuric anhydrides and atmospheric nucleation precursors (sulfuric acid, ammonia and dimethylamine) were further analyzed to identify the potential role of aromatic sulfuric anhydrides in NPF. As the thermodynamic stability of a cluster depends on both the number and strength of hydrogen bonds, the greater stability of the interactions between atmospheric nucleation precursors and aromatic sulfuric anhydrides than with aromatic acids make aromatic sulfuric anhydrides potential participators in the nucleation process of NPF. Moreover, compared with BA, the addition of a -CH₂- functional group in PAA has little influence on the reaction barrier with SO₃ but an inhibitive effect on the thermodynamic stability of clusters. The position of the two -COOH functional groups in PA, mPA and PTA does not have a consistent impact on the reaction barrier with SO₃ or the thermodynamic stability.

New Reactions for the Formation of Organic Nitrate in the Atmosphere

Organic nitrates make an important contribution to the formation of secondary organic aerosols, but the formation mechanisms of organic nitrates are not fully understood at the molecular level. In the present work, we explore a new route for the formation of organic nitrates in the reaction of formaldehyde (HCHO) with nitric acid (HNO₃) catalyzed by water (H₂O), ammonia (NH₃), and dimethylamine ((CH₃)₂NH) using theoretical methods. The present results using CCSD(T)-F12a/cc-pVTZ-F12//M06-2X/MG3S unravel that dimethylamine has a stronger catalytic ability in the reaction of HCHO with HNO₃, reducing the barrier by 21.97 kcal/mol, while water and ammonia only decrease the energy barrier by 7.35 and 13.56 kcal/mol, respectively. In addition, the calculated kinetics combined with the corresponding concentrations of these species show that the HCHO + HNO₃ + (CH₃)₂NH reaction can compete well with the naked HCHO + HNO₃ reaction at 200-240 K, which may make certain contributions to the formation of organic nitrates under some atmospheric conditions.

Quantitative Kinetics of HO2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important

Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO₂ with HCHO, CH₃CHO, and CF₃CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH₃, CF₃), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO₂ + HCHO reaction has a significant dependence on pressure, but the HO₂ + CF₃CHO reaction is almost independent of pressure. Our findings show that the HO₂ + HCHO reaction makes important contribution to the sink for HCHO, and the HO₂ + CF₃CHO reaction is the dominant sink for CF₃CHO in the atmosphere.

Temperature-dependent kinetics of the atmospheric reaction between CH2OO and acetone

Criegee intermediates are important oxidants produced in the ozonolysis of alkenes in the atmosphere. Quantitative kinetics of the reactions of Criegee intermediates are required for atmospheric modeling. However, the experimental studies do not cover the full relevant range of temperature and pressure. Here we report the quantitative kinetics of CH₂OO + CH₃C(O)CH₃ by using our recently developed dual strategy that combines coupled cluster theory with high excitation levels for conventional transition state theory and well validated levels of density functional theory for direct dynamics calculations using canonical variational transition theory including tunneling. We find that the W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ electronic structure method can be used to obtain quantitative kinetics of the CH₂OO + CH₃C(O)CH₃ reaction. Whereas previous investigations considered a one-step mechanistic pathway, we find that the CH₂OO + CH₃C(O)CH₃ reaction occurs in a stepwise manner. This has implications for the modeling of Criegee-intermediate reactions with other ketones and with aldehydes. In the kinetics calculations, we show that recrossing effects of the conventional transition state are negligible for determining the rate constant of CH₂OO + CH₃C(O)CH₃. The present findings reveal that the rate ratio between CH₂OO + CH₃C(O)CH₃ and OH + CH₃C(O)CH₃ has a significant negative dependence on temperature such that the CH₂OO + CH₃C(O)CH₃ reaction can contribute as a significant sink for atmospheric CH₃C(O)CH₃ at low temperature. The present findings should have broad implications in understanding the reactions of Criegee intermediates with carbonyl compounds and ketones in the atmosphere.

A review of secondary organic aerosols formation focusing on organosulfates and organic nitrates

Secondary organic aerosols (SOA) are crucial constitution of fine particulate matter (PM), which are mainly derived from photochemical oxidation products of primary organic matter and volatile organic compounds (VOCs), and can induce terrible impacts to human health, air quality and climate change. As we know, organosulfates (OSs) and organic nitrates (ON) are important contributors for SOA formation, which could be possibly produced through various pathways, resulting in extremely complex formation mechanism of SOA. Although plenty of research has been focused on the origins, spatial distribution and formation mechanisms of SOA, a comprehensive and systematic understanding of SOA formation in the atmosphere remains to be detailed explored, especially the most important OSs and ON dedications. Thus, in this review, we systematically summarize the recent research about origins and formation mechanisms of OSs and ON, and especially focus on their contribution to SOA, so as to have a clearer understanding of the origin, spatial distribution and formation principle of SOA. Importantly, we interpret the complex interaction with coexistence effect of SO_x and NO_x on SOA formation, and emphasize the future insights for SOA research to expect a more comprehensive theory and practice to alleviate SOA burden.

Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure

Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.

Possible atmospheric source of NH2SO3H: the hydrolysis of HNSO2 in the presence of neutral, basic, and acidic catalysts

NH₂SO₃H can directly participate in H₂SO₄-(CH₃)₂NH-based cluster formation, and thereby substantially enhance the cluster formation rate. Herein, the reaction mechanisms and kinetics for the formation of NH₂SO₃H from the hydrolysis of HNSO₂ without and with neutral (H₂O, (H₂O)₂, and (H₂O)₃), basic (NH₃ and CH₃NH₂), and acidic (HCOOH, H₂SO₄, H₂SO₄⋯H₂O, and (H₂SO₄)₂) catalysts were studied theoretically at the CCSD(T)-F12/cc-pVDZ-F12//M06-2X/6-311+G(2df,2pd) level. The calculated results showed that neutral, basic, and acidic catalysts decrease the energy barrier by over 18.1 kcal mol^-1; meanwhile, the product formation of NH₂SO₃H was more strongly bonded to neutral, basic, and acidic catalysts than to the reactants HNSO₂ and H₂O. This reveals that the reported neutral, basic, and acidic catalysts promote the formation of NH₂SO₃H from the hydrolysis of HNSO₂ both kinetically and thermodynamically. Kinetic calculations using the master equation showed that (H₂O)₂ (100% RH) dominate over the other catalysts within the range of 0-10 km altitudes and 230-320 K with its rate ratio larger by at least 2.98 times, whereas HCOOH (3.2 × 10⁹ molecules cm^-3) is the most favorable catalysts at 15 km altitude in the troposphere. Overall, the present results will provide a definitive example that neutral, basic, and acidic catalysts have important influences on atmospheric reactions.

Study of the Reactivity of CH3COOH+• and COOH+ Ions with CH3NH2: Evidence of the Formation of New Peptide-like C(O)-N Bonds

Acetamide, a small organic compound containing a peptide bond, was observed in the interstellar medium, but reaction pathways leading to the formation of this prebiotic molecule remain uncertain. We investigated the possible formation of a peptide-like bond from the reaction between acetic acid (CH₃-COOH) and methylamine (CH₃-NH₂) that were identified in the interstellar medium. From an experimental point of view, a quadrupole/octopole/quadrupole mass spectrometer was used in combination with synchrotron radiation as a tunable source of VUV photons for monitoring the reactivity of selected ions. Acetic acid was photoionized, and the reactivity of CH₃COOH^+• as well as COOH⁺ (produced from either acetic acid or formic acid) ions with neutral CH₃NH₂ was further studied. With no surprise, charge transfer, proton transfer, and concomitant dissociation processes were found to largely dominate the reactivity. However, a C(O)-N bond formation process between the two reactants was also evidenced, with a weak cross section reaction. From a theoretical point of view, results concerning reactivity and barrier heights were obtained using density functional theory, with the LC-ωPBE range-separated functional in combination with the 6-311++G(d,p) Pople basis set and are in perfect agreement with the experimental data.

Atmospheric Ring-Closure and Dehydration Reactions of 1,4-Hydroxycarbonyls in the Gas Phase: The Impact of Catalysts

1,4-Hydroxycarbonyls can potentially undergo sequential reactions involving cyclization followed by dehydration to form dihydrofurans. As dihydrofurans contain a double bond, they are highly reactive toward atmospheric oxidants such as OH, O₃, and NO₃. In the present study, we use ab initio calculations to examine the impact of various atmospheric catalysts on the energetics and kinetics of the gas-phase cyclization and dehydration reaction steps associated with 4-hydroxybutanal, a prototypical 1,4-hydroxycarbonyl molecule. The cyclization step transforms 4-hydroxybutanal into 2-hydroxytetrahydrofuran, which can subsequently undergo dehydration to form 2,3-dihydrofuran. As the barriers associated with the cyclization and dehydration steps for 4-hydroxybutanal are, respectively, 34.8 and 63.0 kcal/mol in the absence of a catalyst, both reaction steps are inaccessible under atmospheric conditions in the gas phase. However, the presence of a suitable catalyst can significantly reduce the reaction barriers, and we have examined the impact of a single molecule of H₂O, HO₂ radical, HC(O)OH, HNO₃, and H₂SO₄ on these reactions. We find that H₂SO₄ reduces the reaction barriers the greatest, with the barrier for the cyclization step being reduced to -13.1 kcal/mol and that for the dehydration step going down to 9.2 kcal/mol, measured relative to their respective separated starting reactants. Interestingly, our kinetic study shows that HNO₃ gives the fastest rate due to the combined effects of a larger atmospheric concentration and a reduced barrier. Thus, our study suggests that, with acid catalysis, the cyclization reaction step can readily occur for 1,4-hydroxycarbonyls in the gas phase. Because the dehydration step exhibits a significant barrier even with acid catalysis, the 2-hydroxytetrahydrofuran products, once formed, are likely lost through their reaction with OH radicals in the atmosphere. We have investigated the reaction pathways and the rate constant for this bimolecular reaction in the presence of excess molecular oxygen (³O₂), as it would occur under tropospheric conditions, using computational chemistry over the 200-300 K temperature range. We find that the main products from these OH-initiated oxidation reactions are succinaldehyde + HO₂ and 2,3-dihydro-2-furanol + HO₂.

Atmospheric Kinetics: Bimolecular Reactions of Carbonyl Oxide by a Triple-Level Strategy

Criegee intermediates in the atmosphere serve as oxidizing agents to initiate aerosol formation, which are particularly important for atmospheric modeling, and understanding their kinetics is one of the current outstanding challenges in climate change modeling. Because experimental kinetics are still limited, we must rely on theory for the complete picture, but obtaining absolute rates from theory is a formidable task. Here, we report the bimolecular reaction kinetics of carbonyl oxide with ammonia, hydrogen sulfide, formaldehyde, and water dimer by designing a triple-level strategy that combines (i) benchmark results close to the complete-basis limit of coupled-cluster theory with the single, double, triple, and quadruple excitations (CCSDTQ/CBS), (ii) a new hybrid meta density functional (M06CR) specifically optimized for reactions of Criegee intermediates, and (iii) variational transition-state theory with both variable rection coordinates and optimized reaction paths, with multidimensional tunneling, and with pressure effects. For (i) we have found that quadruple excitations are required to obtain quantitative reaction barriers, and we designed new composite methods and strategies to reach CCSDTQ/CBS accuracy. The present findings show that (i) the CH₂OO + HCHO reaction can make an important contribution to the sink of HCHO under wide atmospheric conditions in the gas phase and that (ii) CH₂OO + (H₂O)₂ dominates over the CH₂OO + H₂O reaction below 10 km.

Full insights into the roles of pH on hydroxylation of aromatic acids/bases and toxicity evaluation

Advanced oxidation processes (AOPs) based on hydroxyl radicals (•OH) are the most important technologies for the removal of bio-recalcitrant organic pollutants in industrial wastewater. The pH is one of the crucial environmental parameters that affect the removal efficiency of pollutants in AOPs. In this work, the mechanistic and kinetic insights into the roles of pH on the hydroxylation of five aromatic acids and bases in UV/H₂O₂ process have been investigated using theoretical calculation methods. Results show that the reactivity of •OH towards the twelve ionic/neutral species is positively correlated with electron-donating effect of substituents, which contributes to the positively pH-dependent reactivity of aromatic acids and bases towards •OH. The hydroxylation apparent rate constants (k_app, M^-1 s^-1) (at 298 K) increase as the pH values increase from about 1 to 10, but they decrease as the pH values increase from about 10 to 14. However, the best pH values for degradation are not around 10 because the [•OH] decreases continuously with the increasing pH values from 3 to 9.5. Combining the factors of k_app and [•OH], the best degradation pH values are around 5.5~7.5 for p-hydroxybenzoic acid, p-aminophenol, aniline and benzoic acid, 3.0~7.5 for phenol and 5.5~7.5 for mixed pollutants of these aromatic acids/bases in UV/H₂O₂ process. Moreover, a significant number of hydroxylation by-products are more toxic or harmful to aquatic organisms and rat (oral) than their parental pollutants. Altogether, this work provides comprehensive understanding of the roles of pH on •OH-initiated degradation behavior of aromatic acids and bases.

Theoretical Study of Ozonation of Methylparaben and Ethylparaben in Aqueous Solution

Parabens are widely employed in toothpaste, cosmetics, textiles, beverages, and preservatives, causing a serious environmental concern because they are endocrine-disrupting compounds (EDCs). As one of the highly reactive oxidants, ozone has a great effect on EDC removal. To understand the degradation and transformation of parabens in the aquatic environment and their toxicity to aquatic organisms, the degradation reaction of parabens initiated by O₃ was studied meticulously using quantum chemical calculations. The degradation process includes multiple initial reaction channels and consequent degradation pathways of the Criegee intermediates. Through thermodynamic data, the rate constants were computed using the transition state theory (TST). At a temperature of 298 K and a pressure of 1 atm, the calculated rate constants were 3.92 and 3.94 M^-1 s^-1 for methylparaben (MPB) and ethylparaben (EPB), respectively. The rate constants increased as the temperature increased or as the length of the alkyl chain on the benzene ring increased. Through the ecotoxicity assessment procedure, the ecotoxicity of parabens and the products in the degradation process can be assessed. Most degradation byproducts are either less toxic or nontoxic. Some byproducts are still harmful, such as oxalaldehyde (P2) and ethyl 2,3-dioxopropanoate (P10). Furthermore, the ecological toxicity of parabens increased with augmentation of the alkyl chain on the benzene ring. The effect of the alkyl chain length on the benzene ring in the compound cannot be ignored.