-
A methodology of quantifying membrane permeability based on returning probability theory and molecular dynamics simulation
Authors:
Yuya Matsubara,
Ryo Okabe,
Ren Masayama,
Nozomi Morishita Watanabe,
Hiroshi Umakoshi,
Kento Kasahara,
Nobuyuki Matubayasi
Abstract:
We propose a theoretical approach to estimate the permeability coefficient of substrates (permeants) for crossing membranes from donor (D) phase to acceptor (A) phase by means of molecular dynamics (MD) simulation. A fundamental aspect of our approach involves reformulating the returning probability (RP) theory, a rigorous bimolecular reaction theory, to describe permeation phenomena. This reformu…
▽ More
We propose a theoretical approach to estimate the permeability coefficient of substrates (permeants) for crossing membranes from donor (D) phase to acceptor (A) phase by means of molecular dynamics (MD) simulation. A fundamental aspect of our approach involves reformulating the returning probability (RP) theory, a rigorous bimolecular reaction theory, to describe permeation phenomena. This reformulation relies on the parallelism between permeation and bimolecular reaction processes. In the present method, the permeability coefficient is represented in terms of the thermodynamic and kinetic quantities for the reactive (R) phase that exists within the inner region of membranes. One can evaluate these quantities using multiple MD trajectories starting from phase R. We apply the RP theory to the permeation of ethanol and methylamine at different concentrations (infinitely dilute and 1 mol% conditions of permeants). Under the 1 mol% condition, the present method yields a larger permeability coefficient for ethanol ($0.12 \pm 0.01 ~\mathrm{cm~s^{-1}}$) than for methylamine ($0.069\pm 0.006~\mathrm{cm~s^{-1}}$), while the values of the permeability coefficient are satisfactorily close to those obtained from the brute-force MD simulations [$0.18\pm 0.03 ~\mathrm{cm~s^{-1}}$ and $0.052 \pm 0.005 ~\mathrm{cm~s^{-1}}$ for ethanol and methylamine, respectively]. Moreover, upon analyzing the thermodynamic and kinetic contributions to the permeability, we clarify that a higher concentration dependency of permeability for ethanol, as compared to methylamine, arises from the sensitive nature of ethanol's free-energy barrier within the inner region of the membrane against ethanol concentration.
△ Less
Submitted 6 June, 2024; v1 submitted 17 April, 2024;
originally announced April 2024.
-
Influence of cholesterol on hydrogen-bond dynamics of water molecules in lipid-bilayer systems at varying temperatures
Authors:
Kokoro Shikata,
Kento Kasahara,
Nozomi Morishita Watanabe,
Hiroshi Umakoshi,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoyl phosphatidylcholine (DPPC) and palmitoyl sphingomyelin (PSM); they are distinct with respect to…
▽ More
Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoyl phosphatidylcholine (DPPC) and palmitoyl sphingomyelin (PSM); they are distinct with respect to the structures of the hydrogen-bond (H-bond) acceptors. Our investigation focuses on the dynamic properties and H-bonds of water molecules in the lipid-membrane systems, with particular emphasis on the influence of Chol at varying temperatures. Notably, in the gel phase at 303 K, the presence of Chol extends the lifetimes of H-bonds of the oxygen atoms acting as H-bond acceptors within DPPC with water molecules by a factor of 1.5 to 2.5. In the liquid-crystalline phase at 323 K, on the other hand, H-bonding dynamics with lipid membranes remain largely unaffected by Chol. This observed shift in H-bonding states serves as a crucial key to unraveling the subtle control mechanisms governing water dynamics in lipid-membrane systems.
△ Less
Submitted 3 July, 2024; v1 submitted 13 March, 2024;
originally announced March 2024.
-
Unveiling interatomic distances influencing the reaction coordinates in alanine dipeptide isomerization: An explainable deep learning approach
Authors:
Kazushi Okada,
Takuma Kikutsuji,
Kei-ichi Okazaki,
Toshifumi Mori,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
The present work shows that the free energy landscape associated with alanine dipeptide isomerization can be effectively represented by specific interatomic distances without explicit reference to dihedral angles. Conventionally, two stable states of alanine dipeptide in vacuum, i.e., C$7_{\mathrm{eq}}$ ($β$-sheet structure) and C$7_{\mathrm{ax}}$ (left handed $α$-helix structure), have been prima…
▽ More
The present work shows that the free energy landscape associated with alanine dipeptide isomerization can be effectively represented by specific interatomic distances without explicit reference to dihedral angles. Conventionally, two stable states of alanine dipeptide in vacuum, i.e., C$7_{\mathrm{eq}}$ ($β$-sheet structure) and C$7_{\mathrm{ax}}$ (left handed $α$-helix structure), have been primarily characterized using the main chain dihedral angles, $\varphi$ (C-N-C$_α$-C) and $ψ$ (N-C$_α$-C-N). However, our recent deep learning combined with "Explainable AI" (XAI) framework has shown that the transition state can be adequately captured by a free energy landscape using $\varphi$ and $θ$ (O-C-N-C$_α$) [T. Kikutsuji, et al. J. Chem. Phys. 156, 154108 (2022)]. In perspective of extending these insights to other collective variables, a more detailed characterization of transition state is required. In this work, we employ the interatomic distances and bond angles as input variables for deep learning, rather than the conventional and more elaborate dihedral angles. Our approach utilizes deep learning to investigate whether changes in the main chain dihedral angle can be expressed in terms of interatomic distances and bond angles. Furthermore, by incorporating XAI into our predictive analysis, we quantified the importance of each input variable and succeeded in clarifying the specific interatomic distance that affects the transition state. The results indicate that constructing a free energy landscape based on using the identified interatomic distance can clearly distinguish between the two stable states and provide a comprehensive explanation for the energy barrier crossing.
△ Less
Submitted 24 April, 2024; v1 submitted 13 February, 2024;
originally announced February 2024.
-
Unraveling the glass-like dynamic heterogeneity in ring polymer melts: From semi-flexible to stiff chain
Authors:
Shota Goto,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Ring polymers are an intriguing class of polymers with unique physical properties, and understanding their behavior is important for developing accurate theoretical models. In this study, we investigate the effect of chain stiffness and monomer density on static and dynamic behaviors of ring polymer melts using molecular dynamics simulations. Our first focus is on the non-Gaussian parameter of cen…
▽ More
Ring polymers are an intriguing class of polymers with unique physical properties, and understanding their behavior is important for developing accurate theoretical models. In this study, we investigate the effect of chain stiffness and monomer density on static and dynamic behaviors of ring polymer melts using molecular dynamics simulations. Our first focus is on the non-Gaussian parameter of center of mass displacement as a measure of dynamic heterogeneity, which is commonly observed in glass-forming liquids. We find that the non-Gaussianity in the displacement distribution increases with the monomer density and stiffness of the polymer chains, suggesting that excluded volume interactions between centers of mass have a stronger effect on the dynamics of ring polymers. We then analyze the relationship between the radius of gyration and monomer density for semi-flexible and stiff ring polymers. Our results indicate that the relationship between the two varies with chain stiffness, which can be attributed to the competition between repulsive forces inside the ring and from adjacent rings. Finally, we study the dynamics of bond-breakage virtually connected between the centers of mass of rings to analyze the exchanges of inter-molecular networks of bonds. Our results demonstrate that the dynamic heterogeneity of bond-breakage is coupled with the non-Gaussianity in ring polymer melts, highlighting the importance of bond-breaking method in determining the inter-molecular dynamics of ring polymer melts. Overall, our study provides insights into the fundamental mechanism of ring polymers and sheds light on the factors that govern their dynamic behavior.
△ Less
Submitted 8 August, 2023; v1 submitted 25 May, 2023;
originally announced May 2023.
-
Revealing the hidden dynamics of confined water in acrylate polymers: Insights from hydrogen-bond lifetime analysis
Authors:
Kokoro Shikata,
Takuma Kikutsuji,
Nobuhiro Yasoshima,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Polymers contain functional groups that participate in hydrogen bond (H-bond) with water molecules, establishing a robust H-bond network that influences bulk properties. This study utilized molecular dynamics (MD) simulations to examine the H-bonding dynamics of water molecules confined within three poly(meth)acrylates: poly(2-methoxyethyl acrylate) (PMEA), poly(2-hydroxyethyl methacrylate) (PHEMA…
▽ More
Polymers contain functional groups that participate in hydrogen bond (H-bond) with water molecules, establishing a robust H-bond network that influences bulk properties. This study utilized molecular dynamics (MD) simulations to examine the H-bonding dynamics of water molecules confined within three poly(meth)acrylates: poly(2-methoxyethyl acrylate) (PMEA), poly(2-hydroxyethyl methacrylate) (PHEMA), and poly(1-methoxymethyl acrylate) (PMC1A). Results showed that H-bonding dynamics significantly slowed as the water content decreased. Additionally, the diffusion of water molecules and its correlation with H-bond breakage were analyzed. Our findings suggest that when the H-bonds between water molecules and the methoxy oxygen of PMEA are disrupted, those water molecules persist in close proximity and do not diffuse on a picosecond timescale. In contrast, the water molecules H-bonded with the hydroxy oxygen of PHEMA and the methoxy oxygen of PMC1A diffuse concomitantly with the breakage of H-bonds. These results provide an in-depth understanding of the impact of polymer functional groups on H-bonding dynamics.
△ Less
Submitted 21 April, 2023; v1 submitted 2 March, 2023;
originally announced March 2023.
-
Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI)
Authors:
Takuma Kikutsuji,
Yusuke Mori,
Kei-ichi Okazaki,
Toshifumi Mori,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is dif…
▽ More
A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.
△ Less
Submitted 1 April, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
Simulating the nematic-isotropic phase transition of liquid crystal model via generalized replica-exchange method
Authors:
Kengo Takemoto,
Yoshiki Ishii,
Hitoshi Washizu,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl (5CB) was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in gREM, allowing the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyz…
▽ More
The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl (5CB) was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in gREM, allowing the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.
△ Less
Submitted 14 December, 2021; v1 submitted 27 September, 2021;
originally announced September 2021.
-
Effects of chain length on Rouse modes and non-Gaussianity in linear and ring polymer melts
Authors:
Shota Goto,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
The dynamics of ring polymer melts are studied via molecular dynamics simulations of the Kremer-Grest bead-spring model. Rouse mode analysis is performed in comparison with linear polymers by changing the chain length. Rouse-like behavior is observed in ring polymers by quantifying the chain length dependence of the Rouse relaxation time, whereas a crossover from Rouse to reptation behavior is obs…
▽ More
The dynamics of ring polymer melts are studied via molecular dynamics simulations of the Kremer-Grest bead-spring model. Rouse mode analysis is performed in comparison with linear polymers by changing the chain length. Rouse-like behavior is observed in ring polymers by quantifying the chain length dependence of the Rouse relaxation time, whereas a crossover from Rouse to reptation behavior is observed in linear polymers. Furthermore, the non-Gaussian parameters of the monomer bead displacement and chain center-of-mass displacement are analyzed. It is found that the non-Gaussianity of ring polymers is remarkably suppressed with slight growth for the center-of-mass dynamics at long chain length, which is in contrast to the growth in linear polymers both for the monomer bead and center-of-mass dynamics.
△ Less
Submitted 16 September, 2021; v1 submitted 25 June, 2021;
originally announced June 2021.
-
Understanding the scaling of boson peak through insensitivity of elastic heterogeneity to bending rigidity in polymer glasses
Authors:
Naoya Tomoshige,
Shota Goto,
Hideyuki Mizuno,
Tatsuya Mori,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Amorphous materials exhibit peculiar mechanical and vibrational properties, including non-affine elastic responses and excess vibrational states, i.e., the so-called boson peak. For polymer glasses, these properties are considered to be affected by the bending rigidity of the constituent polymer chains. In our recent work [Tomoshige, et al., Sci. Rep. 9 19514 (2019)], we have revealed simple relat…
▽ More
Amorphous materials exhibit peculiar mechanical and vibrational properties, including non-affine elastic responses and excess vibrational states, i.e., the so-called boson peak. For polymer glasses, these properties are considered to be affected by the bending rigidity of the constituent polymer chains. In our recent work [Tomoshige, et al., Sci. Rep. 9 19514 (2019)], we have revealed simple relationships between the variations of vibrational properties and the global elastic properties: the response of the boson peak scales only with that of the global shear modulus. This observation suggests that the spatial heterogeneity of the local shear modulus distribution is insensitive to changes in the bending rigidity. Here, we demonstrate the insensitivity of elastic heterogeneity by directly measuring the local shear modulus distribution. We also study transverse sound wave propagation, which is also shown to scale only with the global shear modulus. Through these analyses, we conclude that the bending rigidity does not alter the spatial heterogeneity of the local shear modulus distribution, which yields vibrational and acoustic properties that are controlled solely by the global shear modulus of a polymer glass.
△ Less
Submitted 30 May, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide
Authors:
Yusuke Mori,
Kei-ichi Okazaki,
Toshifumi Mori,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
We propose a cross-entropy minimization method for finding the reaction coordinate from a large number of collective variables in complex molecular systems. This method is an extension of the likelihood maximization approach describing the committor function with a sigmoid. By design, the reaction coordinate as a function of various collective variables is optimized such that the distribution of t…
▽ More
We propose a cross-entropy minimization method for finding the reaction coordinate from a large number of collective variables in complex molecular systems. This method is an extension of the likelihood maximization approach describing the committor function with a sigmoid. By design, the reaction coordinate as a function of various collective variables is optimized such that the distribution of the committor $p_\mathrm{B}^*$ values generated from molecular dynamics simulations can be described in a sigmoidal manner. We also introduce the $L_2$-norm regularization used in the machine learning field to prevent overfitting when the number of considered collective variables is large. The current method is applied to study the isomerization of alanine dipeptide in vacuum, where 45 dihedral angles are used as candidate variables. The regularization parameter is determined by cross-validation using training and test datasets. It is demonstrated that the optimal reaction coordinate involves important dihedral angles, which are consistent with the previously reported results. Furthermore, the points with $p_\mathrm{B}^*\sim 0.5$ clearly indicate a separatrix distinguishing reactant and product states on the potential of mean force using the extracted dihedral angles.
△ Less
Submitted 29 July, 2020; v1 submitted 29 March, 2020;
originally announced March 2020.
-
Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending
Authors:
Naoya Tomoshige,
Hideyuki Mizuno,
Tatsuya Mori,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increas…
▽ More
The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus $G$. In particular, the boson peak frequency is scaled as $ω_\mathrm{BP} \propto \sqrt{G}$. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.
△ Less
Submitted 5 August, 2019; v1 submitted 25 July, 2019;
originally announced July 2019.
-
Diffusion dynamics of supercooled water modeled with the cage-jump motion and hydrogen-bond rearrangement
Authors:
Takuma Kikutsuji,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
The slow dynamics of glass-forming liquids is generally ascribed to the cage-jump motion. In the cage-jump picture, a molecule remains in a cage formed by neighboring molecules, and after a sufficiently long time, it jumps to escape from the original position by cage-breaking. The clarification of the cage-jump motion is therefore linked to unraveling the fundamental element of the slow dynamics.…
▽ More
The slow dynamics of glass-forming liquids is generally ascribed to the cage-jump motion. In the cage-jump picture, a molecule remains in a cage formed by neighboring molecules, and after a sufficiently long time, it jumps to escape from the original position by cage-breaking. The clarification of the cage-jump motion is therefore linked to unraveling the fundamental element of the slow dynamics. Here, we develop a cage-jump model for the dynamics of supercooled water. The caged and jumping states of a water molecule are introduced with respect to the hydrogen-bond (H-bond) rearrangement process, and describe the motion in supercooled states. It is then demonstrated from the molecular dynamics simulation of the TIP4P/2005 model that the characteristic length and time scales of cage-jump motions provide a good description of the self-diffusion constant that is determined in turn from the long-time behavior of the mean square displacement. Our cage-jump model thus enables to connect between H-bond dynamics and molecular diffusivity.
△ Less
Submitted 30 April, 2019; v1 submitted 14 March, 2019;
originally announced March 2019.
-
Bridging the gap between molecular dynamics and hydrodynamics in nanoscale Brownian motions
Authors:
Keisuke Mizuta,
Yoshiki Ishii,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Through molecular dynamics simulations, we examined hydrodynamic behavior of the Brownian motion of fullerene particles based on molecular interactions. The solvation free energy and the velocity autocorrelation function (VACF) were calculated by using the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) potentials for the solute-solvent and solvent-solvent interactions and by changing the siz…
▽ More
Through molecular dynamics simulations, we examined hydrodynamic behavior of the Brownian motion of fullerene particles based on molecular interactions. The solvation free energy and the velocity autocorrelation function (VACF) were calculated by using the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) potentials for the solute-solvent and solvent-solvent interactions and by changing the size of the fullerene particles. We also measured the diffusion constant of the fullerene particles and the shear viscosity of the host fluid, and then the hydrodynamic radius $a_\mathrm{HD}$ was quantified from the Stokes-Einstein relation. The $a_\mathrm{HD}$ value exceeds that of the gyration radius of the fullerene when the solvation free energy exhibits largely negative values using the LJ potential. In contrast, $a_\mathrm{HD}$ becomes comparable to the size of bare fullerene, when the solvation free energy is positive using the WCA potential. Furthermore, the VACF of the fullerene particles is directly comparable with the analytical expressions utilizing the Navier-Stokes equations both in incompressible and compressible forms. Hydrodynamic long-time tail $t^{-3/2}$ is demonstrated for timescales longer than the kinematic time of the momentum diffusion over the particles' size. However, the VACF in shorter timescales deviates from the hydrodynamic description, particularly for smaller fullerene particles and for the LJ potential. This occurs even though the compressible effect is considered when characterizing the decay of VACF around the sound propagation time scale over the particles' size. These results indicate that the nanoscale Brownian motion is influenced by the solvation structure around the solute particles originating from the molecular interaction.
△ Less
Submitted 23 April, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
How do hydrogen bonds break in supercooled water?: Detecting pathways not going through the saddle point of two-dimensional potential of mean force
Authors:
Takuma Kikutsuji,
Kang Kim,
Nobuyuki Matubayasi
Abstract:
Supercooled water exhibits remarkably slow dynamics similar to the behavior observed for various glass-forming liquids. The local order of tetrahedral structures due to hydrogen-bonds (H-bonds) increases with decreasing temperature. Thus, it is important to clarify the temperature dependence of the H-bond breakage process. This was investigated here using molecular dynamics simulations of TIP4P su…
▽ More
Supercooled water exhibits remarkably slow dynamics similar to the behavior observed for various glass-forming liquids. The local order of tetrahedral structures due to hydrogen-bonds (H-bonds) increases with decreasing temperature. Thus, it is important to clarify the temperature dependence of the H-bond breakage process. This was investigated here using molecular dynamics simulations of TIP4P supercooled water. The two-dimensional (2D) potential of mean force (PMF) is presented using combinations of intermolecular distance and angle between two water molecules. The saddle point of the 2D PMF suggests the presence of the transition state that distinguishes between H-bond and non H-bond states. However, we observed pathways not going through this saddle point particularly at supercooled states, which are due to translational, rather than rotational motions of the molecules. We quantified the characteristic time scales of rotational and translational H-bond breakages. The time scale of translational H-bond breakage shows a non-Arrhenius temperature dependence comparable to that of the H-bond lifetime. This time scale is relevant for the temperature dependence of the transmission coefficient based on the transition state theory. The translational H-bond breakage is also related to cage-jumps observed in glass-forming liquids, which mostly involve spatially correlated molecules. Our findings warrant further exploration of an appropriate free-energy surface or reaction coordinates beyond the geometrical variables of the water dimer to describe a possible saddle point related to collective jump motions.
△ Less
Submitted 16 June, 2018; v1 submitted 4 March, 2018;
originally announced March 2018.