On the influence of nuclear fluctuations on calculated NMR shieldings of benzene and ethylene: a Feynman path integral–ab initio investigation*

The absolute magnetic shieldings of benzene and ethylene have been theoretically studied under the conditions of thermal equilibrium, i.e., under explicit consideration of the nuclear degrees of freedom. For this purpose we have combined the Feynman path integral quantum Monte Carlo (PIMC) formalism with the gauge‐including atomic orbital (GIAO) approach in the Hartree–Fock (HF) approximation. The HF operator has been employed to derive the NMR parameters of the two hydrocarbons via an ensemble averaging over large sets of molecular configurations that are populated in thermal equilibrium. The nuclear fluctuations are responsible for a deshielding of the nuclei relative to the shieldings at the vibrationless minimum of the potential energy surface (PES). The influence of the nuclear degrees of freedom is largest for the isotropic part of the ¹³C shielding tensor. The theoretical results can be explained on the basis of simple geometrical considerations. The bond lengths in thermal equilibrium are larger than the bond lengths at the minimum of the PES. This length enhancement is the prerequisite for a deshielding of the nuclei in thermal equilibrium. The vibrational corrections of the nuclear magnetic resonance (NMR) parameters of benzene and ethylene are quantum driven; classical thermal degrees of freedom of the nuclei are of minor importance. Conceptual problems of theoretical studies of NMR parameters on the basis of a single molecular geometry are emphasized. The influence of the spatial uncertainty of the nuclei becomes decisive in molecules with light atoms. It is pointed out that the combination of the PIMC formalism with electronic Hamiltonians of state‐of‐the‐art quality renders possible accurate determinations of NMR parameters. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem 86: 280–296, 2002     

Viscoelastic response and interlaminar delamination resistance of epoxy/glass fiber/functionalized graphene oxide multi-scale composites

Graphene oxide (GO) was functionalized using three different diamines, namely ethylenediamine (EDA), 4,4′-diaminodiphenyl sulfone (DDS) and p-phenylenediamine (PPD) to reinforce an epoxy/glass fiber (EP/GF) composite laminate, with the aim of improving the overall composite mechanical performance. Different mechanical characterization techniques were used to determine the mechanical performance, including: tensile stress strain, double cantilever beam (DCB) mode-I fracture toughness and dynamic mechanical thermal analysis (DMTA). Scanning electron microscopy (SEM) was used to support the results and conclusions. The results demonstrated remarkable enhancements in the mechanical performance of EP/GF composite laminates by incorporation of functionalized graphene oxide (FGO) nanofiller, whilst the mechanical performance of the GO reinforced composite only improved marginally. Finally, the mechanical performance of the EP/GF/FGO multi-scale composites was found to be dependent on the type of FGO functional groups; of which EDA exhibited the highest performance. These investigations confirmed that the EDA-FGO-reinforced EP/GF composites possess excellent potential to be used as multifunctional engineering materials in industrial applications.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

A new 4-D non-equilibrium fractional-order chaotic system and its circuit implementation

A new 4-D fractional-order chaotic system without equilibrium point is proposed in this paper. There is no chaotic behavior for its corresponding integer-order system. By computer simulations, we find complex dynamical behaviors in this system, and obtain that the lowest order for exhibiting a chaotic attractor is 3.2. We also design an electronic circuit to realize this 4-D fractional-order chaotic system and present some experiment results.     

Modulated nematic structures and chiral symmetry breaking in 2D

ABSTRACT

We have studied the properties of biaxial particles interacting via an anisotropic pair potential, involving second-rank quadrupolar and third-rank octupolar coupling terms, using Monte Carlo simulation. The particles occupy the sites of a 2D square lattice and the interactions are restricted to nearest neighbours. The system exhibits spontaneous chiral symmetry breaking from an isotropic phase to a chiral modulated nematic phase, composed of ambidextrous chiral domains. When twofold axes of quadrupolar and octupolar tensors coincide this modulated phase appears to be the ambidextrous cholesteric phase with pitch comparable to a few lattice spacings. The associated phase transition is first order.     

Simulation of mesogenic diruthenium tetracarboxylates: Development of a force field for coordination polymers of the MMX type

A classical molecular mechanics force field, able to simulate coordination polymers (CP) based on ruthenium carboxylates (Ru₂(O₂CR_eq)₄L_ax) (eq = equatorial group containing aliphatic chains, L_ax= axial ligand), has been developed. New parameters extracted from experimental data and quantum calculations on short aliphatic chains model systems were included in the generalized AMBER force field. The proposed parametrization was evaluated using model systems with known structure, containing either short or long aliphatic chains; experimental results were reproduced satisfactorily. This modified force field, although in a preliminary stage, could then be applied to long chain liquid crystalline compounds. The resulting atomistic simulations allowed assessing the relative influence of the factors determining the CP conformation, determinant for the physical properties of these materials. © 2013 Wiley Periodicals, Inc.     

Binary hairy nanoparticles: Recent progress in theory and simulations

Binary polymer brushes, including mixed homopolymer brushes and diblock copolymer brushes, are an attractive class of environmentally responsive nanostructured materials. Owing to microphase separation of the two chemically distinct components in the brush, multifaceted nanomaterials with functionalized and patterned surfaces can be obtained. This review summarizes recent progress on the theory and simulations related to binary polymer brushes grafted to flat, spherical, and cylindrical substrates, with a focus on patterned morphologies of multifaceted hairy nanoparticles, an intriguing class of hybrid nanostructured particles (e.g., nanospheres and nanorods). In particular, powerful field theory and particle-based simulations suitable for revealing novel structures on these patterned surfaces, including self-consistent field theory and dissipative particle dynamics simulations, are emphasized. The unsolved yet critical issues in this research field, such as dynamic response of binary polymer brushes to environmental stimuli and the hierarchical self-assembly of binary hairy nanoparticles, are briefly discussed. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2014 , 52, 1583–1599     

Dynamical Bifurcation in Gas‐Phase XH− + CH3Y SN2 Reactions: The Role of Energy Flow and Redistribution in Avoiding the Minimum Energy Path

The gas‐phase reactions of XH^? (X=O, S) + CH₃Y (Y=F, Cl, Br) span nearly the whole range of S_N2 pathways, and show an intrinsic reaction coordinate (IRC) (minimum energy path) with a deep well owing to the CH₃XH???Y^? (or CH₃S^????HF) hydrogen‐bonded postreaction complex. MP2 quasiclassical‐type direct dynamics starting at the [HX???CH₃???Y]^? transition‐state (TS) structure reveal distinct mechanistic behaviors. Trajectories that yield the separated CH₃XH+Y^? (or CH₃S^?+HF) products directly are non‐IRC, whereas those that sample the CH₃XH???Y^? (or CH₃S^????HF) complex are IRC. The IRCIRC/non‐IRC ratios of 90:10, 40:60, 25:75, 2:98, 0:100, and 0:100 are obtained for (X, Y)=(S, F), (O, F), (S, Cl), (S, Br), (O, Cl), and (O, Br), respectively. The properties of the energy profiles after the TS cannot provide a rationalization of these results. Analysis of the energy flow in dynamics shows that the trajectories cross a dynamical bifurcation, and that the inability to follow the minimum energy path arises from long vibration periods of the X?C???Y bending mode. The partition of the available energy to the products into vibrational, rotational, and translational energies reveals that if the vibrational contribution is more than 80 %, non‐IRC behavior dominates, unless the relative fraction of the rotational and translational components is similar, in which case a richer dynamical mechanism is shown, with an IRC/non‐IRC ratio that correlates to this relative fraction.     

Multipod structures of lamellae‐forming diblock copolymers in three‐dimensional confinement spaces: Experimental observation and computer simulation

The three‐dimensional (3D) confinement effect on the microphase‐separated structure of a diblock copolymer was investigated both experimentally and computationally. Block copolymer nanoparticles were prepared by adding a poor solvent into a block copolymer solution and subsequently evaporating the good solvent. The 3D structures of the nanoparticles were quantitatively determined with transmission electron microtomography (TEMT). TEMT observations revealed that various complex structures, including tennis‐ball, mushroom‐like, and multipod structures, were formed in the 3D confinement. Detailed structural analysis, showed that one block of the diblock copolymer slightly prefers to segregate into the particle surface compared with the other block. The observed structures were further elaborated using cell dynamics computer simulation. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 1702–1709     

Impact of proliferation strategies on food web viability in a model with closed nutrient cycle

A food web model with a closed nutrient cycle is presented and analyzed via Monte Carlo simulations. The model consists of three trophic levels, each of which is populated by animals of one distinct species. While the species at the intermediate level feeds on the basal species, and is eaten by the predators living at the highest level, the basal species itself uses the detritus of animals from higher levels as the food resource. The individual organisms remain localized, but the species can invade new lattice areas via proliferation. The impact of different proliferation strategies on the viability of the system is investigated. From the phase diagrams generated in the simulations it follows that in general a strategy with the intermediate level species searching for food is the best for the survival of the system. The results indicate that both the intermediate and top level species play a critical role in maintaining the structure of the system.