Non-ergodic fragmentation upon collision-induced activation of cysteine-water cluster cations.

Cysteine-water cluster cations Cys(H2O)3,6+ and Cys(H2O)3,6H+ are assembled in He droplets and probed by tandem mass spectrometry with collision-induced activation. Benchmark experimental data for this biologically important system are complemented with theory to elucidate the details of the collision-induced activation process. Experimental energy thresholds for successive release of water are compared to water dissociation energies from DFT calculations showing that clusters do not only fragment exclusively by sequential emission of single water molecules but also by the release of small water clusters. Release of clustered water is observed also in the ADMP (atom centered density matrix propagation) molecular dynamics model of small Cys(H2O)3+ and Cys(H2O)3H+ clusters. For large clusters Cys(H2O)6+ and Cys(H2O)6H+ the less computationally demanding statistical Microcanonical Metropolis Monte-Carlo method (M3C) is used to model the experimental fragmentation patterns. We are able to detail the energy redistribution in clusters upon collision activation. In the present case, about two thirds of the collision energy redistribute via an ergodic process, while the remaining one third is transferred into a non-ergodic channel leading to ejection of a single water molecule from the cluster. In contrast to molecular fragmentation, which can be well described by statistical models, modelling of collision-induced activation of weakly bound clusters requires inclusion of non-ergodic processes.

A general approach to study molecular fragmentation and energy redistribution after an ionizing event.

We propose to combine quantum chemical calculations, statistical mechanical methods, and photoionization and particle collision experiments to unravel the redistribution of internal energy of the furan cation and its dissociation pathways. This approach successfully reproduces the relative intensity of the different fragments as a function of the internal energy of the system in photoelectron-photoion coincidence experiments and the different mass spectra obtained when ions ranging from Ar+ to Xe25+ or electrons are used in collision experiments. It provides deep insights into the redistribution of the internal energy in the ionized molecule and its influence on the dissociation pathways and resulting charged fragments. The present pilot study demonstrates the efficiency of a statistical exchange of excitation energy among various degrees of freedom of the molecule and proves that the proposed approach is mature to be extended to more complex systems.

Unraveling the structural stability and the electronic structure of ThO2 clusters.

Unraveling the correlations between the geometry, the relative energy and the electronic structure of metal oxide nanostructures is crucial for a better control of their size, shape and properties. In this work, we investigated these correlations for stoichiometric thorium dioxide clusters ranging from ThO2 to Th8O16 using a chemically-driven geometry search algorithm in combination with state-of-the-art first principles calculations. This strategy allows us to homogeneously screen the potential energy surface of actinide oxide clusters for the first time. It is found that the presence of peroxo and superoxo groups tends to increase the total energy of the system by at least 3.5 eV and 7 eV, respectively. For the larger clusters, the presence of terminal oxygen atoms increases the energy by about 0.5 eV. Regarding the electronic structure, it is found that the HOMO-LUMO gap is larger in systems containing only bridging oxygen atoms (∼2-3.5 eV) than for systems containing oxo groups (∼1-3 eV), peroxo groups (∼0-2 eV), and superoxo groups (∼0-1 eV). Furthermore, while the LUMO is always dominated by thorium orbitals, the composition of the HOMO changes in the presence or the absence of oxo, peroxo and/or superoxo groups: in the presence of peroxo groups, it is dominated by thorium orbitals, in all other cases, it is dominated by oxygen orbitals, and is rather localized in the presence of terminal oxo or superoxo groups. These correlations are of great interest for synthesizing clusters with tailored properties, especially for applications in the field of nuclear energy and heterogeneous catalysis.

SCC-DFTB Parameters for Fe-C Interactions.

We present an optimized density-functional tight-binding (DFTB) parameterization for iron-based complexes based on the popular trans3d set of parameters. The transferability of the original and optimized parameterizations is assessed using a set of 50 iron complexes, which include carbonyl, cyanide, polypyridine, and cyclometalated ligands. DFTB-optimized structures predicted using the trans3d parameters show a good agreement with both experimental crystal geometries and density functional theory (DFT)-optimized structures for Fe-N bond lengths. Conversely, Fe-C bond lengths are systematically overestimated. We improve the accuracy of Fe-C interactions by truncating the Fe-O repulsive potential and reparameterizing the Fe-C repulsive potential using a training set of six isolated iron complexes. The new trans3d*-LANLFeC parameter set can produce accurate Fe-C bond lengths in both geometry optimizations and molecular dynamics (MD) simulations, without significantly affecting the accuracy of Fe-N bond lengths. Moreover, the potential energy curves of Fe-C interactions are considerably improved. This improved parameterization may open the door to accurate MD simulations at the DFTB level of theory for large systems containing iron complexes, such as sensitizer-semiconductor assemblies in dye-sensitized solar cells, that are not easily accessible with DFT approaches because of the large number of atoms.

Methyl Paraben and Carbamazepine in Water and Striped Catfish (Pseudoplatystoma magdaleniatum) in the Cauca and Magdalena Rivers.

Between 2017 and 2019, samplings were carried out in the San Jorge, Cauca and Magdalena River basins in Colombia, to determine the presence of methyl paraben and carbamazepine in water and Pseudoplatystoma magdaleniatum. For the analysis of the samples, a validation of the analytical method was performed, following the EPA method 1694 (Pharmaceutical and personal care products in water), with slight modifications. This was done by liquid-chromatography tandem mass spectrometry, for quantification of methyl paraben and carbamazepine, including parameters of linearity, accuracy precision and veracity. Carbamazepine was found in the Magdalena River at 8.03 ± 0.01 µg/L in transition season. In fish samples, methyl paraben and carbamazepine were detected in a range between 32 and 90.80 µg/kg in transition and dry seasons.

Fully versus constrained statistical fragmentation of carbon clusters and their heteronuclear derivatives.

The Microcanonical Metropolis Monte Carlo (MMMC) method has been shown to describe reasonably well fragmentation of clusters composed of identical atomic species. However, this is not so clear in the case of heteronuclear clusters as some regions of phase space might be inaccessible due to the different mobility of the different atomic species, the existence of large isomerization barriers, or the quite different chemical nature of the possible intermediate species. In this paper, we introduce a constrained statistical model that extends the range of applicability of the MMMC method to such mixed clusters. The method is applied to describe fragmentation of isolated clusters with high, moderate, and no heteronuclear character, namely, CnHm, CnN, and Cn clusters for which experimental fragmentation branching ratios are available in the literature. We show that the constrained statistical model describes fairly well fragmentation of CnHm clusters in contrast with the poor description provided by the fully statistical model. The latter model, however, works pretty well for both Cn and CnN clusters, thus showing that the ultimate reason for this discrepancy is the inability of the MMMC method to selectively explore the whole phase space. This conclusion has driven us to predict the fragmentation patterns of the C4N cluster for which experiments are not yet available.

Furan Fragmentation in the Gas Phase: New Insights from Statistical and Molecular Dynamics Calculations.

We present a complete exploration of the different fragmentation mechanisms of furan (C4H4O) operating at low and high energies. Three different theoretical approaches are combined to determine the structure of all possible reaction intermediates, many of them not described in previous studies, and a large number of pathways involving three types of fundamental elementary mechanisms: isomerization, fragmentation, and H/H2 loss processes (this last one was not yet explored). Our results are compared with the existing experimental and theoretical investigations for furan fragmentation. At low energies the first processes to appear are isomerization, which always implies the breaking of one C-O bond and one or several hydrogen transfers; at intermediate energies the fragmentation of the molecular skeleton becomes the most relevant mechanism; and H/H2 loss is the dominant processes at high energy. However, the three mechanisms are active in very wide energy ranges and, therefore, at most energies there is a competition among them.

Extracellular enzymatic activity of two hydrolases in wastewater treatment for biological nutrient removal.

Due to the complex nature of the wastewater (both domestic and non-domestic) composition, biological processes are widely used to remove nutrients, such as carbon (C), nitrogen (N), and phosphorous (P), which cause instability and hence contribute to the damage of water bodies. Systems with different configurations have been developed (including anaerobic, anoxic, and aerobic conditions) for the joint removal of carbon, nitrogen, and phosphorus. The goal of this research is to evaluate the extracellular activity of ß-glucosidase and phosphatase enzymes in a University of Cape Town (UCT) system fed with two synthetic wastewaters of different molecular complexity. Both types of waters have medium strength characteristics similar to those of domestic wastewater with a mean C/N/P ratio of 100:13:1. The operation parameters were hydraulic retention time (HRT) of 10 h, solid retention time (SRT) of 12 days, mean concentration of the influent in terms of chemical oxygen demand (COD), total Kjeldahl nitrogen (TKN), and total phosphorus (TP) of 600, 80, and 6 mg/L, respectively. According to the results obtained, statistically significant differences have been found in the extracellular enzyme activities with the evaluated wastewaters and in the units comprising the treatment system in some of the cases. An analysis of principal components showed that the extracellular enzymatic activity has been correlated to nutrient concentration in wastewater, biomass concentration in the system, and metabolic conditions of treatment phases. Additionally, this research has allowed determining an inverse relationship between wastewater biodegradability and the extracellular enzyme activity of ß-glucosidase and phosphatase. These results highlight the importance of including the analysis of biomass biochemical characteristics as control methods in wastewater treatment systems for the nutrient removal.

Structure, Ionization, and Fragmentation of Neutral and Positively Charged Hydrogenated Carbon Clusters: C(n)H(m)(q+) (n = 1-5, m = 1-4, q = 0-3).

In this work we present a systematic theoretical study of neutral and positively charged hydrogenated carbon clusters (C(n)H(m)(q+) with n = 1­5, m = 1­4, and q = 0­3). A large number of isomers and spin states (1490 in total) was investigated. For all of them, we optimized the geometry and computed the vibrational frequencies at the B3LYP/6-311++G(3df,2dp) level of theory; more accurate values of the electronic energy were obtained at the CCSD(T)/6-311++G(3df,2dp) level over the geometry previously obtained. From these simulations we evaluated several properties such as relative energies between isomers, adiabatic and vertical ionization potentials, and dissociation energies of several fragmentation channels. A new analysis technique is proposed to evaluate a large number of fragmentation channels in a wide energy range.

Communication: unraveling the (4)He droplet-mediated soft-landing from ab initio-assisted and time-resolved density functional simulations: Au@(4)He300/TiO2(110).

An ab-initio-based methodological scheme for He-surface interactions and zero-temperature time-dependent density functional theory for superfluid (4)He droplets motion are combined to follow the short-time collision dynamics of the Au@(4)He300 system with the TiO2(110) surface. This composite approach demonstrates the (4)He droplet-assisted sticking of the metal species to the surface at low landing energy (below 0.15 eV/atom), thus providing the first theoretical evidence of the experimentally observed (4)He droplet-mediated soft-landing deposition of metal nanoparticles on solid surfaces [Mozhayskiy et al., J. Chem. Phys. 127, 094701 (2007) and Loginov et al., J. Phys. Chem. A 115, 7199 (2011)].

Efficient Parameterization of Density Functional Tight-Binding for 5f-Elements: A Th-O Case Study.

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th-O system, which was subsequently applied to the study of ThO2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO2. Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

Solvent states and spectroscopy of doped helium clusters as a quantum-chemistry-like problem.

The Full-Configuration-Interaction Nuclear-Orbital (FCI-NO) approach [J. Chem. Phys., 2009, 131, 19401], as the implementation of the quantum-chemistry ansatz, is overviewed and applied to (He)N-Cl2(X) clusters (N≤ 4). The ground and excited states of both fermionic (3)He and bosonic (4)He [see also, J. Phys. Chem. Lett., 2012, 2, 2145] clusters are studied. It is shown that the FCI-NO approach allows us to overcome three main difficulties: (1) the Fermi-Dirac (Bose-Einstein) nuclear statistics; (2) the wide (highly anharmonic) amplitudes of the He-dopant and He-He motions; and (3) both the weakly attractive (long-range) and the strongly repulsive (short-range) interaction between the helium atoms. Special emphasis is placed on the dependence of the cluster properties on the number of helium atoms, and on the comparison between the two helium isotopes. In particular, we analyze the analogies between quantum rings comprising electrons and (3)He atoms. The synthetic vibro-rotational Raman spectra of Cl2(X) immersed in ((3,4)He)N clusters (N≤ 4) are discussed as a function of the cluster size and the nuclear statistics. It is shown that the Coriolis couplings play a key role in modifying the spectral dopant profile in (3)He. Finally, we point out possible directions for future research using the quantum-chemistry ansatz.

Including nuclear quantum effects into highly correlated electronic structure calculations of weakly bound systems.

An interface between the APMO code and the electronic structure package MOLPRO is presented. The any particle molecular orbital APMO code [González et al., Int. J. Quantum Chem. 108, 1742 (2008)] implements the model where electrons and light nuclei are treated simultaneously at Hartree-Fock or second-order Möller-Plesset levels of theory. The APMO-MOLPRO interface allows to include high-level electronic correlation as implemented in the MOLPRO package and to describe nuclear quantum effects at Hartree-Fock level of theory with the APMO code. Different model systems illustrate the implementation: (4)He2 dimer as a protype of a weakly bound van der Waals system; isotopomers of [He-H-He](+) molecule as an example of a hydrogen bonded system; and molecular hydrogen to compare with very accurate non-Born-Oppenheimer calculations. The possible improvements and future developments are outlined.

Automated software for counting and measuring Hyalella genus using artificial intelligence.

Amphipods belonging to the Hyalella genus are macroinvertebrates that inhabit aquatic environments. They are of particular interest in areas such as limnology and ecotoxicology, where data on the number of Hyalella individuals and their allometric measurements are used to assess the environmental dynamics of aquatic ecosystems. In this study, we introduce HyACS, a software tool that uses a model developed with the YOLOv3's architecture to detect individuals, and digital image processing techniques to extract morphological metrics of the Hyalella genus. The software detects body metrics of length, arc length, maximum width, eccentricity, perimeter, and area of Hyalella individuals, using basic imaging capture equipment. The performance metrics indicate that the model developed can achieve high prediction levels, with an accuracy above 90% for the correct identification of individuals. It can perform up to four times faster than traditional visual counting methods and provide precise morphological measurements of Hyalella individuals, which may improve further studies of the species populations and enhance their use as bioindicators of water quality.

Helium mediated deposition: modeling the He-TiO2(110)-(1×1) interaction potential and application to the collision of a helium droplet from density functional calculations.

This paper is the first of a two-part series dealing with quantum-mechanical (density-functional-based) studies of helium-mediated deposition of catalytic species on the rutile TiO(2)(110)-(1×1) surface. The interaction of helium with the TiO(2)(110)-(1×1) surface is first evaluated using the Perdew-Burke-Ernzerhof functional at a numerical grid dense enough to build an analytical three-dimensional potential energy surface. Three (two prototype) potential models for the He-surface interaction in helium scattering calculations are analyzed to build the analytical potential energy surface: (1) the hard-corrugated-wall potential model; (2) the corrugated-Morse potential model; and (3) the three-dimensional Morse potential model. Different model potentials are then used to study the dynamics upon collision of a (4)He(300) cluster with the TiO(2)(110) surface at zero temperature within the framework of a time-dependent density-functional approach for the quantum fluid [D. Mateo, D. Jin, M. Barranco, and M. Pi, J. Chem. Phys. 134, 044507 (2011)] and classical dynamics calculations. The laterally averaged density functional theory-based potential with an added long-range dispersion interaction term is further applied. At variance with classical dynamics calculations, showing helium droplet splashing out of the surface at impact, the time evolution of the macroscopic helium wave-function predicts that the helium droplet spreads on the rutile surface and leads to the formation of a thin film above the substrate. This work thus provides a basis for simulating helium mediated deposition of metallic clusters embedded within helium nanodroplets.

[The influence of macrophytes on rotifer and microcrustacean assemblage in a tropical floodplain]. / Influencia de las macrófitas sobre la estructura poblacional de rotíferos y microscrustáceos en un plano de inundación tropical.

Most studies on zooplankton ecology have been conducted in open waters. However, it has been considered of great importance to extend such studies to other habitats, such as those generated of macrophytes. We studied the spatial and temporal variation of the microcrustacean and rotifer assemblage structures associated with macrophytes, and compare them with the variation exhibited in open waters. Integrated samples were collected for zooplankton and phytoplankton biomass using a Schindler bottle, in four open water sites and four other sites covered by macrophytes in the floodplain complex of Ayapel (Córdoba, Colombia) during different limnimetric levels. The significant differences in the structure were evaluated using Kruskal & Wallis and discriminant analyses, and the similarity among sampling sites was evaluated using Bray & Curtis analysis. Zooplanktonic richness was favored by macrophytes. However, we did not find a constant spatial pattern for density, and only particular trends apparently conditionated by flood pulses. The presence of Eichhornia azurea (Pontederiaceae) contributed in a significant way with an increase in the diversity and density of benthic taxa. The density of some zooplanktonic groups was related with environmental conditions and phytoplanktonic biomass.

Quantum solvent states and rovibrational spectra of small doped (3)He clusters through the full-configuration-interaction nuclear orbital approach: The ((3)He)(N)-Cl(2)(X) case (N

A full-configuration-interaction nuclear orbital treatment has been recently developed as a benchmark quantum-chemistry-like method to study small doped (3)He clusters [M. P. de Lara-Castells et al., J. Chem. Phys. 125, 221101 (2006)]. Our objective in this paper is to extend our previous study on ((3)He)(N)-Cl(2)(B) clusters, using an enhanced implementation that allows employing very large one-particle basis sets [M. P. de Lara-Castells et al., J. Chem. Phys. 131, 194101 (2009)], and apply the method to the ((3)He)(N)-Cl(2)(X) case, using both a semiempirical T-shaped and an ab initio He-dopant potential with minima at both T-shaped and linear conformations. Calculations of the ground and low-lying excited solvent states stress the key role played by the anisotropy of the He-dopant interaction in determining the global energies and the structuring of the (3)He atoms around the dopant. Whereas (3)He atoms are localized in a broad belt around the molecular axis in ground-state N-sized complexes with N=1-3, irrespective of using the T-shaped or the ab initio He-dopant potential function, the dopant species becomes fully coated by just four (3)He atoms when the He-dopant potential also has a minimum at linear configurations. However, excited solvent states with a central ring-type clustering of the host molecule are found to be very close in energy with the ground state by using the ab initio potential function. A microscopic analysis of this behavior is provided. Additional simulations of the molecular rovibrational Raman spectra, also including excited solvent states, provide further insights into the importance of proper modeling the anisotropy of the He-dopant interaction in these weakly bound systems and of taking into account the low-lying excitations.

Development of Density Functional Tight-Binding Parameters Using Relative Energy Fitting and Particle Swarm Optimization.

We provide a strategy to optimize density functional tight-binding (DFTB) parameterization for the calculation of the structures and properties of organic molecules consisting of hydrogen, carbon, nitrogen, and oxygen. We utilize an objective function based on similarity measurements and the Particle Swarm Optimization (PSO) method to find an optimal set of parameters. This objective function considers not only the common DFTB descriptors of binding energies and atomic forces but also incorporates relative energies of isomers into the fitting procedure for more chemistry-driven results. The quality in the description of the binding energies and atomic forces is measured based on the Ballester similarity index and relative energies through a similarity index induced by the Levenshtein edit distance to quantify the correct energetic order of isomers. Training and testing datasets were created to include all relevant chemical functional groups. The accuracy of this strategy is assessed, and its range of applicability is discussed by comparison against our previous parameterization [A. Krishnapriyan, et al., J. Chem. Theory Comput. 13, 6191 (2017)]. The improved performance of the new DFTB parameterization is validated with respect to the density functional theory large datasets QM-9 [R. Ramakrishnan, et al., Sci. Data 1, 140022 (2014)] and ANI-1 [J. S. Smith, et al., Sci. Data 4, 170193 (2017)], where excellent agreement is found between the structures and properties available in these datasets, and the ones obtained with DFTB.

Dissociation dynamics of the diamondoid adamantane upon photoionization by XUV femtosecond pulses.

This work presents a photodissociation study of the diamondoid adamantane using extreme ultraviolet femtosecond pulses. The fragmentation dynamics of the dication is unraveled by the use of advanced ion and electron spectroscopy giving access to the dissociation channels as well as their energetics. To get insight into the fragmentation dynamics, we use a theoretical approach combining potential energy surface determination, statistical fragmentation methods and molecular dynamics simulations. We demonstrate that the dissociation dynamics of adamantane dications takes place in a two-step process: barrierless cage opening followed by Coulomb repulsion-driven fragmentation.

M3C: A Computational Approach To Describe Statistical Fragmentation of Excited Molecules and Clusters.

The Microcanonical Metropolis Monte Carlo method, based on a random sampling of the density of states, is revisited for the study of molecular fragmentation in the gas phase (isolated molecules, atomic and molecular clusters, complex biomolecules, etc.). A random walk or uniform random sampling in the configurational space (atomic positions) and a uniform random sampling of the relative orientation, vibrational energy, and chemical composition of the fragments is used to estimate the density of states of the system, which is continuously updated as the random sampling populates individual states. The validity and usefulness of the method is demonstrated by applying it to evaluate the caloric curve of a weakly bound rare gas cluster (Ar13), to interpret the fragmentation of highly excited small neutral and singly positively charged carbon clusters (Cn, n = 5,7,9 and Cn+, n = 4,5) and to simulate the mass spectrum of the acetylene molecule (C2H2).