Viscoelasticity of Low-Molecular-Weight Polyelectrolytes.

Shear-thickening fluids that absorb the impact energy of high-velocity projectiles are of great interest for aerospace and body-armor applications. In such a frame, we investigate transient states of neat and aqueous polyelectrolytes (PE) having low molecular weights and containing poly([2-(methacryloyloxy)ethyl]trimethylammonium) as polycations and poly(acrylamide-co-acrylic acid) as polyanions. We compare results with those of bulk water. We employ nonequilibrium molecular dynamics to simulate oscillatory shear, mainly in the linear viscoelastic regime. We find that neat PE exhibits properties of a viscoelastic solid, whereas water and the aqueous mixture of PE conform to viscoelastic liquids with Maxwellian behavior at low angular frequencies. Terminal relaxation times are ∼0.499 and ∼1.385 ps for water and the aqueous mixture of PE, respectively. At high angular frequencies, storage moduli show anomalous behaviors that correspond to transitions between shear thinning and shear thickening in complex shear viscosities. The change in potential energy with the increase of the angular frequency is mainly driven by intramolecular interactions for neat PE, whereas short-range Coulomb interactions are the major contributions for water and the aqueous mixture of PE. Upon observation of the molecular configurations, only the local polyionic structure in the aqueous mixture of PE shows improvement when increasing the angular frequency, whereas the rest remains barely affected. Thus, the water structure in the aqueous mixture of PE allows the storage of energy elastically through the hydrogen-bond network at large angular frequencies, whereas the mechanical contribution of polyions weakens and fully vanishes at the beginning of shear thinning, explaining the superimposed data with data of bulk water. Our method and findings set the path for future molecular simulations in the nonlinear viscoelastic regime with more complex underlying molecular mechanisms.

Neat and Aqueous Polyelectrolytes under a Steady-Shear Flow.

Materials enabling impact-energy absorption of high-velocity projectiles are of great interest for applications like aerospace. In such a frame, shear thickening fluids were found very useful. Here, we investigated nanorheological properties of neat and aqueous polyelectrolytes of low molecular weights containing poly([2-(methacryloyloxy) ethyl] trimethyl ammonium) as polycations and poly(acrylamide-co-acrylic acid) as polyanions. Results were compared with pure water. We employed nonequilibrium molecular dynamics with the SLLOD algorithm to compute the viscosity at various shear rates. Systems containing polyelectrolytes exhibit shear thickening. The analysis of molecular configurations revealed a strong disruption of the ionic structure and more clusters with smaller sizes on increasing the shear rate. Potential energies showed that shear thickening originates from an increase in intramolecular and van der Waals interactions resulting from the increasing difficultly of polyelectrolyte-based systems to relax at high shear rates. Our method and findings underscore the importance of accounting for the molecular scale in the design of materials absorbing the impact energy efficiently.

Isotope Effect in the Liquid Properties of Water Confined in 100 nm Nanofluidic Channels.

Liquids confined in 10-100 nm spaces show different liquid properties from those in the bulk. Proton transfer plays an essential role in liquid properties. The Grotthuss mechanism, in which charge transfer occurs among neighboring water molecules, is considered to be dominant in bulk water. However, the rotational motion and proton transfer kinetics have not been studied well, which makes further analysis difficult. In this study, an isotope effect was used to study the kinetic effect of rotational motion and proton hopping processes by measurement of the viscosity, proton diffusion coefficient, and the proton hopping activation energy. As a result, a significant isotope effect was observed. These results indicate that the rotational motion is not significant, and the decrease of the proton hopping activation energy enhances the apparent proton diffusion coefficient.

A Theoretical study on the charge and discharge states of Na-ion battery cathode material, Na1+x FePO4 F.

Na2 FePO4 F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na2 FePO4 F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na1.5 FePO4 F has the P21 /b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P21 /b11 Na1.5 FePO4 F and Pbcn Na2 FePO4 F or Na1 FePO4 F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na2 FePO4 F has diffusion paths in the a- and c-axis directions, Na2 FePO4 F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.

First-Principles Study of the Band Diagrams and Schottky-Type Barrier Heights of Aqueous Ta3N5 Interfaces.

The photo(electro)chemical production of hydrogen by water splitting is an efficient and sustainable method for the utilization of solar energy. To improve photo(electro)catalytic activity, a Schottky-type barrier is typically useful to separate excited charge carriers in semiconductor electrodes. Here, we focused on studying the band diagrams and the Schottky-type barrier heights of Ta3N5, which is one of the most promising materials as a photoanode for water splitting. The band alignments of the undoped and n-type Ta3N5 with adsorbents in a vacuum were examined to determine how impurities and adsorbents affect the band positions and Fermi energies. The band edge positions as well as the density of surface states clearly depended on the density of ON impurities in the bulk and surface regions. Finally, the band diagrams of the n-type Ta3N5/water interfaces were calculated with an improved interfacial model to include the effect of electrode potential with explicit water molecules. We observed partial Fermi level pinning in our calculations at the Ta3N5/water interface, which affects the driving force for charge separation.

Redox-responsive molecular helices with highly condensed π-clouds.

Helices have long attracted the attention of chemists, both for their inherent chiral structure and their potential for applications such as the separation of chiral compounds or the construction of molecular machines. As a result of steric forces, polymeric o-phenylenes adopt a tight helical conformation in which the densely packed phenylene units create a highly condensed π-cloud. Here, we show an oligomeric o-phenylene that undergoes a redox-responsive dynamic motion. In solution, the helices undergo a rapid inversion. During crystallization, however, a chiral symmetry-breaking phenomenon is observed in which each crystal contains only one enantiomeric form. Crystals of both handedness are obtained, but in a non-racemic mixture. Furthermore, in solution, the dynamic motion of the helical oligomer is dramatically suppressed by one-electron oxidation. X-ray crystallography of both the neutral and oxidized forms indicated that a hole, generated upon oxidation, is shared by the repeating o-phenylene units. This enables conformational locking of the helix, and represents a long-lasting chiroptical memory.

Hydrogen-bonding dynamics in aqueous solutions of amides and acids: monomer, dimer, trimer, and polymer.

Hydrogen-bonding dynamics in aqueous solutions of series of amides and acids have been investigated by means of femtosecond Raman-induced Kerr effect spectroscopy and ab initio quantum chemistry calculation. The amides and acids studied here are acetamide, 1,3-propanedicarboxamide, 1,3,5-pentanetricaroxamide, polyacrylamide with Mw=1500, acetic acid, 1,3-propanedicarboxylic acid, 1,3,5-pentanetricarboxylic acid, and poly(acrylic acid) with Mw=2000. The femtosecond damped transient feature for aqueous amide solutions, which arises from the intermolecular hydrogen bonds of amide and water, becomes clearer with the larger molecular weight of amide. A characteristic vibrational band at about 100 cm(-1) is assigned as the hydrogen-bonding vibrational mode and the ab initio quantum chemistry calculation result indicates that at least two waters, which make up the hydrogen-bonding network with amide, are necessary for this mode. The hydrogen-bonding vibrational mode at about 100 cm(-1) in aqueous amide solutions shifts to the higher frequency with the larger molecular weight amide in consequence of the stronger intermolecular interaction between amide and water. The evidence likely comes from the stronger hydrophobic interaction for polymer than oligomers and monomer. In the picosecond time region, an extra slow relaxation process with a time constant of about 60 ps has been found in the aqueous polymer solutions. The relaxation is assigned as a local motion of the constitutional repeat unit of polymers from comparison with monomer and oligomers.

Real-time observation of intramolecular proton transfer in the electronic ground state of chloromalonaldehyde: an ab initio study of time-resolved photoelectron spectra.

The authors report on studies of time-resolved photoelectron spectra of intramolecular proton transfer in the ground state of chloromalonaldehyde, employing ab initio photoionization matrix elements and effective potential surfaces of reduced dimensionality, wherein the couplings of proton motion to the other molecular vibrational modes are embedded by averaging over classical trajectories. In the simulations, population is transferred from the vibrational ground state to vibrationally hot wave packets by pumping to an excited electronic state and dumping with a time-delayed pulse. These pump-dump-probe simulations demonstrate that the time-resolved photoelectron spectra track proton transfer in the electronic ground state well and, furthermore, that the geometry dependence of the matrix elements enhances the tracking compared with signals obtained with the Condon approximation. Photoelectron kinetic energy distributions arising from wave packets localized in different basins are also distinguishable and could be understood, as expected, on the basis of the strength of the optical couplings in different regions of the ground state potential surface and the Franck-Condon overlaps of the ground state wave packets with the vibrational eigenstates of the ion potential surface.

Time-resolved photoelectron spectroscopy of proton transfer in the ground state of chloromalonaldehyde: wave-packet dynamics on effective potential surfaces of reduced dimensionality.

We report on a simple but widely useful method for obtaining time-independent potential surfaces of reduced dimensionality wherein the coupling between reaction and substrate modes is embedded by averaging over an ensemble of classical trajectories. While these classically averaged potentials with their reduced dimensionality should be useful whenever a separation between reaction and substrate modes is meaningful, their use brings about significant simplification in studies of time-resolved photoelectron spectra in polyatomic systems where full-dimensional studies of skeletal and photoelectron dynamics can be prohibitive. Here we report on the use of these effective potentials in the studies of dump-probe photoelectron spectra of intramolecular proton transfer in chloromalonaldehyde. In these applications the effective potentials should provide a more realistic description of proton-substrate couplings than the sudden or adiabatic approximations commonly employed in studies of proton transfer. The resulting time-dependent photoelectron signals, obtained here assuming a constant value of the photoelectron matrix element for ionization of the wave packet, are seen to track the proton transfer.

Very fast tunneling in the early stage of reaction dynamics.

The difference between quantum and classical survival probabilities for molecular dissociation dynamics in the time domain, which arises mainly from quantum mechanical tunneling, has interesting characteristics that are not noticed through the counterpart in energy domain. It is shown that the early stage undergoes a fast tunneling, while the later stage is characterized with a long-lasting slow tunneling. The mechanism of this behavior is analyzed in terms of a quasi-semiclassical theory featuring the geometrical distribution of the so-called tunneling points. In particular, the role of dynamical tunneling is discussed as a phenomenon that typifies the time dependence of tunneling dynamics. It is predicted that these tunneling characteristics will be reflected in the isotope effect and should be experimentally observable.

Extended quantization condition for constructive and destructive interferences and trajectories dominating molecular vibrational eigenstates.

The role of destructive quantum interference in semiclassical quantization of molecular vibrational states is studied. This aspect is crucial for correct quantization, since failure in the appropriate treatment of destructive interference quite often results in many spurious peaks and broad background to hide the true peaks. We first study the time-Fourier transform of the autocorrelation function without performing summation over the trajectories. The resultant quantity, the prespectrum which is a function of individual classical trajectories, provides a clear view about how destructive interference among the trajectories should function. It turns out that the prespectrum is oscillatory but never a random noise. On the contrary, it bears a systematic and regular structure, which is sometimes characterized in terms of very sharp and high peaks in the energy space of the sampled classical trajectories. We have found an extended quantization condition that is responsible for generating these peaks in the prespectrum, which we call the prior quantization condition. Integration of the prespectrum over the trajectory space is supposed to give "zero" (practically a small value of the order of the Planck constant) at a noneigenvalue energy, which is actually a materialization of the destructive interference. Besides, certain finite peaks in the prespectrum survive after the integration to form the true spikes (eigenvalues) in the final spectrum, if they satisfy an additional resonance condition. For these resonance components, the prior quantization condition is reduced to the Einstein-Brillouin-Keller quantization condition. Based on these analyses, we propose a rather conventional filtering technique to efficiently handle tedious computation for destructive interference, and numerically verify that it works well even for multidimensional chaotic systems. This filtering technique is further utilized to extract a few trajectories that dominate an eigenstate of molecular vibration.

Time-dependent probability of quantum tunneling in terms of the quasisemiclassical method.

In view of the rapid progress in experiments of the tunneling dynamics in the time domain, we develop a quasisemiclassical method that is aimed at a study of the proton-transfer dynamics in a large system such as tropolone and its interesting derivatives, to which not only full quantum mechanics, but even a standard semiclassical theory is never easy to apply. In our very tractable method for multidimensional systems, the tunneling paths are generated in terms of the generalized classical mechanics, but the quantum phases arising from the action integral, the Maslov index, and the semicalssical amplitude factor as well in the semiclassical kernels are entirely neglected. This approach is called the quasisemiclassical method. One of the technical issues involved in the general semiclassical scheme is how to locate points from which a tunneling path emanates. Hence the studies of such tunneling points and the quasisemiclassical method should be examined collectively. We test several ways of determining the tunneling point, including those already proposed in the literature and a newly proposed one. It is shown numerically that the quasisemiclassical method with an appropriate choice of tunneling points reproduces the full quantum mechanical tunneling probability reasonably well. This case study indicates that the present conventional approach is promising to the study of large systems. The role of tunneling points in the initial process of tunneling is also discussed.