Fast, Scalable, and Interactive Software for Landau-de Gennes Numerical Modeling of Nematic Topological Defects

Numerical modeling of nematic liquid crystals using the tensorial Landau-de Gennes (LdG) theory provides detailed insights into the structure and energetics of the enormous variety of possible topological defect configurations that may arise when the liquid crystal is in contact with colloidal inclusions or structured boundaries. However, these methods can be computationally expensive, making it challenging to predict (meta)stable configurations involving several colloidal particles, and they are often restricted to system sizes well below the experimental scale. Here we present an open-source software package that exploits the embarrassingly parallel structure of the lattice discretization of the LdG approach. Our implementation, combining CUDA/C++ and OpenMPI, allows users to accelerate simulations using both CPU and GPU resources in either single- or multiple-core configurations. We make use of an efficient minimization algorithm, the Fast Inertial Relaxation Engine (FIRE) method, that is well-suited to large-scale parallelization, requiring little additional memory or computational cost while offering performance competitive with other commonly used methods. In multi-core operation we are able to scale simulations up to supra-micron length scales of experimental relevance, and in single-core operation the simulation package includes a user-friendly GUI environment for rapid prototyping of interfacial features and the multifarious defect states they can promote. To demonstrate this software package, we examine in detail the competition between curvilinear disclinations and point-like hedgehog defects as size scale, material properties, and geometric features are varied. We also study the effects of an interface patterned with an array of topological point-defects.Comment: 16 pages, 6 figures, 1 youtube link. The full catastroph

Comment on "Quantum waveguide array generator for performing Fourier transforms: Alternate route to quantum computing"

A Comment on the paper "Quantum waveguide array generator for performing Fourier transforms: Alternate route to quantum computing" by R. Akis and D.K. Ferry, Appl. Phys. Lett. 79, 2823 (2001). The authors reply in Appl. Phys. Lett. 80, 2420 (2002).Comment: 1 pag

B-modes and the Nature of Inflation

Observations of the cosmic microwave background do not yet determine whether inflation was driven by a slowly-rolling scalar field or involved another physical mechanism. In this paper we discuss the prospects of using the power spectra of scalar and tensor modes to probe the nature of inflation. We focus on the leading modification to the slow-roll dynamics, which entails a sound speed $c_s$ for the scalar fluctuations. We derive analytically a lower bound on $c_s$ in terms of a given tensor-to-scalar ratio $r$, taking into account the difference in the freeze-out times between the scalar and tensor modes. We find that any detection of primordial B-modes with $r > 0.01$ implies a lower bound on $c_s$ that is stronger than the bound derived from the absence of non-Gaussianity in the Planck data. For $r \gtrsim 0.1$, the bound would be tantalizingly close to a critical value for the sound speed, $(c_s)_\star = 0.47$ (corresponding to $(f_{\rm NL}^{\rm equil})_\star = -0.93$), which we show serves as a threshold for non-trivial dynamics beyond slow-roll. We also discuss how an order-one level of equilateral non-Gaussianity is a natural observational target for other extensions of the canonical paradigm.Comment: 25+7 pages, 9 figures. Published versio