Characterizing moisture uptake and plasticization effects of water on amorphous amylose starch models using molecular dynamics methods.

Dynamics and thermophysical properties of amorphous starch were explored using molecular dynamics (MD) simulations. Using the OPLS3e force field, simulations of short amylose chains in water were performed to determine force field accuracy. Using well-tempered metadynamics, a free energy map of the two glycosidic angles of an amylose molecule was constructed and compared with other modern force fields. Good agreement of torsional sampling for both solvated and amorphous amylose starch models was observed. Using combined grand canonical Monte Carlo (GCMC)/MD simulations, a moisture sorption isotherm curve is predicted along with temperature dependence. Concentration-dependent activation energies for water transport agree quantitatively with previous experiments. Finally, the plasticization effect of moisture content on amorphous starch was investigated. Predicted glass transition temperature (Tg) depression as a function of moisture content is in line with experimental trends. Further, our calculations provide a value for the dry Tg for amorphous starch, a value which no experimental value is available.

Molecular Design Based on Donor-Weak Donor Scaffold for Blue Thermally-Activated Delayed Fluorescence Designed by Combinatorial DFT Calculations.

Quantum chemical calculations are necessary to develop advanced emitter materials showing thermally-activated delayed fluorescence (TADF) for organic light-emitting diodes (OLEDs). However, calculation costs become problematic when more accurate functionals were used, therefore it is judicious to use a multimethod approach for efficiency. Here we employed combinatorial chemistry in silico to develop the deep blue TADF materials with a new concept of homo-junction design. The homo-junction materials containing TADF candidates designed by calculation were synthesized and analyzed. We found that these materials showed the emission from charge transfer (CT) state, and the clear delayed emission was provided in solid state. Because the homo-junction TADF materials showed three exponential decayed emission in solid state, we employed novel four-state kinetic analysis.

Massive Theoretical Screen of Hole Conducting Organic Materials in the Heteroacene Family by Using a Cloud-Computing Environment.

Materials exhibiting higher mobilities than conventional organic semiconducting materials such as fullerenes and fused thiophenes are in high demand for applications in printed electronics. To discover new molecules in the heteroacene family that might show improved charge mobility, a massive theoretical screen of hole conducting properties of molecules was performed by using a cloud-computing environment. Over 7 000 000 structures of fused furans, thiophenes and selenophenes were generated and 250 000 structures were randomly selected to perform density functional theory (DFT) calculations of hole reorganization energies. The lowest hole reorganization energy calculated was 0.0548 eV for a fused thioacene having 8 aromatics rings. Hole mobilities of compounds with the lowest 130 reorganization energy were further processed by applying combined DFT and molecular dynamics (MD) methods. The highest mobility calculated was 1.02 and 9.65 cm2/(V s) based on percolation and disorder theory, respectively, for compounds containing selenium atoms with 8 aromatic rings. These values are about 20 times higher than those for dinaphthothienothiophene (DNTT).

Structural convergence properties of amorphous InGaZnO4 from simulated liquid-quench methods.

The study of structural properties of amorphous structures is complicated by the lack of long-range order and necessitates the use of both cutting-edge computer modeling and experimental techniques. With regards to the computer modeling, many questions on convergence arise when trying to assess the accuracy of a simulated system. What cell size maximizes the accuracy while remaining computationally efficient? More importantly, does averaging multiple smaller cells adequately describe features found in bulk amorphous materials? How small is too small? The aims of this work are: (1) to report a newly developed set of pair potentials for InGaZnO4 and (2) to explore the effects of structural parameters such as simulation cell size and numbers on the structural convergence of amorphous InGaZnO4. The total number of formula units considered over all runs is found to be the critical factor in convergence as long as the cell considered contains a minimum of circa fifteen formula units. There is qualitative agreement between these simulations and X-ray total scattering data - peak trends and locations are consistently reproduced while intensities are weaker. These new IGZO pair potentials are a valuable starting point for future structural refinement efforts.

Group additivity-Pourbaix diagrams advocate thermodynamically stable nanoscale clusters in aqueous environments.

The characterization of water-based corrosion, geochemical, environmental and catalytic processes rely on the accurate depiction of stable phases in a water environment. The process is aided by Pourbaix diagrams, which map the equilibrium solid and solution phases under varying conditions of pH and electrochemical potential. Recently, metastable or possibly stable nanometric aqueous clusters have been proposed as intermediate species in non-classical nucleation processes. Herein, we describe a Group Additivity approach to obtain Pourbaix diagrams with full consideration of multimeric cluster speciation from computations. Comparisons with existing titration results from experiments yield excellent agreement. Applying this Group Additivity-Pourbaix approach to Group 13 elements, we arrive at a quantitative evaluation of cluster stability, as a function of pH and concentration, and present compelling support for not only metastable but also thermodynamically stable multimeric clusters in aqueous solutions.

Catalytic Efficiency Is a Function of How Rhodium(I) (5 + 2) Catalysts Accommodate a Conserved Substrate Transition State Geometry: Induced Fit Model for Explaining Transition Metal Catalysis.

The origins of differential catalytic reactivities of four Rh(I) catalysts and their derivatives in the (5 + 2) cycloaddition reaction were elucidated using density functional theory. Computed free energy spans are in excellent agreement with known experimental rates. For every catalyst, the substrate geometries in the transition state remained constant (<0.1 Å RMSD for atoms involved in bond-making and -breaking processes). Catalytic efficiency is shown to be a function of how well the catalyst accommodates the substrate transition state geometry and electronics. This shows that the induced fit model for explaining biological catalysis may be relevant to transition metal catalysis. This could serve as a general model for understanding the origins of efficiencies of catalytic reactions.

Mechanism and the origins of stereospecificity in copper-catalyzed ring expansion of vinyl oxiranes: a traceless dual transition-metal-mediated process.

Density functional theory computations of the Cu-catalyzed ring expansion of vinyloxiranes is mediated by a traceless dual Cu(I)-catalyst mechanism. Overall, the reaction involves a monomeric Cu(I)-catalyst, but a single key step, the Cu migration, requires two Cu(I)-catalysts for the transformation. This dual-Cu step is found to be a true double Cu(I) transition state rather than a single Cu(I) transition state in the presence of an adventitious, spectator Cu(I). Both Cu(I) catalysts are involved in the bond forming and breaking process. The single Cu(I) transition state is not a stationary point on the potential energy surface. Interestingly, the reductive elimination is rate-determining for the major diastereomeric product, while the Cu(I) migration step is rate-determining for the minor. Thus, while the reaction requires dual Cu(I) activation to proceed, kinetically, the presence of the dual-Cu(I) step is untraceable. The diastereospecificity of this reaction is controlled by the Cu migration step. Suprafacial migration is favored over antarafacial migration due to the distorted Cu π-allyl in the latter.