Showing 1–32 of 32 results for author: Mori, Y

Detections of He-3 in Ni-based binary metal nanocomposites with Cu in zirconia exposed to hydrogen gas at elevated temperatures

Authors: Tomoya Yamauchi, Yutaka Mori, Shuto Higashi, Hayato Seiichi, Masahiko Hasegawa, Akito Takahashi, Akira Taniike, Masato Kanasaki

Abstract: The present study aims to detect helium-3 in nickel-based metal nano-composites doped with zirconia, which exhibited anomalous heat generation when exposed to hydrogen gas at approximately 450�C. Two complementary analytical techniques were employed: Nuclear Reaction Analysis (NRA) utilizing 1.4 MeV deuteron beams from a tandem accelerator, and Thermal Desorption Spectrometry (TDS) using a quadrup… ▽ More The present study aims to detect helium-3 in nickel-based metal nano-composites doped with zirconia, which exhibited anomalous heat generation when exposed to hydrogen gas at approximately 450�C. Two complementary analytical techniques were employed: Nuclear Reaction Analysis (NRA) utilizing 1.4 MeV deuteron beams from a tandem accelerator, and Thermal Desorption Spectrometry (TDS) using a quadrupole mass spectrometer. Both methods successfully detected helium-3 in the samples. Given the extreme rarity of this isotope, its presence strongly suggests the occurrence of nuclear reactions within the nickel-containing materials. These findings lend support to the 4H/TSC (4 Hydrogen/Tetrahedral Symmetric Condensate) model, which uniquely predicts helium-3 as one of the primary reaction products. △ Less

Submitted 9 September, 2024; originally announced September 2024.

Comments: 18 pages, 7 figures. This is the version of the article before peer review, as submitted by an author to Japanese Journal of Applied Physics

Beam Stacking Experiment at a Fixed Field Alternating Gradient Accelerator

Authors: T. Uesugi, Y. Ishi, Y. Kuriyama, Y. Mori, C. Jolly, D. J. Kelliher, J. -B. Lagrange, A. P. Letchford, S. Machida, D. W. Poshuma de Boer, C. T. Rogers, E. Yamakawa, M. Topp-Mugglestone

Abstract: A key challenge in particle accelerators is to achieve high peak intensity. Space charge is particularly strong at lower energy such as during injection and typically limits achievable peak intensity. The beam stacking technique can overcome this limitation by accumulating a beam at high energy where space charge is weaker. In beam stacking, a bunch of particles is injected and accelerated to high… ▽ More A key challenge in particle accelerators is to achieve high peak intensity. Space charge is particularly strong at lower energy such as during injection and typically limits achievable peak intensity. The beam stacking technique can overcome this limitation by accumulating a beam at high energy where space charge is weaker. In beam stacking, a bunch of particles is injected and accelerated to high energy. This bunch continues to circulate, while a second and subsequent bunches are accelerated to merge into the first. It also allows the user cycle and acceleration cycles to be separated which is often valuable. Beam stacking is not possible in a time varying magnetic field, but a fixed field machine such as an Fixed Field Alternating Gradient Accelerator (FFA) does not sweep the magnetic field. In this paper, we describe experimental demonstration of beam stacking of two beams at KURNS FFA in Kyoto University. The momentum spread and intensity of the beam was analysed by study of the Schottky signal, demonstrating stacking with only a slight increase of momentum spread of the combined beams. The intensity of the first beam was, however, significantly reduced. RF knock-out is the suspected source of the beam loss. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 14 pages, 16 figures

TiO2 doping effect on reflective coating mechanical loss for gravitational wave detection at low temperature

Authors: Yukino Mori, Yota Nakayama, Kazuhiro Yamamoto, Takafumi Ushiba, Dani`ele Forest, Christophe Michel, Laurent Pinard, Julien Teillo, Gianpietro Cagnoli

Abstract: We measured the mechanical loss of a dielectric multilayer reflective coating (ion-beam-sputtered SiO2 and Ta2O5) with and without TiO2 on sapphire disks between 6 and 77 K. The measured loss angle exhibited a temperature dependence, and the local maximum was found at approximately 20 K. This maximum was 7.0*10^(-4) (with TiO2) and 7.7*10^(-4) (without TiO2), although the previous measurement for… ▽ More We measured the mechanical loss of a dielectric multilayer reflective coating (ion-beam-sputtered SiO2 and Ta2O5) with and without TiO2 on sapphire disks between 6 and 77 K. The measured loss angle exhibited a temperature dependence, and the local maximum was found at approximately 20 K. This maximum was 7.0*10^(-4) (with TiO2) and 7.7*10^(-4) (without TiO2), although the previous measurement for the coating on sapphire disks showed almost no temperature dependence (Phys. Rev. D 74 022002 (2006)). We evaluated the coating thermal noise in KAGRA and discussed future investigation strategies. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Finite population effects on optimal communication for social foragers

Authors: Hyunjoong Kim, Yoichiro Mori, Joshua B Plotkin

Abstract: Foraging is crucial for animals to survive. Many species forage in groups, as individuals communicate to share information about the location of available resources. For example, eusocial foragers, such as honey bees and many ants, recruit members from their central hive or nest to a known foraging site. However, the optimal level of communication and recruitment depends on the overall group size,… ▽ More Foraging is crucial for animals to survive. Many species forage in groups, as individuals communicate to share information about the location of available resources. For example, eusocial foragers, such as honey bees and many ants, recruit members from their central hive or nest to a known foraging site. However, the optimal level of communication and recruitment depends on the overall group size, the distribution of available resources, and the extent of interference between multiple individuals attempting to forage from a site. In this paper, we develop a discrete-time Markov chain model of eusocial foragers, who communicate information with a certain probability. We compare the stochastic model and its corresponding infinite-population limit. We find that foraging efficiency tapers off when recruitment probability is too high -- a phenomenon that does not occur in the infinite-population model, even though it occurs for any finite population size. The marginal inefficiency at high recruitment probability increases as the population increases, similar to a boundary layer. In particular, we prove there is a significant gap between the foraging efficiency of finite and infinite population models in the extreme case of complete communication. We also analyze this phenomenon by approximating the stationary distribution of foragers over sites in terms of mean escape times from multiple quasi-steady states. We conclude that for any finite group of foragers, an individual who has found a resource should only sometimes recruit others to the same resource. We discuss the relationship between our analysis and multi-agent multi-arm bandit problems. △ Less

Submitted 1 August, 2023; originally announced August 2023.

A new scheme of negative muon beam generation with MuCF

Authors: Yoshiharu Mori

Abstract: A new scheme of low energy negative muon source with the muon catalyzed fusion (MuCF) is described. In the MuCF reaction process, muonic helium ions ($μ$He$^+$) are created. By re-accelerating and stripping $μ$He$^+$ ions, a low emittance negative muon beam is generated. A new scheme of low energy negative muon source with the muon catalyzed fusion (MuCF) is described. In the MuCF reaction process, muonic helium ions ($μ$He$^+$) are created. By re-accelerating and stripping $μ$He$^+$ ions, a low emittance negative muon beam is generated. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Input optics systems of the KAGRA detector during O3GK

Authors: T. Akutsu, M. Ando, K. Arai, Y. Arai, S. Araki, A. Araya, N. Aritomi, H. Asada, Y. Aso, S. Bae, Y. Bae, L. Baiotti, R. Bajpai, M. A. Barton, K. Cannon, Z. Cao, E. Capocasa, M. Chan, C. Chen, K. Chen, Y. Chen, C-I. Chiang, H. Chu, Y-K. Chu, S. Eguchi , et al. (228 additional authors not shown)

Abstract: KAGRA, the underground and cryogenic gravitational-wave detector, was operated for its solo observation from February 25th to March 10th, 2020, and its first joint observation with the GEO 600 detector from April 7th -- 21st, 2020 (O3GK). This study presents an overview of the input optics systems of the KAGRA detector, which consist of various optical systems, such as a laser source, its intensit… ▽ More KAGRA, the underground and cryogenic gravitational-wave detector, was operated for its solo observation from February 25th to March 10th, 2020, and its first joint observation with the GEO 600 detector from April 7th -- 21st, 2020 (O3GK). This study presents an overview of the input optics systems of the KAGRA detector, which consist of various optical systems, such as a laser source, its intensity and frequency stabilization systems, modulators, a Faraday isolator, mode-matching telescopes, and a high-power beam dump. These optics were successfully delivered to the KAGRA interferometer and operated stably during the observations. The laser frequency noise was observed to limit the detector sensitivity above a few kHz, whereas the laser intensity did not significantly limit the detector sensitivity. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Well-posedness and applications of classical elastohydrodynamics for a swimming filament

Authors: Yoichiro Mori, Laurel Ohm

Abstract: We consider a classical elastohydrodynamic model of an inextensible filament undergoing planar motion in $\mathbb{R}^3$. The hydrodynamics are described by resistive force theory, and the fiber elasticity is governed by Euler-Bernoulli beam theory. Our aim is twofold: (1) Serve as a starting point for developing the mathematical analysis of filament elastohydrodynamics, particularly the analytical… ▽ More We consider a classical elastohydrodynamic model of an inextensible filament undergoing planar motion in $\mathbb{R}^3$. The hydrodynamics are described by resistive force theory, and the fiber elasticity is governed by Euler-Bernoulli beam theory. Our aim is twofold: (1) Serve as a starting point for developing the mathematical analysis of filament elastohydrodynamics, particularly the analytical treatment of an inextensibility constraint, and (2) As an application, prove conditions on internal fiber forcing that allow a free-ended filament to swim. Our analysis of fiber swimming speed is supplemented with a numerical optimization of the internal fiber forcing, as well as a novel numerical method for simulating an inextensible swimmer. △ Less

Submitted 12 August, 2022; v1 submitted 5 August, 2022; originally announced August 2022.

Comments: 38 pages, many figures toward the end. Some typos fixed

Determination of the electron trap level in Fe-doped GaN by phonon-assisted conduction phenomenon

Authors: Hiroki Fukuda, Akira Nagakubo, Shigeyoshi Usami, Masashi Ikeda, Masayuki Imanishi, Masashi Yoshimura, Yusuke Mori, Kanta Adachi, Hirotsugu Ogi

Abstract: We acoustically measured the energy level for thermally activated conduction (TAC) in high-resistivity Fe-doped GaN using the non-contacting antenna-transmission acoustic-resonance method. The acoustic attenuation takes a maximum at a specific temperature, where the TAC is accelerated with the help of phonon energy. The Debye type relaxation is thus observed for acoustic attenuation, and its activ… ▽ More We acoustically measured the energy level for thermally activated conduction (TAC) in high-resistivity Fe-doped GaN using the non-contacting antenna-transmission acoustic-resonance method. The acoustic attenuation takes a maximum at a specific temperature, where the TAC is accelerated with the help of phonon energy. The Debye type relaxation is thus observed for acoustic attenuation, and its activation energy (0.54$\pm$0.04 eV) was determined with attenuation measurements at various frequencies and temperatures. This value agrees with the E3 level in GaN, indicating that thermally associated conduction originates from the E3 trap level. We also measured the five independent elastic constants at high temperatures. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Plasmon-molecule remote coupling via column-structured silica layer for enhancing biophotonic analysis

Authors: Takeo Minamikawa, Reiko Sakaguchi, Yoshinori Harada, Hideharu Hase, Yasuo Mori, Tetsuro Takamatsu, Yu Yamasaki, Yukihiro Morimoto, Masahiro Kawasaki, Mitsuo Kawasaki

Abstract: We demonstrated remote plasmonic enhancement (RPE) by a dense random array of Ag nanoislands (AgNIs) that were partially gold-alloyed and attached with column-structured silica (CSS) overlayer of more than 100 nm in thickness. The physical and chemical protection of the CSS layer could lead to reducing the mutual impact between analyte molecules and metal nanostructures. RPE plate was realized jus… ▽ More We demonstrated remote plasmonic enhancement (RPE) by a dense random array of Ag nanoislands (AgNIs) that were partially gold-alloyed and attached with column-structured silica (CSS) overlayer of more than 100 nm in thickness. The physical and chemical protection of the CSS layer could lead to reducing the mutual impact between analyte molecules and metal nanostructures. RPE plate was realized just by sputtering and chemical immersion processes, resulting in high productivity. We found a significant enhancement on the order of 10$^7$-fold for Raman scattering and 10$^2$-fold for fluorescence by RPE even without the proximity of metal nanostructures and analyte molecules. We confirmed the feasibility of RPE for biophotonic analysis. RPE worked for dye molecules in cells cultured on the CSS layer, enabling the enhanced fluorescence biosensing of intracellular signaling dynamics in HeLa cells. RPE also worked for biological tissues, enhancing Raman histological imaging of esophagus tissues with esophageal adventitia of a Wistar rat attached atop the CSS layer. We also investigated the wavelength dependency of RPE on the on- or off-resonant with the dye molecular transition dipoles with various molecular concentrations. The results suggested that the RPE occurred by remote resonant coupling between the localized surface plasmon of AgNIs and the molecular transition dipole of the analyte via the CSS structure. The RPE plate affords practical advantages for potential biophotonic analyses such as high productivity and biocompatibility. We thus anticipate that RPE will advance to versatile analytical tools in chemistry, biology, and medicine. △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: 34 pages, 11 figures, 1 table

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.

Comments: 9 pages, 5 figures, 2 tables, to appear in J. Chem. Phys

Journal ref: J. Chem. Phys. 156, 154108 (2022)

Evolution of prosocial behavior in multilayer populations

Authors: Qi Su, Alex McAvoy, Yoichiro Mori, Joshua B. Plotkin

Abstract: Human societies include diverse social relationships. Friends, family, business colleagues, and online contacts can all contribute to one's social life. Individuals may behave differently in different domains, but success in one domain may engender success in another. Here, we study this problem using multilayer networks to model multiple domains of social interactions, in which individuals experi… ▽ More Human societies include diverse social relationships. Friends, family, business colleagues, and online contacts can all contribute to one's social life. Individuals may behave differently in different domains, but success in one domain may engender success in another. Here, we study this problem using multilayer networks to model multiple domains of social interactions, in which individuals experience different environments and may express different behaviors. We provide a mathematical analysis and find that coupling between layers tends to promote prosocial behavior. Even if prosociality is disfavored in each layer alone, multilayer coupling can promote its proliferation in all layers simultaneously. We apply this analysis to six real-world multilayer networks, ranging from the socio-emotional and professional relationships in a Zambian community, to the online and offline relationships within an academic University. We discuss the implications of our results, which suggest that small modifications to interactions in one domain may catalyze prosociality in a different domain. △ Less

Submitted 25 October, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 68 pages; final version

Journal ref: Nature Human Behaviour (2022)

Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid

Authors: Yoichiro Mori, Laurel Ohm

Abstract: We consider the accuracy of slender body theory in approximating the force exerted by a thin fiber on the surrounding viscous fluid when the fiber velocity is prescribed. We term this the slender body inverse problem, as it is known that slender body theory converges to a well-posed PDE solution when the force is prescribed and the fiber velocity is unknown. From a PDE perspective, the slender bod… ▽ More We consider the accuracy of slender body theory in approximating the force exerted by a thin fiber on the surrounding viscous fluid when the fiber velocity is prescribed. We term this the slender body inverse problem, as it is known that slender body theory converges to a well-posed PDE solution when the force is prescribed and the fiber velocity is unknown. From a PDE perspective, the slender body inverse problem is simply the Dirichlet problem for the Stokes equations, but from an approximation perspective, nonlocal slender body theory exhibits high wavenumber instabilities which complicate analysis. Here we consider two methods for regularizing the slender body approximation: spectral truncation and the $δ$-regularization of Tornberg and Shelley (2004). For a straight, periodic fiber with constant radius $ε>0$, we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We show that the spectrum of the original slender body approximation agrees closely with the PDE solution at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff $\sim1/ε$. For both the truncated and $δ$-regularized approximations, we obtain similar convergence results to the PDE solution as $ε\to0$: a fiber velocity with $H^1$ regularity gives $O(ε)$ convergence, while a fiber velocity with at least $H^2$ regularity yields $O(ε^2)$ convergence. Moreover, we determine the dependence of the $δ$-regularized error estimate on the regularization parameter $δ$. △ Less

Submitted 15 August, 2020; originally announced August 2020.

Spin-projection for quantum computation: A low-depth approach to strong correlation

Authors: Takashi Tsuchimochi, Yuto Mori, Seiichiro L. Ten-no

Abstract: Although spin is a core property in fermionic systems, its symmetry can be easily violated in a variational simulation, especially when strong correlation plays a vital role therein. In this study, we will demonstrate that the broken spin-symmetry can be restored exactly in a quantum computer, with little overhead in circuits, while delivering additional strong correlation energy with the desired… ▽ More Although spin is a core property in fermionic systems, its symmetry can be easily violated in a variational simulation, especially when strong correlation plays a vital role therein. In this study, we will demonstrate that the broken spin-symmetry can be restored exactly in a quantum computer, with little overhead in circuits, while delivering additional strong correlation energy with the desired spin quantum number. The proposed scheme permits drastic reduction of a potentially large number of measurements required to ensure spin-symmetry by employing a superposition of only a few rotated quantum states. Our implementation is universal, simple, and, most importantly, straightforwardly applicable to any ansatz proposed to date. △ Less

Submitted 20 December, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: 9 pages, 7 figures, 2 tables

Journal ref: Phys. Rev. Research 2, 043142 (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.

Comments: 8 pages, 5 figures, 1 table for main text, 5 pages for supplementary materials, to appear in J. Chem. Phys. for "Classical Molecular Dynamics (MD) Simulations: Codes, Algorithms, Force Fields, and Applications"

Journal ref: J. Chem. Phys. 153, 054115 (2020)

Thermonuclear Fusion Triggered by Collapsing Standing Whistler Waves in Magnetized Overdense Plasmas

Authors: Takayoshi Sano, Shinsuke Fujioka, Yoshitaka Mori, Kunioki Mima, Yasuhiko Sentoku

Abstract: Thermal fusion plasmas initiated by standing whistler waves are investigated numerically by two- and one-dimensional Particle-in-Cell simulations. When a standing whistler wave collapses due to the wave breaking of ion plasma waves, the energy of the electromagnetic waves transfers directly to the ion kinetic energy. Here we find that the ion heating by the standing whistler wave is operational ev… ▽ More Thermal fusion plasmas initiated by standing whistler waves are investigated numerically by two- and one-dimensional Particle-in-Cell simulations. When a standing whistler wave collapses due to the wave breaking of ion plasma waves, the energy of the electromagnetic waves transfers directly to the ion kinetic energy. Here we find that the ion heating by the standing whistler wave is operational even in multi-dimensional simulations of multi-ion species targets, such as deuterium-tritium (DT) ices and solid ammonia borane (H$_6$BN). The energy conversion efficiency to ions becomes as high as 15% of the injected laser energy, which depends significantly on the target thickness and laser pulse duration. The ion temperature could reach a few tens of keV or much higher if appropriate laser-plasma conditions are selected. DT fusion plasmas generated by this method must be useful as efficient neutron sources. Our numerical simulations suggest that the neutron generation efficiency exceeds 10$^9$ n/J per steradian, which is beyond the current achievements of the state-of-the-art laser experiments. The standing whistler wave heating would expand the experimental possibility for an alternative ignition design of magnetically confined laser fusion, and also for more difficult fusion reactions including the aneutronic proton-boron reaction. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 16 pages, 11 figures, accepted for publication in PRE

Journal ref: Phys. Rev. E 101, 013206 (2020)

An error bound for the slender body approximation of a thin, rigid fiber sedimenting in Stokes flow

Authors: Yoichiro Mori, Laurel Ohm

Abstract: We investigate the motion of a thin rigid body in Stokes flow and the corresponding slender body approximation used to model sedimenting fibers. In particular, we derive a rigorous error bound comparing the rigid slender body approximation to the classical PDE for rigid motion in the case of a closed loop with constant radius. Our main tool is the slender body PDE framework established by the auth… ▽ More We investigate the motion of a thin rigid body in Stokes flow and the corresponding slender body approximation used to model sedimenting fibers. In particular, we derive a rigorous error bound comparing the rigid slender body approximation to the classical PDE for rigid motion in the case of a closed loop with constant radius. Our main tool is the slender body PDE framework established by the authors and D. Spirn in [18,19], which we adapt to the rigid setting. △ Less

Submitted 15 March, 2020; v1 submitted 16 April, 2019; originally announced April 2019.

Theoretical justification and error analysis for slender body theory with free ends

Authors: Yoichiro Mori, Laurel Ohm, Daniel Spirn

Abstract: Slender body theory is a commonly used approximation in computational models of thin fibers in viscous fluids, especially in simulating the motion of cilia or flagella in swimming microorganisms. In [23], we developed a PDE framework for analyzing the error introduced by the slender body approximation for closed-loop fibers with constant radius $ε$, and showed that the difference between our close… ▽ More Slender body theory is a commonly used approximation in computational models of thin fibers in viscous fluids, especially in simulating the motion of cilia or flagella in swimming microorganisms. In [23], we developed a PDE framework for analyzing the error introduced by the slender body approximation for closed-loop fibers with constant radius $ε$, and showed that the difference between our closed-loop PDE solution and the slender body approximation is bounded by an expression proportional to $ε|\logε|$. Here we extend the slender body PDE framework to the free endpoint setting, which is more physically relevant from a modeling standpoint but more technically demanding than the closed loop analysis. The main new difficulties arising in the free endpoint setting are defining the endpoint geometry, identifying the extent of the 1D slender body force density, and determining how the well-posedness constants depend on the non-constant fiber radius. Given a slender fiber satisfying certain geometric constraints at the filament endpoints and a one-dimensional force density satisfying an endpoint decay condition, we show a bound for the difference between the solution to the slender body PDE and the slender body approximation in the free endpoint setting. The bound is a sum of the same $ε|\logε|$ term appearing in the closed loop setting and an endpoint term proportional to $ε$, where $ε$ is now the maximum fiber radius. △ Less

Submitted 31 January, 2019; originally announced January 2019.

COMET Phase-I Technical Design Report

Authors: The COMET Collaboration, R. Abramishvili, G. Adamov, R. R. Akhmetshin, A. Allin, J. C. Angélique, V. Anishchik, M. Aoki, D. Aznabayev, I. Bagaturia, G. Ban, Y. Ban, D. Bauer, D. Baygarashev, A. E. Bondar, C. Cârloganu, B. Carniol, T. T. Chau, J. K. Chen, S. J. Chen, Y. E. Cheung, W. da Silva, P. D. Dauncey, C. Densham, G. Devidze , et al. (170 additional authors not shown)

Abstract: The Technical Design for the COMET Phase-I experiment is presented in this paper. COMET is an experiment at J-PARC, Japan, which will search for neutrinoless conversion of muons into electrons in the field of an aluminium nucleus ($μ-e$ conversion, $μ^- N \to e^- N$); a lepton flavor violating process. The experimental sensitivity goal for this process in the Phase-I experiment is… ▽ More The Technical Design for the COMET Phase-I experiment is presented in this paper. COMET is an experiment at J-PARC, Japan, which will search for neutrinoless conversion of muons into electrons in the field of an aluminium nucleus ($μ-e$ conversion, $μ^- N \to e^- N$); a lepton flavor violating process. The experimental sensitivity goal for this process in the Phase-I experiment is $3.1\times10^{-15}$, or 90 % upper limit of branching ratio of $7\times 10^{-15}$, which is a factor of 100 improvement over the existing limit. The expected number of background events is 0.032. To achieve the target sensitivity and background level, the 3.2 kW 8 GeV proton beam from J-PARC will be used. Two types of detectors, CyDet and StrECAL, will be used for detecting the \mue conversion events, and for measuring the beam-related background events in view of the Phase-II experiment, respectively. Results from simulation on signal and background estimations are also described. △ Less

Submitted 19 May, 2020; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: A minor correction applied in Eq. 3

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2020, Issue 3, March 2020, 033C01

Theoretical justification and error analysis for slender body theory

Authors: Yoichiro Mori, Laurel Ohm, Daniel Spirn

Abstract: Slender body theory facilitates computational simulations of thin fibers immersed in a viscous fluid by approximating each fiber using only the geometry of the fiber centerline curve and the line force density along it. However, it has been unclear how well slender body theory actually approximates Stokes flow about a thin but truly three-dimensional fiber, in part due to the fact that simply pres… ▽ More Slender body theory facilitates computational simulations of thin fibers immersed in a viscous fluid by approximating each fiber using only the geometry of the fiber centerline curve and the line force density along it. However, it has been unclear how well slender body theory actually approximates Stokes flow about a thin but truly three-dimensional fiber, in part due to the fact that simply prescribing data along a one-dimensional curve does not result in a well-posed boundary value problem for the Stokes equations in $\mathbb{R}^3$. Here, we introduce a PDE problem to which slender body theory (SBT) provides an approximation, thereby placing SBT on firm theoretical footing. The slender body PDE is a new type of boundary value problem for Stokes flow where partial Dirichlet and partial Neumann conditions are specified everywhere along the fiber surface. Given only a 1D force density along a closed fiber, we show that the flow field exterior to the thin fiber is uniquely determined by imposing a fiber integrity condition: the surface velocity field on the fiber must be constant along cross sections orthogonal to the fiber centerline. Furthermore, a careful estimation of the residual, together with stability estimates provided by the PDE well-posedness framework, allow us to establish error estimates between the slender body approximation and the exact solution to the above problem. The error is bounded by an expression proportional to the fiber radius (up to logarithmic corrections) under mild regularity assumptions on the 1D force density and fiber centerline geometry. △ Less

Submitted 31 January, 2019; v1 submitted 30 June, 2018; originally announced July 2018.

From Electrodiffusion Theory to the Electrohydrodynamics of Leaky Dielectrics through the Weak Electrolyte Limit

Authors: Yoichiro Mori, Yuan-Nan Young

Abstract: The Taylor-Melcher (TM) model is the standard model for describing the dynamics of poorly conducting leaky dielectric fluids under an electric field. The TM model treats the fluids as Ohmic conductors, without modeling the underlying ion dynamics. On the other hand, electrodiffusion models, which have been successful in describing electrokinetic phenomena, incorporate ionic concentration dynamics.… ▽ More The Taylor-Melcher (TM) model is the standard model for describing the dynamics of poorly conducting leaky dielectric fluids under an electric field. The TM model treats the fluids as Ohmic conductors, without modeling the underlying ion dynamics. On the other hand, electrodiffusion models, which have been successful in describing electrokinetic phenomena, incorporate ionic concentration dynamics. Mathematical reconciliation of the electrodiffusion picture and the TM model has been a major issue for electrohydrodynamic theory. Here, we derive the TM model from an electrodiffusion model in which we explicitly model the electrochemistry of ion dissociation. We introduce salt dissociation reaction terms in the bulk electrodiffusion equations and take the limit in which the salt dissociation is weak; the assumption of weak dissociation corresponds to the fact that the TM model describes poor conductors. Together with the assumption that the Debye length is small, we derive the TM model with or without the surface charge convection term depending upon the scaling of relevant dimensionless parameters. An important quantity that emerges is the Galvani potential (GP), the jump in voltage across the liquid-liquid interface between the two leaky dielectric media; the GP arises as a natural consequence of the interfacial boundary conditions for the ionic concentrations, and is absent under certain parametric conditions. When the GP is absent, we recover the TM model. Our analysis also reveals the structure of the Debye layer at the liquid-liquid interface, which suggests how interfacial singularities may arise under strong imposed electric fields. In the presence of a non-zero GP, our model predicts that the liquid droplet will drift under an imposed electric field, the velocity of which is computed explicitly to leading order. △ Less

Submitted 10 July, 2018; v1 submitted 7 September, 2017; originally announced September 2017.

MSC Class: 76W05

MuSIC: delivering the world's most intense muon beam

Authors: S. Cook, R. D'Arcy, A. Edmonds, M. Fukuda, K. Hatanaka, Y. Hino, Y. Kuno, M. Lancaster, Y. Mori, T. Ogitsu, H. Sakamoto, A. Sato, N. H. Tran, N. M. Truong, M. Wing, A. Yamamoto, M. Yoshida

Abstract: A new muon beamline, muon science innovative channel (MuSIC), was set up at the Research Centre for Nuclear Physics (RCNP), Osaka University, in Osaka, Japan, using the 392 MeV proton beam impinging on a target. The production of an intense muon beam relies on the efficient capture of pions, which subsequently decay to muons, using a novel superconducting solenoid magnet system. After the pion-cap… ▽ More A new muon beamline, muon science innovative channel (MuSIC), was set up at the Research Centre for Nuclear Physics (RCNP), Osaka University, in Osaka, Japan, using the 392 MeV proton beam impinging on a target. The production of an intense muon beam relies on the efficient capture of pions, which subsequently decay to muons, using a novel superconducting solenoid magnet system. After the pion-capture solenoid the first $36^\circ$ of the curved muon transport line was commissioned and the muon flux was measured. In order to detect muons, a target of either copper or magnesium was placed to stop muons at the end of the muon beamline. Two stations of plastic scintillators located upstream and downstream from the muon target were used to reconstruct the decay spectrum of muons. In a complementary method to detect negatively-charged muons, the X-ray spectrum yielded by muonic atoms in the target were measured in a germanium detector. Measurements, at a proton beam current of 6 pA, yielded $(10.4 \pm 2.7) \times 10^5$ muons per Watt of proton beam power ($μ^+$ and $μ^-$), far in excess of other facilities. At full beam power (400 W), this implies a rate of muons of $(4.2 \pm 1.1) \times 10^8$ muons s$^{-1}$, amongst the highest in the world. The number of $μ^-$ measured was about a factor of 10 lower, again by far the most efficient muon beam produced. The set up is a prototype for future experiments requiring a high-intensity muon beam, such as a muon collider or neutrino factory, or the search for rare muon decays which would be a signature for phenomena beyond the Standard Model of particle physics. Such a muon beam can also be used in other branches of physics, nuclear and condensed matter, as well as other areas of scientific research. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 15 pages, 5 figures, submitted to Phys. Rev. ST-Acc. Beams

Rhythmomimetic drug delivery: modeling, analysis and numerical simulation

Authors: M. Carme Calderer, Yoichiro Mori, Ronald A. Siegel, Lingxing Yao

Abstract: We develop, analyse and numerically simulate a model of a prototype, glucose-driven, rhythmic drug delivery device, aimed at hormone therapies, and based on chemomechanical interaction in a polyelectrolyte gel membrane. The pH-driven interactions trigger volume phase transitions between the swollen and collapsed states of the gel.For a robust set of material parameters, we find a class of solution… ▽ More We develop, analyse and numerically simulate a model of a prototype, glucose-driven, rhythmic drug delivery device, aimed at hormone therapies, and based on chemomechanical interaction in a polyelectrolyte gel membrane. The pH-driven interactions trigger volume phase transitions between the swollen and collapsed states of the gel.For a robust set of material parameters, we find a class of solutions of the governing system that oscillate between such states, and cause the membrane to rhythmically swell, allowing for transport of the drug, fuel and reaction products across it, and collapse, hampering all transport across it. The frequency of the oscillations can be adjusted so that it matches the natural frequency of the hormone to be released. The work is linked to extensive laboratory experimental studies of the device built by Siegel's team. The thinness of the membrane and its clamped boundary, together with the homogeneously held conditions in the experimental apparatus, justify neglecting spatial dependence on the fields of the problem. Upon identifying the forces and energy relevant to the system, and taking into account its dissipative properties, we apply Rayleigh's variational principle to obtain the governing equations. The material assumptions guarantee the monotonicity of the system and lead to the existence of a three dimensional limit cycle. By scaling and asymptotic analysis, this limit cycle is found to be related to a two-dimensional one that encodes the volume phase transitions of the model. The identification of the relevant parameter set of the model is aided by a Hopf bifurcation study of steady state solutions. △ Less

Submitted 5 April, 2016; v1 submitted 6 January, 2015; originally announced December 2015.

MSC Class: 34C12; 37G15; 74B20; 82B26; 94C45

Characterization techniques for fixed-field alternating gradient accelerators and beam studies using the KURRI 150 MeV proton FFAG

Authors: S. L. Sheehy, D. J. Kelliher, S. Machida, C. Rogers, C. R. Prior, L. Volat, M. Haj Tahar, Y. Ishi, Y. Kuriyama, M. Sakamoto, T. Uesugi, Y. Mori

Abstract: In this paper we describe the methods and tools used to characterize a 150 MeV proton scaling Fixed Field Alternating Gradient (FFAG) accelerator at Kyoto University Research Reactor Institute. Many of the techniques used are unique to this class of machine and are thus of relevance to any future FFAG accelerator. For the first time we detail systematic studies undertaken to improve the beam quali… ▽ More In this paper we describe the methods and tools used to characterize a 150 MeV proton scaling Fixed Field Alternating Gradient (FFAG) accelerator at Kyoto University Research Reactor Institute. Many of the techniques used are unique to this class of machine and are thus of relevance to any future FFAG accelerator. For the first time we detail systematic studies undertaken to improve the beam quality of the FFAG. The control of beam quality in this manner is crucial to demonstrating high power operation of FFAG accelerators in future. △ Less

Submitted 6 May, 2016; v1 submitted 26 October, 2015; originally announced October 2015.

A Multidomain Model for Ionic Electrodiffusion and Osmosis with an Application to Cortical Spreading Depression

Authors: Yoichiro Mori

Abstract: Ionic electrodiffusion and osmotic water flow are central processes in many physiological systems. We formulate a system of partial differential equations that governs ion movement and water flow in biological tissue. A salient feature of this model is that it satisfies a free energy identity, ensuring the thermodynamic consistency of the model. A numerical scheme is developed for the model in one… ▽ More Ionic electrodiffusion and osmotic water flow are central processes in many physiological systems. We formulate a system of partial differential equations that governs ion movement and water flow in biological tissue. A salient feature of this model is that it satisfies a free energy identity, ensuring the thermodynamic consistency of the model. A numerical scheme is developed for the model in one spatial dimension and is applied to a model of cortical spreading depression, a propagating breakdown of ionic and cell volume homeostasis in the brain. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: submitted for publication, Aug. 28, 2014

MSC Class: 92C20; 92C30; 78A70

Accelerator system for the PRISM based muon to electron conversion experiment

Authors: A. Alekou, R. Appleby, M. Aslaninejad, R. J. Barlow, R. Chudzinski K. M. Hock, J. Garland, L. J. Jenner, D. J. Kelliher, Y. Kuno, A. Kurup, J-B. Lagrange, M. Lancaster, S. Machida, Y. Mori, B. Muratori, C. Ohmori, H. Owen, J. Pasternak, T. Planche, C. Prior, A. Sato, Y. Shi, S. Smith, Y. Uchida, H. Witte , et al. (1 additional authors not shown)

Abstract: The next generation of lepton flavor violation experiments need high intensity and high quality muon beams. Production of such beams requires sending a short, high intensity proton pulse to the pion production target, capturing pions and collecting the resulting muons in the large acceptance transport system. The substantial increase of beam quality can be obtained by applying the RF phase rotatio… ▽ More The next generation of lepton flavor violation experiments need high intensity and high quality muon beams. Production of such beams requires sending a short, high intensity proton pulse to the pion production target, capturing pions and collecting the resulting muons in the large acceptance transport system. The substantial increase of beam quality can be obtained by applying the RF phase rotation on the muon beam in the dedicated FFAG ring, which was proposed for the PRISM project.This allows to reduce the momentum spread of the beam and to purify from the unwanted components like pions or secondary protons. A PRISM Task Force is addressing the accelerator and detector issues that need to be solved in order to realize the PRISM experiment. The parameters of the required proton beam, the principles of the PRISM experiment and the baseline FFAG design are introduced. The spectrum of alternative designs for the PRISM FFAG ring are shown. Progress on ring main systems like injection and RF are presented. The current status of the study and its future directions are discussed. △ Less

Submitted 2 October, 2013; originally announced October 2013.

Comments: Studies performed within the PRISM Task Force initiative

nuSTORM - Neutrinos from STORed Muons: Proposal to the Fermilab PAC

Authors: D. Adey, S. K. Agarwalla, C. M. Ankenbrandt, R. Asfandiyarov, J. J. Back, G. Barker, E. Baussan, R. Bayes, S. Bhadra, V. Blackmore, A. Blondel, S. A. Bogacz, C. Booth, S. B. Boyd, A. Bravar, S. J. Brice, A. D. Bross, F. Cadoux, H. Cease, A. Cervera, J. Cobb, D. Colling, P. Coloma, L. Coney, A. Dobbs , et al. (88 additional authors not shown)

Abstract: The nuSTORM facility has been designed to deliver beams of electron neutrinos and muon neutrinos (and their anti-particles) from the decay of a stored muon beam with a central momentum of 3.8 GeV/c and a momentum acceptance of 10%. The facility is unique in that it will: 1. Allow searches for sterile neutrinos of exquisite sensitivity to be carried out; 2. Serve future long- and short-baseline neu… ▽ More The nuSTORM facility has been designed to deliver beams of electron neutrinos and muon neutrinos (and their anti-particles) from the decay of a stored muon beam with a central momentum of 3.8 GeV/c and a momentum acceptance of 10%. The facility is unique in that it will: 1. Allow searches for sterile neutrinos of exquisite sensitivity to be carried out; 2. Serve future long- and short-baseline neutrino-oscillation programs by providing definitive measurements of electron neutrino and muon neutrino scattering cross sections off nuclei with percent-level precision; and 3. Constitutes the crucial first step in the development of muon accelerators as a powerful new technique for particle physics. The document describes the facility in detail and demonstrates its physics capabilities. This document was submitted to the Fermilab Physics Advisory Committee in consideration for Stage I approval. △ Less

Submitted 31 July, 2013; originally announced August 2013.

Neutrinos from Stored Muons nuSTORM: Expression of Interest

Authors: D. Adey, S. K. Agarwalla, C. M. Ankenbrandt, R. Asfandiyarov, J. J. Back, G. Barker, E. Baussan, R. Bayes, S. Bhadra, V. Blackmore, A. Blondel, S. A. Bogacz, C. Booth, S. B. Boyd, A. Bravar, S. J. Brice, A. D. Bross, F. Cadoux, H. Cease, A. Cervera, J. Cobb, D. Colling, L. Coney, A. Dobbs, J. Dobson , et al. (84 additional authors not shown)

Abstract: The nuSTORM facility has been designed to deliver beams of electron and muon neutrinos from the decay of a stored muon beam with a central momentum of 3.8 GeV/c and a momentum spread of 10%. The facility is unique in that it will: serve the future long- and short-baseline neutrino-oscillation programmes by providing definitive measurements of electron-neutrino- and muon-neutrino-nucleus cross sect… ▽ More The nuSTORM facility has been designed to deliver beams of electron and muon neutrinos from the decay of a stored muon beam with a central momentum of 3.8 GeV/c and a momentum spread of 10%. The facility is unique in that it will: serve the future long- and short-baseline neutrino-oscillation programmes by providing definitive measurements of electron-neutrino- and muon-neutrino-nucleus cross sections with percent-level precision; allow searches for sterile neutrinos of exquisite sensitivity to be carried out; and constitute the essential first step in the incremental development of muon accelerators as a powerful new technique for particle physics. Of the world's proton-accelerator laboratories, only CERN and FNAL have the infrastructure required to mount nuSTORM. Since no siting decision has yet been taken, the purpose of this Expression of Interest (EoI) is to request the resources required to: investigate in detail how nuSTORM could be implemented at CERN; and develop options for decisive European contributions to the nuSTORM facility and experimental programme wherever the facility is sited. The EoI defines a two-year programme culminating in the delivery of a Technical Design Report. △ Less

Submitted 7 May, 2013; originally announced May 2013.

Comments: 59 pages; 24 figures; 5 tables

Report number: CERN-SPSC-2013-015 / SPSC-EOI-009

nuSTORM: Neutrinos from STORed Muons

Authors: P. Kyberd, D. R. Smith, L. Coney, S. Pascoli, C. Ankenbrandt, S. J. Brice, A. D. Bross, H. Cease, J. Kopp, N. Mokhov, J. Morfin, D. Neuffer, M. Popovic, P. Rubinov, S. Striganov, A. Blondel, A. Bravar, E. Noah, R. Bayes, F. J. P. Soler, A. Dobbs, K. Long, J. Pasternak, E. Santos, M. O. Wascko , et al. (13 additional authors not shown)

Abstract: The results of LSND and MiniBooNE, along with the recent papers on a possible reactor neutrino flux anomaly give tantalizing hints of new physics. Models beyond the neutrino-SM have been developed to explain these results and involve one or more additional neutrinos that are non-interacting or "sterile." Neutrino beams produced from the decay of muons in a racetrack-like decay ring provide a power… ▽ More The results of LSND and MiniBooNE, along with the recent papers on a possible reactor neutrino flux anomaly give tantalizing hints of new physics. Models beyond the neutrino-SM have been developed to explain these results and involve one or more additional neutrinos that are non-interacting or "sterile." Neutrino beams produced from the decay of muons in a racetrack-like decay ring provide a powerful way to study this potential new physics. In this Letter of Intent, we describe a facility, nuSTORM, "Neutrinos from STORed Muons," and an appropriate far detector for neutrino oscillation searches at short baseline. We present sensitivity plots that indicated that this experimental approach can provide over 10 sigma confirmation or rejection of the LSND/MinBooNE results. In addition we indicate how the facility can be used to make precision neutrino interaction cross section measurements important to the next generation of long-baseline neutrino oscillation experiments. △ Less

Submitted 1 June, 2012; originally announced June 2012.

Report number: Fermilab P-1028

MICE: the Muon Ionization Cooling Experiment. Step I: First Measurement of Emittance with Particle Physics Detectors

Authors: U. Bravar, M. Bogomilov, Y. Karadzhov, D. Kolev, I. Russinov, R. Tsenov, L. Wang, F. Y. Xu, S. X. Zheng, R. Bertoni, M. Bonesini, R. Mazza, V. Palladino, G. Cecchet, A. de Bari, M. Capponi, A. Iaciofano, D. Orestano, F. Pastore, L. Tortora, S. Ishimoto, S. Suzuki, K. Yoshimura, Y. Mori, Y. Kuno , et al. (123 additional authors not shown)

Abstract: The Muon Ionization Cooling Experiment (MICE) is a strategic R&D project intended to demonstrate the only practical solution to providing high brilliance beams necessary for a neutrino factory or muon collider. MICE is under development at the Rutherford Appleton Laboratory (RAL) in the United Kingdom. It comprises a dedicated beamline to generate a range of input muon emittances and momenta, with… ▽ More The Muon Ionization Cooling Experiment (MICE) is a strategic R&D project intended to demonstrate the only practical solution to providing high brilliance beams necessary for a neutrino factory or muon collider. MICE is under development at the Rutherford Appleton Laboratory (RAL) in the United Kingdom. It comprises a dedicated beamline to generate a range of input muon emittances and momenta, with time-of-flight and Cherenkov detectors to ensure a pure muon beam. The emittance of the incoming beam will be measured in the upstream magnetic spectrometer with a scintillating fiber tracker. A cooling cell will then follow, alternating energy loss in Liquid Hydrogen (LH2) absorbers to RF cavity acceleration. A second spectrometer, identical to the first, and a second muon identification system will measure the outgoing emittance. In the 2010 run at RAL the muon beamline and most detectors were fully commissioned and a first measurement of the emittance of the muon beam with particle physics (time-of-flight) detectors was performed. The analysis of these data was recently completed and is discussed in this paper. Future steps for MICE, where beam emittance and emittance reduction (cooling) are to be measured with greater accuracy, are also presented. △ Less

Submitted 30 July, 2013; v1 submitted 9 October, 2011; originally announced October 2011.

Comments: Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011

A Model of Electrodiffusion and Osmotic Water Flow and its Energetic Structure

Authors: Yoichiro Mori, Chun Liu, Robert S. Eisenberg

Abstract: We introduce a model for ionic electrodiffusion and osmotic water flow through cells and tissues. The model consists of a system of partial differential equations for ionic concentration and fluid flow with interface conditions at deforming membrane boundaries. The model satisfies a natural energy equality, in which the sum of the entropic, elastic and electrostatic free energies are dissipated th… ▽ More We introduce a model for ionic electrodiffusion and osmotic water flow through cells and tissues. The model consists of a system of partial differential equations for ionic concentration and fluid flow with interface conditions at deforming membrane boundaries. The model satisfies a natural energy equality, in which the sum of the entropic, elastic and electrostatic free energies are dissipated through viscous, electrodiffusive and osmotic flows. We discuss limiting models when certain dimensionless parameters are small. Finally, we develop a numerical scheme for the one-dimensional case and present some simple applications of our model to cell volume control. △ Less

Submitted 26 January, 2011; originally announced January 2011.

MSC Class: 92C30 (primary) 35Q92 (secondary)

Enhanced Sampling Algorithms

Authors: Ayori Mitsutake, Yoshiharu Mori, Yuko Okamoto

Abstract: In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped in states of these energy local minima. Enhanced conformational sampling techniques are thus in great de… ▽ More In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped in states of these energy local minima. Enhanced conformational sampling techniques are thus in great demand. A simulation in generalized ensemble performs a random walk in potential energy space and can overcome this difficulty. From only one simulation run, one can obtain canonical-ensemble averages of physical quantities as functions of temperature by the single-histogram and/or multiple-histogram reweighting techniques. In this article we review uses of the generalized-ensemble algorithms in biomolecular systems. Three well-known methods, namely, multicanonical algorithm, simulated tempering, and replica-exchange method, are described first. Both Monte Carlo and molecular dynamics versions of the algorithms are given. We then present various extensions of these three generalized-ensemble algorithms. The effectiveness of the methods is tested with short peptide and protein systems. △ Less

Submitted 28 December, 2010; originally announced December 2010.

Comments: 42 pages, 7 figures, to be published in "Biomolecular simulations: methods and protocols" edited by Luca Monticelli and Emppu Salonen (Humana Press, 2011)

From Three-Dimensional Electrophysiology to the Cable Model: an Asymptotic Study

Authors: Yoichiro Mori

Abstract: Cellular electrophysiology is often modeled using the cable equations. The cable model can only be used when ionic concentration effects and three dimensional geometry effects are negligible. The Poisson model, in which the electrostatic potential satisfies the Poisson equation and the ionic concentrations satisfy the drift-diffusion equation, is a system of equations that can incorporate such e… ▽ More Cellular electrophysiology is often modeled using the cable equations. The cable model can only be used when ionic concentration effects and three dimensional geometry effects are negligible. The Poisson model, in which the electrostatic potential satisfies the Poisson equation and the ionic concentrations satisfy the drift-diffusion equation, is a system of equations that can incorporate such effects. The Poisson model is unfortunately prohibitively expensive for numerical computation because of the presence of thin space charge layers at internal membrane boundaries. As a computationally efficient and biophysically natural alternative, we introduce the electroneutral model in which the Poisson equation is replaced by the electroneutrality condition and the presence of the space charge layer is incorporated in boundary conditions at the membrane interfaces. We use matched asymptotics and numerical computations to show that the electroneutral model provides an excellent approximation to the Poisson model. Further asymptotic calculations illuminate the relationship of the electroneutral or Poisson models with the cable model, and reveal the presence of a hierarchy of electrophysiology models. △ Less

Submitted 25 January, 2009; originally announced January 2009.