Autonomous intelligent agents for accelerated materials discovery

We present an end-to-end computational system for autonomous materials discovery. The system aims for cost-effective optimization in large, high-dimensional search spaces of materials by adopting a sequential, agent-based approach to deciding which experiments to carry out. In choosing next experiments, agents can make use of past knowledge, surrogate models, logic, thermodynamic or other physical constructs, heuristic rules, and different exploration–exploitation strategies. We show a series of examples for (i) how the discovery campaigns for finding materials satisfying a relative stability objective can be simulated to design new agents, and (ii) how those agents can be deployed in real discovery campaigns to control experiments run externally, such as the cloud-based density functional theory simulations in this work. In a sample set of 16 campaigns covering a range of binary and ternary chemistries including metal oxides, phosphides, sulfides and alloys, this autonomous platform found 383 new stable or nearly stable materials with no intervention by the researchers.

We present an end-to-end computational system for autonomous materials discovery.     

Kinetic investigation of the reactions of S-4-nitrophenyl 4-substituted thiobenzoates with secondary alicyclic amines in aqueous ethanol

The reactions of S-4-nitrophenyl 4-X-substituted thiobenzoates (X = H, Cl, and NO(2): 1, 2, and 3, respectively) with a series of secondary alicyclic amines (SAA) were subjected to a kinetic investigation in 44 wt % ethanol-water, at 25.0 degrees C and an ionic strength of 0.2 M (KCl). The reactions were followed spectrophotometrically by monitoring the release of 4-nitrobenzenethiolate anion at 420-425 nm. Under excess amine, pseudo-first-order rate constants (k(obsd)) are obtained for all reactions. The plots of k(obsd) vs [SAA] at constant pH are linear with the slope (k(N)) independent of pH. The statistically corrected Br?nsted-type plots (log k(N)/q vs pK(a) + log p/q) for the reactions of 1 and 2 are nonlinear with slopes at high pK(a), beta(1) = 0.27 and 0.10, respectively, and slopes at low pK(a), beta(2) = 0.86 and 0.84, respectively. The Br?nsted curvature is centered at pK(a) (pK(a)(0)) 10.0 and 10.4, respectively. The reactions of SAA with 3 exhibit a linear Br?nsted-type plot of slope 0.81. These results are consistent with a stepwise mechanism, through a zwitterionic tetrahedral intermediate (T(+/-)). For the reactions of 1 and 2, there is a change in rate-determining step with amine basicity, from T(+/-) breakdown to products at low pK(a), to T(+/-) formation at high pK(a). For the reactions of 3, breakdown to products of T(+/-) is rate limiting for all the SAA series (pK(a)(0) > 11). The increasing pK(a)(0) value as the substituent in the acyl group becomes more electron withdrawing is attributed to an increasing nucleofugality of SAA from T(+/-). The greater pK(a)(0) value for the reactions of SAA with 1, relative to that found in the pyridinolysis of 2,4-dinitrophenyl benzoate (pK(a)(0) = 9.5), is explained by the greater nucleofugality from T(+/-) of the former amines, compared to isobasic pyridines, and the greater leaving ability from T(+/-) of 2,4-dinitrophenoxide relative to 4-nitrobenzenethiolate.     

Kinetics and mechanism of the reactions of S-2,4-dinitrophenyl 4-substituted thiobenzoates with secondary alicyclic amines

[reaction: see text] The title reactions, in 44 wt % ethanol-water at 25.0 degrees C, exhibit slightly curved Br?nsted-type plots (log kN versus pKa of amines) with slopes beta1 = 0.1-0.44 (at high pKa) and beta2 ca. 0.7 (at low pKa). The magnitude of some of these slopes, together with the fact that the curvature center (pKa(0) = 9.5-10.8) does not change with the electronic effects of the benzoyl substituent, suggests that these reactions are not stepwise, but concerted.     

An efficient coarse-grained approach for the electron transport through large molecular systems under dephasing environment

Dephasing effects in electron transport in molecular systems connected between contacts average out the quantum characteristics of the system, forming a bridge to the classical behavior as the size of the system increases. For the evaluation of the conductance of the molecular systems which have sizes within this boundary domain, it is necessary to include these dephasing effects. These effects can be calculated by using the D’Amato-Pastawski model. However, this method is computationally demanding for large molecular systems since transmission functions for all pairs of atomic orbitals need to be calculated. To overcome this difficulty, we develop an efficient coarse-grained model for the calculation of conductance of molecular junctions including decoherence. By analyzing the relationship between chemical potential and inter-molecular coupling, we find that the chemical potential drops stepwise in the systems with weaker inter-unit coupling. Using this property, an efficient coarse-grained algorithm which can reduce computational costs considerably without losing the accuracy is derived and applied to one-dimensional organic systems as a demonstration. This model can be used for the study of the orientation dependence of conductivity in various phases (amorphous, crystals, and polymers) of large molecular systems such as organic semiconducting materials.     

Dual Nucleophilic Substitution Reactions of O,O‐Diethyl 2,4‐dinitrophenyl Phosphate and Thionophosphate Triesters

The reactions of the title compounds with phenoxides, secondary alicyclic (SA) amines, and pyridines, in 44 wt% ethanol–water, at 25°C and an ionic strength of 0.2 M, were subjected to kinetic and product studies. From analytical techniques (HPLC and NMR), two pathways were detected (nucleophilic attack at the phosphoryl center and at the C‐1 aromatic carbon) for the reactions of all the nucleophiles with the phosphate ( 2 ) and for the pyridinolysis of the thionophosphate ( 1 ). Only aromatic nucleophilic substitution was found for the reactions of 1 with phenoxides and SA amines. For the dual reactions, the nucleophilic rate constants (k_N) were separated in two terms: $k_{\rm N}^{\rm P}$ and $k_{\rm N}^{{\rm Ar}}$, which are the rate constants for the corresponding electrophilic centers. The absence of a break in the Brønsted‐type plots for the attack at P is consistent with concerted mechanisms. The Brønsted slopes, β_Ar 0.32–0.71, for the attack at the aromatic C‐1, are in agreement with stepwise mechanisms where formation of a Meisenheimer complex is the rate‐determining step. © 2013 Wiley Periodicals, Inc. Int J Chem Kinet 45: 202–211, 2013     

Catalytic asymmetric synthesis of β-triazolyl amino alcohols by asymmetric transfer hydrogenation of α-triazolyl amino alkanones

The synthesis of optically active β-triazolyl amino alcohols was carried out via ruthenium catalyzed asymmetric transfer hydrogenation of α-triazolyl amino alkanones. This reaction proceeds under mild reaction conditions with up to 99% yield and 99.9% enantiomeric excess (ee). This protocol was applied to the synthesis of an enantiopure antitubercular agent and its arylated product with retention in enantiomeric purity. The absolute configuration at the stereogenic center of the chiral product as found to be (S).     

Electron spin resonance in pyrolyzed lignin

An investigation has been made into electron spin resonance (ESR) occurring in the ranges of carbons prepared from periodate lignin by vacuum pyrolysis. Line shape, line width, g-factors, saturation effects, and free spin concentration have been measured and the effect of oxygen on these parameters has been examined. A g-shift and lines with inhomogeneous broadening have been observed. In the presence of oxygen the lower-temperature carbons showed an unexpected increase in the detectable free spin concentration. The discussion includes factors affecting the nature of the free spins in the pyrolyzed material.     

Novel Capsular Aggregates from Flexible Tripodal Triureas with Cs Symmetry

Tris(2‐ and 3‐ureidobenzyl)amines with C_s symmetry self‐assemble in solution forming mixtures of regioisomeric capsular aggregates, one of which is chiral and the other centrosymmetric. Under certain conditions, a predominance of the centrosymmetric regioisomer is found before equilibrium, that is, a mixture close to the statistical ratio of the two species is reached. In the solid state, there is a preference for the centrosymmetric capsules. Molecular models of both regioisomeric aggregates have been built and analyzed for comparison. Guests inside capsules formed by self‐assembly of desymmetrized tris(3‐ureidobenzyl)amines feel different magnetic environments, depending on whether they are inside a chiral or an achiral regioisomeric container. Of special significance are the experiments with a more flexible triurea endowed with an ureidopropylic arm, which self‐assembles with the same efficiency as the more rigid tris(ureidobenzyl)amines.