Search

We present the simulation of the photochemical dynamics of cyclobutanone induced by the excitation of the 3 s Rydberg state. For this purpose, we apply the complete active space self-consistent field method together with the spin–orbit multireference configuration interaction singles treatment, combined with the trajectory surface hopping for the inclusion of nonadiabatic effects. The simulations were performed in the spin-adiabatic representation, including nine electronic states derived from three singlet and two triplet spin-diabatic states. Our simulations reproduce the two previously observed primary dissociation channels: the C2 pathway yielding C2H4 + CH2CO and the C3 pathway producing c-C3H6 + CO. In addition, two secondary products, CH2 + CO from the C2 pathway and C3H6 from the C3 pathway, both of them previously reported, are also observed in our simulation. We determine the ratio of the C3:C2 products to be about 2.8. Our findings show that most of the trajectories reach their electronic ground state within 200 fs, with dissociation events finished after 300 fs. We also identify the minimum energy conical intersections that are responsible for the relaxation and provide an analysis of the photochemical reaction mechanism based on multidimensional scaling. Furthermore, we demonstrate a minimal impact of triplet states on the photodissociation mechanism within the observed timescale. In order to provide a direct link to experiments, we simulate the gas phase ultrafast electron diffraction patterns and connect their features to the underlying structural dynamics. [ABSTRACT FROM AUTHOR]

We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels–Alder reaction between furan and ethene as an example. [ABSTRACT FROM AUTHOR]

Understanding the dynamics of polymers in confined environments is pivotal for diverse applications ranging from polymer upcycling to bioseparations. In this study, we develop an entropic barrier model using self-consistent field theory that considers the effect of attractive surface interactions, solvation, and confinement on polymer kinetics. In this model, we consider the translocation of a polymer from one cavity into a second cavity through a single-segment-width nanopore. We find that, for a polymer in a good solvent (i.e., excluded volume, u0 > 0), there is a nonmonotonic dependence of mean translocation time (τ) on surface interaction strength, ɛ. At low ɛ, excluded volume interactions lead to an energetic penalty and longer translocation times. As ɛ increases, the surface interactions counteract the energetic penalty imposed by excluded volume and the polymer translocates faster through the nanopore. However, as ɛ continues to increase, an adsorption transition occurs, which leads to significantly slower kinetics due to the penalty of desorption from the first cavity. The ɛ at which this adsorption transition occurs is a function of the excluded volume, with higher u0 leading to an adsorption transition at higher ɛ. Finally, we consider the effect of translocation across different size cavities. We find that the kinetics for translocation into a smaller cavity speeds up while translocation to a larger cavity slows down with increasing ɛ due to higher surface contact under stronger confinement. [ABSTRACT FROM AUTHOR]

Block copolymer self-assembly in conjunction with nonsolvent-induced phase separation (SNIPS) has been increasingly leveraged to fabricate integral-asymmetric membranes. The large number of formulation and processing parameters associated with SNIPS, however, has prevented the reliable construction of high performance membranes. In this study, we apply dynamical self-consistent field theory to model the SNIPS process and investigate the effect of various parameters on the membrane morphology: solvent selectivity, nonsolvent selectivity, initial film composition, and glass transition composition. We examine how solvent selectivity and concentration of polymers in the film impact the structure of micelles that connect to form the membrane matrix. In particular, we find that preserving the order in the surface layer and forming a connection between the supporting and surface layer are nontrivial and sensitive to each parameter studied. The effect of each parameter is discussed, and suggestions are made for successfully fabricating viable block copolymer membranes. [ABSTRACT FROM AUTHOR]

The phase behavior of binary blends of AB diblock copolymers of compositions f and 1 − f is examined using field-theoretic simulations. Highly asymmetric compositions (i.e., f ≈ 0) behave like homopolymer blends macrophase separating into coexisting A- and B-rich phases as the segregation is increased, whereas more symmetric diblocks (i.e., f ≈ 0.5) microphase separate into an ordered lamellar phase. In self-consistent field theory, these behaviors are separated by a Lifshitz critical point at f = 0.2113. However, its lower critical dimension is believed to be four, which implies that the Lifshitz point should be destroyed by fluctuations. Consistent with this, it is found to transform into a tricritical point. Furthermore, the highly swollen lamellar phase near the mean-field Lifshitz point disorders into a bicontinuous microemulsion (BμE), consisting of large interpenetrating A- and B-rich microdomains. BμE has been previously reported in ternary blends of AB diblock copolymer with its parent A- and B-type homopolymers, but in that system the homopolymers have a tendency to macrophase separate. Our alternative system for creating BμE is free of this macrophase separation. [ABSTRACT FROM AUTHOR]

Despite the wide existence of vesicles in living cells as well as their important applications like drug delivery, the underlying mechanism of vesicle fusion/fission remains under debate. Classical models cannot fully explain recent observations in experiments and simulations. Here, we develop a constrained self-consistent field theory that allows tracking the shape evolution and free energy as a function of center-of-mass separation distance. Fusion and fission are described in a unified framework. Both the kinetic pathway and the mechanical response can be simultaneously captured. By taking vesicles formed by polyelectrolytes as a model system, we predict discontinuous transitions between the three morphologies: parent vesicle with a single cavity, hemifission/hemifusion, and two separated child vesicles, as a result of breaking topological isomorphism. With the increase in inter-vesicle repulsion, we observe a great reduction in the cleavage energy, indicating that vesicle fission can be achieved without hemifission, in good agreement with simulation results. The force–extension relationship elucidates typical plasticity for separating two vesicles. The super extensibility in the mechanical response of vesicle is in stark contrast to soft particles with other morphologies, such as cylinder and sphere. Our work elucidates the fundamental physical chemistry based on intrinsic topological features of vesicle fusion/fission, which provides insights into various phenomena observed in experiments and simulations. [ABSTRACT FROM AUTHOR]

Pathways of two-body fragmentation of BrCNq+ (q = 2, 3) have been explored by combined experimental and theoretical studies. In the experiment, the BrCN molecule is ionized by 1 keV electron impact and the created fragment ions are detected using an ion momentum imaging spectrometer. Six two-body fragmentation channels are identified. By measuring the momentum vectors of the fragment ions, the kinetic energy release (KER) distributions for these channels have been determined. Theoretically, the potential energy curves of BrCNq+ (q = 2, 3) as a function of Br–C and C–N internuclear distances are calculated by the complete active space self-consistent field method. By comparing the measured KER and theoretical predictions, pathways for the fragmentation channels are assigned. The relative branching ratios of the channels are also determined. [ABSTRACT FROM AUTHOR]

Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a "sign problem" associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a "partial saddle point approximation" (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems. [ABSTRACT FROM AUTHOR]

The first finite basis set implementation of the real-time time-dependent self-consistent field method in a dynamic (time-dependent) magnetic field using London atomic orbitals (LAOs) is presented. The accuracy of the finite basis approach using LAOs is benchmarked against numerical results from the literature for the hydrogen atom and H2 in the presence of rapidly oscillating magnetic fields. This comparison is used to inform the choice of appropriate basis sets for studies under such conditions. Remarkably, relatively modest compact LAO basis sets are sufficient to obtain accurate results. Analysis of electron dynamics in the hydrogen atom shows that LAO calculations correctly capture the time evolution of orbital occupations. The Fourier transformation of the autocorrelation function yields a power spectrum exhibiting harmonics associated with coherent emission, which closely matches the literature and further confirms the accuracy of this approach. The dynamical response of the electron density in H2 for a magnetic field parallel to the internuclear axis shows similar behavior to benchmark studies. The flexibility of this implementation is then demonstrated by considering how the dynamical response changes as a function of the orientation of the molecule relative to the applied field. At non-parallel orientations, the symmetry of the system is lowered and numerical benchmark data, which exploit cylindrical symmetry, are no-longer readily available. The present study demonstrates the utility of LAO-based calculations for extreme dynamic magnetic fields, providing a stress-test on the choice of basis. Future applications of this approach for less extreme dynamic magnetic fields are briefly discussed. [ABSTRACT FROM AUTHOR]

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick's theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence. [ABSTRACT FROM AUTHOR]

10.1063/5.0102669 5 A. Szulc, M. Mungan, and I. Regev, "Cooperative effects driving the multi-periodic dynamics of cyclically sheared amorphous solids", J. Chem. Phys. In the work of Szulc I et al. i ,[5] the authors consider an intriguing aspect of cyclically sheared amorphous solids, which is the presence of limit cycles that are multi-periodic. Due to this "melt memory", the crystallization temperature is higher than in the preceding cycle of crystallization. A polymer melt obtained by heating from its semicrystalline state undergoes crystallization at a higher crystallization temperature upon cooling than in the previous cycle of crystallization. [Extracted from the article]

A computational framework that leverages data from self-consistent field theory simulations with deep learning to accelerate the exploration of parameter space for block copolymers is presented. This is a substantial two-dimensional extension of the framework introduced in the work of Xuan et al. [J. Comput. Phys. 443, 110519 (2021)]. Several innovations and improvements are proposed. (1) A Sobolev space-trained, convolutional neural network is employed to handle the exponential dimension increase of the discretized, local average monomer density fields and to strongly enforce both spatial translation and rotation invariance of the predicted, field-theoretic intensive Hamiltonian. (2) A generative adversarial network (GAN) is introduced to efficiently and accurately predict saddle point, local average monomer density fields without resorting to gradient descent methods that employ the training set. This GAN approach yields important savings of both memory and computational cost. (3) The proposed machine learning framework is successfully applied to 2D cell size optimization as a clear illustration of its broad potential to accelerate the exploration of parameter space for discovering polymer nanostructures. Extensions to three-dimensional phase discovery appear to be feasible. [ABSTRACT FROM AUTHOR]

Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction‐based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self‐consistent field Theory (VSCF) and post‐VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode‐coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode‐mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI‐PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post‐VSCF methods to generate more accurate representations of computed infrared spectra. [ABSTRACT FROM AUTHOR]

A model potential previously developed for the ammonia molecule is treated in a single-center partial-wave approximation in analogy with a self-consistent field method developed by Moccia. The latter was used in a number of collision studies. The model potential is used to calculate DC Stark resonance parameters, i.e., resonance positions and shifts using the exterior complex scaling method for the radial coordinate. Three molecular valence orbitals are investigated for fields along the three Cartesian coordinates, i.e., along the molecular axis and in two perpendicular directions. The work extends previous work on the planar-geometry water molecule for which non-monotonic shifts were observed. We find such non-monotonic shifts for fields along the molecular axis. For perpendicular fields, we report the splitting of the 1 e orbitals into a fast- and a slow-ionizing orbital. [ABSTRACT FROM AUTHOR]

A self-consistent field method based on a density matrix functional theory scheme is presented to compute the potential energy curve of the hydrogen molecule. A functional to recover the so-called cumulant correlation energy is introduced. The form of such functional is very simple being dependent explicitly only on the square of the overlap matrix element between the unrestricted non-orthogonal occupied orbitals of the calculation. The functional is a sum of two contributions referring to non-dynamic and dynamic correlation. The agreement with the full-CI potential energy curve is within chemical accuracy. [ABSTRACT FROM AUTHOR]

The structure and scaling properties of inwardly curved polymer brushes, tethered under good solvent conditions to the inner surface of spherical shells such as membranes and vesicles, are studied by extensive molecular dynamics simulations and compared with earlier scaling and self-consistent field theory predictions for different molecular weights of the polymer chains N and grafting densities σg in the case of strong surface curvature, R−1. We examine the variation of the critical radius R*(σg), separating the regimes of weak concave brushes and compressed brushes, predicted earlier by Manghi et al. [Eur. Phys. J. E 5, 519–530 (2001)], as well as various structural properties such as the radial monomer- and chain-end density profiles, orientation of bonds, and brush thickness. The impact of chain stiffness, κ, on concave brush conformations is briefly considered as well. Eventually, we present the radial profiles of the local pressure normal, PN, and tangential, PT, to the grafting surface, and the surface tension γ(σg), for soft and rigid brushes, and find a new scaling relationship P N (R) ∝ σ g 4 , independent of the degree of chain stiffness. [ABSTRACT FROM AUTHOR]

An alternative approach to density functional theory based on self-consistent field theory for ring polymers is applied to neutral atoms hydrogen to neon in their ground-states. The spontaneous emergence of an atomic shell structure and spherical symmetry-breaking of the total electron density are predicted by the model using the ideas of polymer excluded-volume between pairs of electrons to enforce the Pauli-exclusion principle and an exact electron self-interaction correction. The Pauli potential is approximated by neglecting inter-atomic correlations along with other types of correlations, and comparisons to Hartree–Fock theory are made, which also ignores correlations. The model shows excellent agreement with Hartree–Fock theory to within the standards of orbital-free density functional theory for the atomic binding energies and density profiles of the first six elements, providing exact matches for the elements hydrogen and helium. The predicted shell structure starts to deviate significantly past the element neon, and spherical symmetry-breaking is first predicted to occur at carbon instead of boron. The self-consistent field theory energy functional that describes the model is decomposed into thermodynamic components to trace the origin of spherical symmetry-breaking. It is found to arise from the electron density approaching closer to the nucleus in non-spherical distributions, which lowers the energy despite resulting in frustration between the quantum kinetic energy, electron–electron interaction, and the Pauli exclusion interaction. The symmetry-breaking effect is found to have a minimal impact on the binding energies, which suggests that the spherical-averaging approximation used in previous work is physically reasonable when investigating atomic systems. The pair density contour plots display behavior similar to polymer macro-phase separation, where individual electron pairs occupy single lobe structures that together form a dumbbell shape analogous to the 2p orbital shape. It is further shown that the predicted densities satisfy known constraints and produce the same total electronic density profile that is predicted by other formulations of quantum mechanics. [ABSTRACT FROM AUTHOR]

Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS. [ABSTRACT FROM AUTHOR]

This study examines the ultraviolet (UV) divergence in field-theoretic simulations (FTSs) of block copolymer melts, which causes an unphysical dependence on the grid resolution, Δ, used to represent the fields. Our FTSs use the discrete Gaussian–chain model and a partial saddle-point approximation to enforce incompressibility. Previous work has demonstrated that the UV divergence can be accounted for by defining an effective interaction parameter, χ = z ∞ χ b + c 2 χ b 2 + c 3 χ b 3 + ⋯ , in terms of the bare interaction parameter, χb, used in the FTSs, where the coefficients of the expansion are determined by a Morse calibration. However, the need to use different grid resolutions for different ordered phases generally restricts the calibration to the linear approximation, χ ≈ z∞χb, and prevents the calculation of order–order transitions. Here, we resolve these two issues by showing how the nonlinear calibration can be translated between different grids and how the UV divergence can be removed from free energy calculations. By doing so, we confirm previous observations from particle-based simulations. In particular, we show that the free energy closely matches self-consistent field theory (SCFT) predictions, even in the region where fluctuations disorder the periodic morphologies, and similarly, the periods of the ordered phases match SCFT predictions, provided the SCFT is evaluated with the nonlinear χ. [ABSTRACT FROM AUTHOR]

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model's predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane. [ABSTRACT FROM AUTHOR]

A synchrotron-based photoionization spectrum up to 27 eV represents a considerable improvement in resolution over early He(I) and He(II) spectra. Symmetry-adapted coupled cluster calculations of the ionic state sequence give the sequence of state vertical ionization energies (VIE) as 12B2 < 12B1 < 12A2 < 22B1 < 12A1. Generally, these symmetry-adapted cluster configuration interactions VIE match reasonably well with the experimental spectrum over this wide energy range. Density functional calculations of the corresponding adiabatic terms (AIE) were also performed. Higher energy ionic states were determined by complete active space self-consistent field methods; these include all π-ionizations and some σ-ionic states. These were analyzed by Franck–Condon (FC) procedures and compared with an experiment. The spectral onset is complex, where two states, later shown to be the 12B2 and 12B1 states, are strongly overlapping. The superposition of the FC vibrational structure in the 12B2 and 12B1 states accounts for most of the peaks arising at the onset of the photoelectron spectra. However, the small separation between these two ionic states makes vibronic interaction fairly inevitable. In the absence of Herzberg–Teller analyses for ionic states, we have sought and determined a transition state between the 12B2 and 12B1 states, showing that vibronic coupling does occur. The lack of degradation in the vibrational envelope of the higher of the two states contrasts with our previous work on the halogenobenzenes, where overlapping state envelopes led to considerable widening of the line width at half-height of the higher energy states. [ABSTRACT FROM AUTHOR]

We study the effect of chain rigidity on tailoring the nanoparticle locations for neutral and selective particles embedded in the lamellar morphology formed by semiflexible diblock copolymer chains using self-consistent field calculations. The nanoparticles are modeled through a cavity function, and the semiflexible chains are represented by the continuous Kratsky-Porod chain model. In general situation, the nanoparticles prefer to stay at the interface in order to reduce the interface areas and thus the system free energy. However, the particle distribution at the domain center is subtle, and the underlying physics is intrinsically different depending on the polymer flexibility. In the case of flexible chains, the entropy just contributes a constant shift to the free energy when the nanoparticles move around the domain center indicating that the local metastable state if appears at the domain center is wholly attributed to the local minimum in the enthalpy. If the polymers are rigid, the variation of the particle distribution at the domain center has a close relation with the polymer rigidity and nanoparticle size. In the case of strongly rigid polymers with small nanoparticles, a nearly uniform particle distribution at the domain center is observed, while in other cases, a local enhancement of particle distribution there is found. In contrast to the case of flexible chains, further analysis reveals the crucial role of entropy in controlling the shape of particle distributions at the phase domain. Specifically, the local metastable state appears in the domain center is determined by the large entropy there which arises from the weak coupling of bond orientations that allows the polymer chains to be relatively relaxed. When the particle becomes selective, its distribution in the phase domain exhibits a shift almost uniformly rather than changes its profile, and the underlying physics still holds. In all, our study establishes a strong coupling between the chain rigidity and effect of entropy. [ABSTRACT FROM AUTHOR]

The effect of colloidal nanoparticles on the phase changes of the amphiphilic AB linear diblock, A1A2B, and A2B heteroarm star copolymers confined between two polymer brush substrates was investigated by using a real-space self-consistent field theory. By changing the concentrations of nanoparticles and polymer brushes, the phase structure of the amphiphilic AB copolymer transforms from lamellar to core-shell hexagonal phase to cylinder phase. The pattern of A2B heteroarm star copolymer changes from core-shell hexagonal phases to lamellar phases and the layer decreases when increasing the density of the polymer brushes. The results showed that the phase behavior of the system is strongly influenced by the polymer brush architecture and the colloidal nanoparticle numbers. The colloidal nanoparticles and the soft confined surface of polymer brushes make amphiphilic AB copolymers easier to form ordered structures. The dispersion of the nanoparticles was also investigated in detail. The soft surfaces of polymer brushes and the conformation of the block copolymers work together to force the nanoparticles to disperse evenly. It will give helpful guidance for making some new functional materials by nano etching technology, nano photoresist, and nanoprinting. [ABSTRACT FROM AUTHOR]

The unrestricted Hartree‐Fock method is extended to correlation calculation within the density‐matrix functional theory. The method is derived from an entropic cumulant functional for the correlation energy. The eigenvalue equations for the spin‐orbitals are modified by the orbital occupation numbers. The Euler equation for the occupation numbers results in the Fermi‐Dirac distribution, which is very efficient to update as soon as the orbital eigenvalue equations are solved. The method is demonstrated on the ground state of O2. [ABSTRACT FROM AUTHOR]

We constructed phase diagrams of conformationally asymmetric diblock copolymer A-B melts using the polymer self-consistent field (SCF) calculations of both the dissipative particle dynamics chain (DPDC) model (i.e., compressible melts of discrete Gaussian chains with the DPD non-bonded potential) and the "standard" model (i.e., incompressible melts of continuous Gaussian chains with the Dirac δ-function non-bonded potential) in the χN-ε plane, where χN and ε characterize, respectively, the repulsion and conformational asymmetry between the A and B blocks, at the A-block volume fraction f = 0.2 and 0.3. Consistent with previous SCF calculations of the "standard" model, σ and A15 are the only stable Frank–Kasper (FK) phases among the five FK (i.e., σ, A15, C14, C15 and Z) phases considered. The stability of σ and A15 is due to their delicate balance between the energetic and entropic contributions to the Helmholtz free energy per chain of the system, which, within our parameter range, increases in the order of σ/A15, Z, and C14/C15. While in general the SCF phase diagrams of these two models are qualitatively consistent, A15 is not stable for the DPDC model at the copolymer chain length N = 10 and f = 0.3; any differences in the SCF phase diagrams are solely due to the differences between these two models. [ABSTRACT FROM AUTHOR]

Aiming at the difficult problem of solving the conformation statistics of complex polymers, this study presents a novel and concise conformation statistics theoretical approach based on Monte Carlo and Neural Network method. This method offers a new research idea for investigating the conformation statistics of complex polymers, characterized by its simplicity and practicality. It can be applied to more complex topological structure, more higher degree of freedom polymer systems with higher dimensions, theory research on dynamic self-consistent field theory and polymer field theory, as well as the analysis of scattering experimental data. The conformation statistics of complex polymers determine the structure and response properties of the system. Using the new method proposed in this study, taking the semiflexible ring diblock copolymer as an example, Monte Carlo simulation is used to sample this ring conformation to construct the dataset of polymer. The structure factor describing conformation statistics are expressed as continuous functions of structure parameters by neural network supervised learning. This is the innovation of this work. As an application, the structure factors represented by neural networks were introduced into the random phase approximation theory to study the microphase separation of semiflexible ring diblock copolymers. The influence of the ring's topological properties on the phase transition behavior was pointed out. [ABSTRACT FROM AUTHOR]

The structure of a thin film of symmetric ABA triblock copolymer melts in an external in‐plane DC or AC electric field is studied theoretically. The situation is considered when the triblock copolymer forms a hexagonal morphology of standing cylinders in bulk in the absence of an external field. Self‐consistent field theory calculations are carried out to determine the most thermodynamically favorable thin film structure among the cylindrical phases perpendicular and parallel to the substrate and the lamellar phase perpendicular to the substrate. The results are presented as phase diagrams with the film thickness and electric field energy on the axes and as distributions of the local composition, which serve as an order parameter in the system. It is confirmed that electric fields only weakly affect the spinodal curves of block copolymers but they can reorient or markedly modify microphase‐separated morphologies in those systems. Such restructuring is consistent with a considerable modulation of the free surface of a copolymer film and it can lead to the coexistence of different phases that appear in the film areas with different local thicknesses. [ABSTRACT FROM AUTHOR]

Recently, we proposed a new orbital analysis method, natural reaction orbital (NRO), which automatically extracts orbital pairs that characterize electron transfer in reaction processes by singular value decomposition of the first-order orbital response matrix to the nuclear coordinate displacements [Ebisawa et al., Phys. Chem. Chem. Phys. 24, 3532 (2022)]. NRO analysis along the intrinsic reaction coordinate (IRC) for several typical chemical reactions demonstrated that electron transfer occurs mainly in the vicinity of transition states and in regions where the energy profile along the IRC shows shoulder features, allowing the reaction mechanism to be explained in terms of electron motion. However, its application has been limited to single configuration theories such as Hartree–Fock theory and density functional theory. In this work, the concept of NRO is extended to multiconfigurational wavefunctions and formulated as the multiconfiguration NRO (MC-NRO). The MC-NRO method is applicable to various types of electronic structure theories, including multiconfigurational theory and linear response theory, and is expected to be a practical tool for extracting the essential qualitative features of a broad range of chemical reactions, including covalent bond dissociation and chemical reactions in electronically excited states. In this paper, we calculate the IRC for five basic chemical reaction processes at the level of the complete active space self-consistent field theory and discuss the phenomenon of electron transfer by performing MC-NRO analysis along each IRC. Finally, issues and future prospects of the MC-NRO method are discussed. [ABSTRACT FROM AUTHOR]

This paper combines the valence bond block diabatization approach with the idea of orbital breathing. With highly compact wave functions, the breathing orbital valence bond (BOVB) method is applied to investigate several atomic and molecular properties, including the electron affinity of F, the adiabatic and diabatic potential energy curves and the dipole moment curves of the two lowest-lying 1Σ+ states, the electronic coupling curve and the crossing distance of the two diabatic states, and the spectroscopic constants of the ground states for LiF. The configuration selection scheme proposed in this work is quite general, requiring only the selection of several de-excitation and excitation orbitals in a sense like the restricted active space self-consistent field method. Practically, this is also the first time that BOVB results are extrapolated to complete basis set limit. Armed with the chemical intuition provided by valence bond theory, the classic but challenging covalent-ionic interaction in the title molecule is not only conceptually interpreted but also accurately computed. [ABSTRACT FROM AUTHOR]

We examine athermal binary blends composed of conformationally asymmetric polymers of equal molecular volume next to a surface of width ξ. The self-consistent field theory (SCFT) of Gaussian chains predicts that the more compact polymer with the shorter average end-to-end length, R0, is entropically favored at the surface. Here, we extend the SCFT to worm-like chains with small persistence lengths, ℓp, relative to their contour lengths, ℓc, for which R 0 ≈ 2 ℓ p ℓ c . In the limit of ℓp ≪ ξ, we recover the Gaussian-chain prediction where the segregation depends only on the product ℓpℓc, but for realistic polymer/air surfaces with ξ ∼ ℓp, the segregation depends separately on the two quantities. Although the surface continues to favor flexible polymers with smaller ℓp and bulky polymers with shorter ℓc, the effect of bulkiness is more pronounced. This imbalance can, under specific conditions, lead to anomalous surface segregation of the more extended polymer. For this to happen, the polymer must be bulkier and stiffer, with a stiffness that is sufficient to produce a larger R0 yet not so rigid as to reverse the surface affinity that favors bulky polymers. [ABSTRACT FROM AUTHOR]

Block polymers are an attractive platform for uncovering the factors that give rise to selfassembly in soft matter owing to their relatively simple thermodynamic description, as captured in self-consistent field theory (SCFT). SCFT historically has found great success explaining experimental data, allowing one to construct phase diagrams from a set of candidate phases, and there is now strong interest in deploying SCFT as a screening tool to guide experimental design. However, using SCFT for phase discovery leads to a conundrum: How does one discover a new morphology if the set of candidate phases needs to be specified in advance? This long-standing challenge was surmounted by training a deep convolutional generative adversarial network (GAN) with trajectories from converged SCFT solutions, and then deploying the GAN to generate input fields for subsequent SCFT calculations. The power of this approach is demonstrated for network phase formation in neat diblock copolymer melts via SCFT. A training set of only five networks produced 349 candidate phases spanning known and previously unexplored morphologies, including a chiral network. This computational pipeline, constructed here entirely from open-source codes, should find widespread application in block polymer phase discovery and other forms of soft matter. [ABSTRACT FROM AUTHOR]

The potential energy curves for molecular ions up to N 2 4 + are calculated in an ab initio manner using the multi configurational self-consistent field method. Specifically, we implement in an automatic way a previously used double loop optimisation scheme within the multi configurational self-consisted field method. We obtain the potential energy curves up to N 2 4 + ions with any combination of core, inner valence, and outer valence holes. Finally, we provide the code used to generate these potential energy curves. [ABSTRACT FROM AUTHOR]

The theory of special relativity is a topic not addressed in theoretical chemistry courses. Einstein's theory of relativity is mainly applied to high‐energy phenomena governed by physical laws not used to describe the chemical bond. Nevertheless, the electrons of heavy elements orbit around the nucleus with relativistic velocities. Here it is, then, the laws that appeared so unusual for chemistry became fundamental tools for investigating the properties of these atoms. Studying the electronic structure of heavy elements is the best place for chemists to learn to tame the pitfalls of a little‐known subject, such as Einstein's theory, and become familiar with its elegant formalism. In this paper, the main relativistic approaches to quantum chemistry are revisited. The problem of the finite nuclear size, the approximation methods used to solve the relativistic equations for k‐shell electrons, and the problem of the negative energy continuum are investigated. Some ansatzes to solve divergence's problem of relativistic corrections near the nucleus are proposed. Finally, a possible connection between the self‐consistent field method and bondonic chemistry is also investigated. [ABSTRACT FROM AUTHOR]

Using self-consistent field theory, we study the effect of reversible cross-linking on the nucleation behavior of a binary polymer blend where only one of the components is able to form cross-links. To control the total number of cross-links and their distribution, we introduce a position-dependent cross-linking probability function that is characterized mainly by two parameters, the magnitude and the width. In the weakly cross-linked region, where the product of the magnitude and width, I, is small, the nucleation behavior is classical-like and the profile of the free energy excess is unimodal. In contrast, in the strongly cross-linked region, the profile of the free energy excess becomes bimodal, and the free energy minimum specifies a metastable nucleus. In a certain I, the free energy barrier for the metastable nucleus turns to be negative, which means it becomes more stable. In both cases, the free energy barrier of the critical nucleus is lower than that without cross-linking, indicating that cross-linking always facilitates nucleation although the dynamic behavior may be different when a metastable nucleus is involved during the nucleation process. The free energy analysis demonstrates that the interaction energy rather than the entropy is responsible for the properties of the critical nucleus. Our study provides an easy alternative way for the control of the nucleation behavior and may attract practical interest. [ABSTRACT FROM AUTHOR]

Using analytical considerations and particle-based simulations of a coarse-grained model, we study the relaxation of a density modulation in a polymer system without nonbonded interactions. We demonstrate that shallow density modulations with identical amplitudes and wavevectors that have been prepared by different processes exhibit different nonexponential decay behaviors. Thus, in contrast to the popular assumption of dynamic self-consistent field theory, the density alone does not suffice to characterize the configuration of the polymer system. We provide an analytic description within Linear-Response Theory (LRT) and the Rouse model that quantitatively agree with the results of the particle-based simulations. LRT is equivalent to a generalized model-B dynamics with an Onsager coefficient that is nonlocal in space and time. Alternatively, the Rouse description can be cast into a dynamic density-functional theory that uses the full probability distribution of single-chain configurations as a dynamic variable and yields a memory-free description of the dynamics that quantitatively accounts for the dependence on the preparation process. An approximate scheme that only considers the joint distribution of the first two Rouse modes—the ellipsoid model—is also explored. [ABSTRACT FROM AUTHOR]

Many macromolecules of biological and technological interest are both chiral and semi-flexible. DNA and collagen are good examples. Such molecules often form chiral nematic (or cholesteric) phases, as is well-documented in collagen and chitin. This work presents a method for studying cholesteric phases in the highly successful self-consistent field theory of worm-like chains, offering a new way of studying many biologically relevant molecules. The method involves an effective Hamiltonian with a chiral term inspired by the Oseen–Frank (OF) model of liquid crystals. This method is then used to examine the formation of cholesteric phases in chiral-nematic worm-like chains as a function of polymer flexibility, as well as the optimal cholesteric pitch and distribution of polymer segment orientations. Our approach not only allows for the determination of the isotropic–cholesteric transition and segment distributions, beyond what the OF model promises, but also explicitly incorporates polymer flexibility into the study of the cholesteric phase, offering a more complete understanding of the behavior of semiflexible chiral-nematic polymers. [ABSTRACT FROM AUTHOR]

The self-assembly of block copolymer melts and solutions with two-dimensional density inhomogeneity is studied using modified inhomogeneous statistical associating fluid theory (iSAFT). A real-space combinatorial screening method under density functional theory formalism is proposed and used to map out the phase diagram of block copolymer melts including order–disorder transitions and order–order transitions. The predicted phase diagram agrees well with molecular dynamics simulation and self-consistent field theory. The compressibility effect on order–disorder transition temperature for block copolymer melts is modeled using iSAFT. The pressure induced temperature change by theory has a similar trend to experimental studies. Then, the lyotropic and thermotropic self-assembly phase behavior of block copolymer solutions is investigated. Detailed density distributions by iSAFT provide insight into the lyotropic properties of the block copolymer solutions at the molecular level. The effect of the block copolymer molecular architecture is studied by comparing block copolymers with different molecular packing parameters. Block copolymer solutions in the inverted hexagonal phase are predicted by theory for the block copolymer having a large molecular packing parameter. Finally, solvent selectivity is studied by modeling the block copolymers in a neutral good solvent. The enhanced local solvent concentration predicted by theory explains the reason for fewer ordered phases found in experiments. [ABSTRACT FROM AUTHOR]

Real-time time-dependent density functional theory (RT-TDDFT) is an attractive tool to model quantum dynamics by real-time propagation without the linear response approximation. Sharing the same technical framework of RT-TDDFT, imaginary-time time-dependent density functional theory (it-TDDFT) is a recently developed robust-convergence ground state method. Presented here are high-precision all-electron RT-TDDFT and it-TDDFT implementations within a numerical atom-centered orbital (NAO) basis function framework in the FHI-aims code. We discuss the theoretical background and technical choices in our implementation. First, RT-TDDFT results are validated against linear-response TDDFT results. Specifically, we analyze the NAO basis sets' convergence for Thiel's test set of small molecules and confirm the importance of the augmentation basis functions for adequate convergence. Adopting a velocity-gauge formalism, we next demonstrate applications for systems with periodic boundary conditions. Taking advantage of the all-electron full-potential implementation, we present applications for core level spectra. For it-TDDFT, we confirm that within the all-electron NAO formalism, it-TDDFT can successfully converge systems that are difficult to converge in the standard self-consistent field method. We finally benchmark our implementation for systems up to ∼500 atoms. The implementation exhibits almost linear weak and strong scaling behavior. [ABSTRACT FROM AUTHOR]

Polymer-mediated colloidal interactions control the stability and phase properties of colloid–polymer mixtures that are critical for a wide range of important applications. In this work, we develop a versatile self-consistent field theory (SCFT) approach to study this type of interaction based on a continuum confined polymer solution model with explicit solvent and confining walls. The model is formulated in the grand canonical ensemble, and the potential of mean force for the polymer-mediated interaction is computed from grand potentials. We focus on the case of non-adsorbing linear polymers and present a systematic investigation on depletion effects using SCFT. The properties of confined polymer solutions are probed, and mean-field profiles of induced interactions are shown across different physical regimes. We expose a detailed parametric dependence of the interaction, concerning both attractive and repulsive parts, on polymer concentration, chain length, and solvent quality and explore the effect of wall surface roughness, demonstrating the versatility of the proposed approach. Our findings show good agreement with previous numerical studies and experiments, yet extend prior work to new regimes. Moreover, the mechanisms of depletion attraction and repulsion, along with the influence of individual control factors, are further discussed. We anticipate that this study will provide useful insights into depletion forces and can be readily extended to examine more complex colloid–polymer mixtures. [ABSTRACT FROM AUTHOR]

Proteins achieve their functional, active, and operative three dimensional native structures by overcoming the possibility of being trapped in non-native energy minima present in the energy landscape. The enormous and intricate interactions that play an important role in protein folding also determine the stability of the proteins. The large number of stabilizing/destabilizing interactions makes proteins to be only marginally stable as compared to the other competing structures. Therefore, there are some possibilities that they become trapped in the non-native conformations and thus get misfolded. These misfolded proteins lead to several debilitating diseases. This work performs a comparative study of some existing foldability criteria in the computational design of misfold resistant protein sequences based on self-consistent mean field theory. The foldability criteria selected for this study are Ef, Δ, and Φ that are commonly used in protein design procedures to determine the most efficient foldability criterion for the design of misfolding resistant proteins. The results suggest that the foldability criterion Δ is significantly better in designing a funnel energy landscape stabilizing the target state. The results also suggest that inclusion of negative design features is important for designing misfolding resistant proteins, but more information about the non-native conformations in terms of Φ leads to worse results compared to even simple positive design. The sequences designed using Δ show better resistance to misfolding in the Monte Carlo simulations performed in the study. [ABSTRACT FROM AUTHOR]

As shown theoretically earlier via both weak segregation and self-consistent field theories, ordering of confined molten di- and tri-block copolymer morphologies in the presence of a proper 1D patterned substrate could induce the formation of 3D bicontinuous (in particular, diamond-like) morphologies (DLMs). The purpose of the present paper is to study, unlike the previous studies, how the stable DLMs are formed not in a melt but in a solution of symmetric diblock copolymers with a nonselective solvent that wets the thin film on the patterned substrate. It is shown, via a straightforward self-consistent field calculation of the total solution free energy for various competing phases, that the DLM could be formed in the solutions (with the solvent volume fraction of 0.5), which provides much faster thermodynamic equilibration of the solution as compared to the melt. The last circumstance can ease the production of stable DLMs in thin films of copolymers. The phase diagram describing the stable phases on the plane "the pattern period—the film thickness" is built. [ABSTRACT FROM AUTHOR]

[Display omitted] Hypothesis The assembly of core-corona supraparticles in aqueous dispersions has been regularly assisted by auxiliary monomers/oligomers which modify the individual particles with, e.g. , surface grafting of polyethylene oxide (PEO) chains or other hydrophilic monomers. However, this modification complicates the preparation and purification procedures and increases potential upscaling efforts. Hybrid polymer-silica core-corona supracolloids could be more simply assembled if the PEO chains from surfactants, typically used by default as polymer stabilizers, concomitantly act as assembly promotors. The supracolloids assembly could therefore be more easily achieved without requiring particles functionalization or post-purification steps. Methods The self-assembly of supracolloidal particles prepared with PEO-surfactant stabilized (Triton X-405) and/or PEO-grafted polymer particles is compared to differentiate the roles of the PEO chains in the assembly of core-corona supraparticles. Using time-resolved dynamic light scattering (DLS) and cryogenic transmission electron microscopy(cryo-TEM), the effect of concentration of PEO chains (from surfactant) on the kinetics and dynamics of supracolloids assembly is investigated. Self-consistent field (SCF) lattice theory was used to numerically study the distribution of PEO chains at the interfaces present in the supracolloidal dispersions. Findings The PEO based surfactant can be used as assembly promoter of core-corona hybrid supracolloids due to its amphiphilic nature and via establishing hydrophobic interactions. The concentration of the PEO surfactant, and especially the PEO chains distribution over the different interfaces, crucially affect the supracolloids assembly. A simplified pathway for preparing hybrid supracolloidal particles with a well-controlled corona coverage over polymer cores is presented. [ABSTRACT FROM AUTHOR]

The article discusses the use of attosecond and strong-field spectroscopy to investigate correlated electronic and nuclear dynamics in atoms and molecules. It emphasizes the importance of isotopic substitution to study the effects of nuclear motion on molecular responses and provides experimental data on photoelectron spectrograms generated in methane and deuteromethane, revealing differences in the oscillations of photoelectron yield between the two molecules.

The kinetic paths of structural evolution and formation of block copolymer (BCP) particles are explored using dynamic self‐consistent field theory (DSCFT). It is shown that the process‐directed self‐assembly of BCP immersed in a poor solvent leads to the formation of striped ellipsoids, onion‐like particles and double‐spiral lamellar particles. The theory predicts a reversible path of shape transition between onion‐like particles and striped ellipsoidal ones by regulating the temperature (related to the Flory–Huggins parameter between the two components of BCP, χAB) and the selectivity of solvent toward one of the two BCP components. Furthermore, a kinetic path of shape transition from onion‐like particles to double‐spiral lamellar particles, and then back to onion‐like particles is demonstrated. By investigating the inner‐structural evolution of a BCP particle, it is identified that changing the intermediate bi‐continuous structure into a layered one is crucial for the formation of striped ellipsoidal particles. Another interesting finding is that the formation of onion‐like particles is characterized by a two‐stage microphase separation. The first is induced by the solvent preference, and the second is controlled by the thermodynamics. The findings lead to an effective way of tailoring nanostructure of BCP particles for various industrial applications. [ABSTRACT FROM AUTHOR]

Combining Brownian dynamics simulations and self-consistent field theory, we demonstrate that the architectural parameters of core-shell comb-like chains have a decisive influence on their final self-assembled conformations. When the ratio of hydrophilic and hydrophobic segments is Ng,g:Ng,r = 3:1, no core-shell structure is observed, and unsegregated chains or clusters are substituted. When Ng,g:Ng,r = 3:2 or 3:3, the core-shell comb-like chains can be assembled into single micelles or structures where several small micelles are strung together by the backbone, which is similar to the pearl-necklace structure formed by the polyelectrolytes in poor solvent. With the increase of backbone length or grafting density, the probability of forming single micelles becomes lower, but the structures of two, three or more small micelles strung together are more observed. Our results indicate how to obtain the desired self-assembled structures of core-shell comb-like chains by regulating the architectural parameters, which could provide insights for the optimization of molecular design in various applications. [ABSTRACT FROM AUTHOR]

Designing the kinetic pathways of assembling macromolecules such as block copolymers and DNA strands is crucial not only for an achievement of thermodynamically equilibrium nanostructures over macroscopic areas, but also for a better understanding of formation process of higher-level superstructures where well-tailored assemblies act as mesoscopic building units. Theoretical analysis and computer simulations provide excellent opportunities to microscopically reveal the kinetics and mechanism of structural evolution as well as the collective behaviors of building units. In this perspective, we summarize our efforts of theoretical and computational modelling to understand the long-range ordering mechanisms and the organization kinetics of assembling macromolecules along designable pathways. First, we present the computational modelling and recent strategies of designable pathways for the achievement of long-range ordering. Then, from the computational views, we give the applications of pathway-designed strategies to explore the ordering mechanism and kinetics in the course of structural evolution, covering the block copolymers and their nanocomposites under zone annealing as well as the hierarchical self-assembly of mesoscopic building units (e.g., patchy micelles and DNA-functionalized nanoparticles). Finally, we outlook future directions in the field of designable pathways for the achievement of long-range ordered nanostructures. This perspective could promote further efforts towards the wide applications of theoretical and computational modelling in the construction of soft hybrid metamaterials. [ABSTRACT FROM AUTHOR]

Dynamic self-consistent field theory (DSCFT) is a fruitful approach for modeling the structural evolution and collective kinetics for a wide variety of multicomponent polymers. However, solving a set of DSCFT equations remains daunting because of high computational demand. Herein, a machine learning method, integrating low-dimensional representations of microstructures and long short-term memory neural networks, is used to accelerate the predictions of structural evolution of multicomponent polymers. It is definitively demonstrated that the neural-network-trained surrogate model has the capability to accurately forecast the structural evolution of homopolymer blends as well as diblock copolymers, without the requirement of "on-the-fly" solution of DSCFT equations. Importantly, the data-driven method can also infer the latent growth laws of phase-separated microstructures of multicomponent polymers through simply using a few of time sequences from their past, without the prior knowledge of the governing dynamics. Our study exemplifies how the machine-learning-accelerated method can be applied to understand and discover the physics of structural evolution in the complex polymer systems. [ABSTRACT FROM AUTHOR]

A trace ratio optimization problem over the Stiefel manifold is investigated from the perspectives of both theory and numerical computations. Necessary conditions in the form of nonlinear eigenvalue problem with eigenvector dependency (NEPv) are established and a numerical method based on the self-consistent field (SCF) iteration with a postprocessing step is designed to solve the NEPv and the method is proved to be always convergent. As an application to multi-view subspace learning, a new framework and its instantiated concrete models are proposed and demonstrated on real world data sets. Numerical results show that the efficiency of the proposed numerical methods and effectiveness of the new orthogonal multi-view subspace learning models. [ABSTRACT FROM AUTHOR]

The phase behavior of ABCBA linear pentablock terpolymers is investigated by using the 3-dimensional self-consistent field theory. In this study, phase diagrams are constructed and used to discuss how the self-assembled morphologies are influenced by the compositions and the Flory–Huggins interaction parameters (χ) among the three components by decreasing symmetric χN value (segregation strength, where N is the total degree of polymerization) from 80 to 30. In the segregation regime of χN from 80 to 50, the microstructures formed by ABCBA linear pentablock terpolymers according to the compositions of the components are very similar. In particular, diverse complex network structures (e.g. diamond, hexagonally perforated lamellae, Fddd, and gyroid) and binary crystalline phases of cylinders and spheres can be observed. This is mainly due to the fact that the two free ends of the A block in the ABCBA linear copolymer allow the macromolecules to relieve packing frustrations. However, in the intermediate system with symmetric χN < 50, the aggregation of each components becomes weaker so that behavior of the pentablock chains is similar to triblock and diblock chains. Accordingly, the diamond and hexagonally perforated lamellae tend to transfer to gyroid and Fddd observed in linear diblock and triblock cases. Moreover, by altering the composition ratio of A/C and the length of the B block, alternatively arranged A/C spheres resemble ionic and metallic crystals (e.g. NaCl, CsCl, Li3Bi, and Nb3Sn) and alternating A/C cylinders with coordination numbers of A/C (equal to 4/4, 6/3, and 4/2) can be still observed by decreasing symmetric χN to 40 in the intermediate regime. [ABSTRACT FROM AUTHOR]

Recently, a few AB‐type multiblock copolymers have been successfully designed to form stable square cylinder phase based on self‐consistent field theory (SCFT) calculations. The previous works only identify the stability region of the square phase but not analyzing its stability, which is closely related to the free‐energy landscape. In this work, we have reexamined the stability of the square phase in B1A1B2A2B3 ${{{\rm B}}_{1}{{\rm A}}_{1}{{\rm B}}_{2}{{\rm A}}_{2}{{\rm B}}_{3}}$ linear pentablock and (B1AB2)​5 ${({{\rm B}}_{1}{\rm A}{{\rm B}}_{2}{)}_{5}}$ star triblock copolymers by drawing the free‐energy landscape with respect to the two dimensions of a rectangular unit cell. Our results demonstrate that the square phase continuously transfers to the rectangular phase as the degree of packing frustration is gradually released. Moreover, the prolate contour lines of the free‐energy landscape indicate the weak stability of the square phase in the B1A1B2A2B3 ${{{\rm B}}_{1}{{\rm A}}_{1}{{\rm B}}_{2}{{\rm A}}_{2}{{\rm B}}_{3}}$ copolymer. In contrast, the stability of the square phase is notably improved in the (B1AB2)​5 ${({{\rm B}}_{1}{\rm A}{{\rm B}}_{2}{)}_{5}}$ copolymer due to its enhanced concentration of bridging configurations. Our work sheds light on the understanding of the stability of the square cylinder phase in block copolymers. Accordingly, we propose some possible strategies for further designing new AB‐type block copolymer systems to obtain more stable square phase. [ABSTRACT FROM AUTHOR]