Search

Aiming to provide a very accurate, efficient, and robust quartic equation solver for physical applications, we have proposed an algorithm that builds on the previous works of P. Strobach and S. L. Shmakov. It is based on the decomposition of the quartic polynomial into two quadratics, whose coefficients are first accurately estimated by handling carefully numerical errors and afterward refined through the use of the Newton-Raphson method. Our algorithm is very accurate in comparison with other state-of-the-art solvers that can be found in the literature, but (most importantly) it turns out to be very efficient according to our timing tests. A crucial issue for us is the robustness of the algorithm, i.e., its ability to cope with the detrimental effect of round-off errors, no matter what set of quartic coefficients is provided in a practical application. In this respect, we extensively tested our algorithm in comparison to other quartic equation solvers both by considering specific extreme cases and by carrying out a statistical analysis over a very large set of quartics. Our algorithm has also been heavily tested in a physical application, i.e., simulations of hard cylinders, where it proved its absolute reliability as well as its efficiency.

Double-stranded DNA (dsDNA) fragments exhibit noncovalent attractive interactions between their tips. It is still unclear how DNA liquid crystal self-assembly is affected by such blunt-end attractions. It is demonstrated that stiff dsDNA fragments with moderate aspect ratio can specifically self-assemble in concentrated aqueous solutions into different types of smectic mesophases on the basis of selectively screening of blunt-end DNA stacking interactions. To this end, this type of attractions are engineered at the molecular level by constructing DNA duplexes where the attractions between one or both ends are screened by short hairpin caps. All-DNA bilayer and monolayer smectic-A type of phases, as well as a columnar phase, can be stabilized by controlling attractions strength. The results imply that the so far elusive smectic-A in DNA rod-like liquid crystals is a thermodynamically stable phase. The existence of the bilayer smectic phase is confirmed by Monte-Carlo simulations of hard cylinders decorated with one attractive terminal site. This work demonstrates that DNA blunt-ends behave as well-defined monovalent attractive patches whose strength and position can be potentially precisely tuned, highlighting unique opportunities concerning the stabilization of nonconventional DNA-based lyotropic liquid crystal phases assembled by all-DNA patchy particles with arbitrary geometry and composition.

Liquid-crystalline phases in all-DNA systems have been extensively studied in the past and although nematic, cholesteric and columnar mesophases have been observed, the smectic phase remained elusive. Recently, it has been found evidence of a smectic-A ordering in an all-DNA system, where the constituent particles, which are gapped DNA duplexes, resemble chain-sticks. It has been argued that in the smectic-A phase these DNA chain-sticks should be folded as a means to suppress aggregate polydispersity and excluded volume. Nevertheless, if initial crystalline configurations are prepared in silico with gapped DNA duplexes either fully unfolded or fully folded by carrying out computer simulations one can end up with two different phases having at the same concentration and temperature the majority of gapped DNA duplexes either folded or unfolded. This result suggests that these two phases have a small free energy difference, since no transition is observed from one to the other within the simulation time span. In the present manuscript, we assess which of these two phases is thermodynamically stable through a suitable protocol based on thermodynamic integration. Our method is rather general and it can be used to discriminate stable states from metastable ones of comparable free energy.

The hard cylinder model decorated with attractive patches proved to be very useful recently in studying physical properties of several colloidal systems. Phase diagram, elastic constants and cholesteric properties obtained from computer simulations based on a simple hard cylinder model have been all successfully and quantitatively compared to experimental results. Key to these simulations is an efficient algorithm to check the overlap between hard cylinders. Here, we propose two algorithms to check the hard cylinder overlap and we assess their efficiency through a comparison with an existing method available in the literature and with the well-established algorithm for simulating hard spherocylinders. In addition, we discuss a couple of optimizations for performing computer simulations of patchy anisotropic particles and we estimate the speed-up which they can provide in the case of patchy hard cylinders.

Multicenter transition metal complexes are the key moieties of many processes in chemistry, biochemistry, and materials science such as in the active sites of enzymes, molecular catalysts, and biological electron carriers. Their electronic structure, often characterized by high-spin-polarized metal sites, is a challenge for theoretical chemists because of their high degree of dynamical and static correlation. Static correlation is necessary both for the appropriate description of the metal-ligand bonding and for a correct description of the multideterminant character arising from the magnetic interactions between spin centers. Density functional theory (DFT) is usually applied using a single-determinant broken-symmetry state that is lacking the correct spin symmetry when the ground state has total low-spin character. To alleviate this drawback, we use the extended broken-symmetry (EBS) approach to derive approximate ground-state energies and, for the first time, forces for the correctly symmetric ground state of an arbitrary number of spin centers within the framework of the Heisenberg-Dirac-van Vleck Hamiltonian. Remarkably, the proposed procedure supplies relaxed geometries that are fully consistent with the calculated J-coupling constants. We apply the method to investigate the relaxed geometrical structure of the low-spin ground state of iron-sulfur clusters with two, three, and four iron centers. We observed significant differences in both geometrical parameters and coupling constant J between the symmetrized ground state, the high-spin, and the broken-symmetry optimized structures. These changes are often comparable with the differences observed by using different functionals, and the use of EBS always improves the description of the studied systems. It will be therefore important to include it in any DFT attempt to quantitatively describe multicenter transition metal complexes in the future.