-
Diffuse optical tomography in time domain with the inverse Rytov series
Authors:
Chi Zhang,
Manabu Machida
Abstract:
The Rytov approximation has been commonly used to obtain reconstructed images for optical tomography. However, the method requires linearization of the nonlinear inverse problem. Here, we demonstrate nonlinear Rytov approximations by developing the inverse Rytov series for the time-dependent diffusion equation. The method is experimentally verified.
The Rytov approximation has been commonly used to obtain reconstructed images for optical tomography. However, the method requires linearization of the nonlinear inverse problem. Here, we demonstrate nonlinear Rytov approximations by developing the inverse Rytov series for the time-dependent diffusion equation. The method is experimentally verified.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Exact Calculation of Nonideal Fields Demonstrates Dominance of Injection in Relativistic Reconnection
Authors:
Samuel R. Totorica,
Seiji Zenitani,
Shuichi Matsukiyo,
Mami Machida,
Kazuhiro Sekiguchi,
Amitava Bhattacharjee
Abstract:
Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-…
▽ More
Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-consistently in kinetic particle-in-cell simulations of reconnection, we conclusively show the dominant role of the nonideal field for the early stage of energization known as injection. The importance of the nonideal field increases with magnetization, guide field, and in three-dimensions, indicating its general importance for reconnection in natural astrophysical systems. We obtain the statistical properties of the injection process from the simulations, paving the way for the development of extended MHD models capable of accurately modeling particle acceleration in large-scale systems. The novel analysis method developed in this study can be applied broadly to give new insight into a wide range of processes in plasma physics.
△ Less
Submitted 15 September, 2023;
originally announced September 2023.
-
Coarse-graining of the perplexity for the spatial distribution of molecules
Authors:
Lili Xu,
Chi Zhang,
Soho Oyama,
Manabu Machida,
Tomoaki Kahyo,
Mitsutoshi Setou
Abstract:
Biological tissue consists of various molecules. Instead of focusing on a particular molecule, we consider the Shannon entropy which is calculated from the abundance of different molecules at each spot in the tissue. The spatial distribution of the Shannon entropy is of interest. In this paper, we first obtain the heat map of perplexity, whose logarithm is the entropy. To characterize the spatial…
▽ More
Biological tissue consists of various molecules. Instead of focusing on a particular molecule, we consider the Shannon entropy which is calculated from the abundance of different molecules at each spot in the tissue. The spatial distribution of the Shannon entropy is of interest. In this paper, we first obtain the heat map of perplexity, whose logarithm is the entropy. To characterize the spatial variety of molecules, we propose a scalar $k$ that is concerned with the coarse-graining of the perplexity heat map. To verify the usefulness of the number, experiments with MSI (mass spectrometry imaging) were performed for mouse kidneys. We found that $k$ has large values in the renal pelvis area, cortex area, veins, and arteries in the mouse kidney, whereas fractal dimensions fail to distinguish those regions.
△ Less
Submitted 10 December, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
The depth of the banana and the impulse stripe illumination for diffuse optical tomography
Authors:
Manabu Machida,
Keita Osada,
Keiichiro Kagawa
Abstract:
The stripe illumination lies between the illumination in the spatial-frequency domain and the point illumination. Although the stripe illumination has a periodic structure as the illumination in the spatial-frequency domain, light from the stripe illumination can reach deep regions in biological tissue since it can be regarded as an array of point illuminations. For a pair of a source and a detect…
▽ More
The stripe illumination lies between the illumination in the spatial-frequency domain and the point illumination. Although the stripe illumination has a periodic structure as the illumination in the spatial-frequency domain, light from the stripe illumination can reach deep regions in biological tissue since it can be regarded as an array of point illuminations. For a pair of a source and a detector, the shape of light paths which connect the source and detector is called the banana shape. First, we investigate the depth of the banana. In the case of the zero boundary condition, we found that the depth of the center of the banana is about $0.2d_{\rm SD}$ for typical optical parameters, where $d_{\rm SD}$ is the distance between the source and detector on the boundary. In general, the depth depends on the absorption and diffusion coefficients, and the ratio of refractive indices on the boundary. Next, we perform diffuse optical tomography for the stripe illumination against forward data taken by Monte Carlo simulation. We consider an impulse illumination of the shape of a stripe. By this time-resolved measurement, the absorption coefficient of a target is reconstructed.
△ Less
Submitted 13 May, 2024; v1 submitted 4 August, 2022;
originally announced August 2022.
-
Shape-constrained reconstruction in diffuse optical tomography by simulated annealing
Authors:
Tetsuya Mimura,
Yu Jiang,
Norikazu Todoroki,
Manabu Machida
Abstract:
When the inverse problem of diffuse optical tomography (DOT) is solved with the Born or Rytov approximation, the size of the matrix of the linear inverse problem becomes large if the volume (or area) of the domain in biological tissue used for reconstruction is large. The number of unknown parameters in DOT is reduced when constraints about the shape of a target are imposed for the inverse problem…
▽ More
When the inverse problem of diffuse optical tomography (DOT) is solved with the Born or Rytov approximation, the size of the matrix of the linear inverse problem becomes large if the volume (or area) of the domain in biological tissue used for reconstruction is large. The number of unknown parameters in DOT is reduced when constraints about the shape of a target are imposed for the inverse problem. Due to such constraints, the inverse problem becomes nonlinear even when the (first) Born or Rytov approximation is employed. We solve this nonlinear inverse problem by the simulated annealing, which is not trapped by local minima of the cost function.
△ Less
Submitted 7 July, 2021;
originally announced July 2021.
-
Simulation analysis of the Compton-to-peak method for quantifying radiocesium deposition quantities
Authors:
Alex Malins,
Kotaro Ochi,
Masahiko Machida,
Yukihisa Sanada
Abstract:
Compton-to-peak analysis is a method for selecting coefficients for converting count rates measured with in situ gamma ray spectrometry to radioactivity concentrations of $^{134}$Cs and $^{137}$Cs in the environment. Compton-to-peak analysis is based on the count rate ratio between the spectral regions containing scattered gamma rays to the primary $^{134}$Cs and $^{137}$Cs photopeaks (known as th…
▽ More
Compton-to-peak analysis is a method for selecting coefficients for converting count rates measured with in situ gamma ray spectrometry to radioactivity concentrations of $^{134}$Cs and $^{137}$Cs in the environment. Compton-to-peak analysis is based on the count rate ratio between the spectral regions containing scattered gamma rays to the primary $^{134}$Cs and $^{137}$Cs photopeaks (known as the Compton-to-peak ratio - RCP). RCP changes with the vertical distribution of $^{134}$Cs and $^{137}$Cs within the ground. Inferring this distribution enables the selection of appropriate count rate to activity concentration conversion coefficients. Here PHITS was used to simulate the dependency of RCP on different vertical distributions of $^{134}$Cs and $^{137}$Cs within the ground. A model was created of a LaBr$_3$(Ce) detector used in drone helicopter aerial surveys in Fukushima Prefecture. The model was verified by comparing simulated gamma ray spectra to measurements from test sources. Simulations were performed for the infinite half-space geometry to calculate the dependency of RCP on the mass depth distribution (exponential or uniform) of $^{134}$Cs and $^{137}$Cs within the ground, and on the altitude of the detector above the ground. The calculations suggest that the sensitivity of the Compton-to-peak method is greatest for the initial period following nuclear fallout when $^{134}$Cs and $^{137}$Cs are located close to the ground surface, and for aerial surveys conducted at low altitudes. This is because the relative differences calculated between RCP with respect to changes in the mass depth distribution were largest for these two cases. Data on the measurement height above and on the $^{134}$Cs to $^{137}$Cs activity ratio is necessary for applying the Compton-to-peak method to determine the distribution and radioactivity concentration of $^{134}$Cs and $^{137}$Cs within the ground.
△ Less
Submitted 8 November, 2020;
originally announced November 2020.
-
Homogeneity of neutron transmission imaging over a large sensitive area with a four-channel superconducting detector
Authors:
The Dang Vu,
Hiroaki Shishido,
Kenji M. Kojima,
Tomio Koyama,
Kenichi Oikawa,
Masahide Harada,
Shigeyuki Miyajima,
Takayuki Oku,
Kazuhiko Soyama,
Kazuya Aizawa,
Mutsuo Hidaka,
Soh Y. Suzuki,
Manobu M. Tanaka,
Alex Malins,
Masahiko Machida,
Takekazu Ishida
Abstract:
We previously proposed a method to detect neutrons by using a current-biased kinetic inductance detector (CB-KID), where neutrons are converted into charged particles using a 10B conversion layer. The charged particles are detected based on local changes in kinetic inductance of X and Y superconducting meanderlines under a modest DC bias current. The system uses a delay-line method to locate the p…
▽ More
We previously proposed a method to detect neutrons by using a current-biased kinetic inductance detector (CB-KID), where neutrons are converted into charged particles using a 10B conversion layer. The charged particles are detected based on local changes in kinetic inductance of X and Y superconducting meanderlines under a modest DC bias current. The system uses a delay-line method to locate the positions of neutron-10B reactions by acquiring the four arrival timestamps of signals that propagate from hot spots created by a passing charged particle to the end electrodes of the meanderlines. Unlike conventional multi-pixel imaging systems, the CB-KID system performs high spatial resolution imaging over a 15 mm x 15 mm sensitive area using only four channel readouts. Given the large sensitive area, it is important to check the spatial homogeneity and linearity of detected neutron positions when imaging with CB-KID. To this end we imaged a pattern of 10B dot absorbers with a precise dot pitch to investigate the spatial homogeneity of the detector. We confirmed the spatial homogeneity of detected dot positions based on the distribution of measured dot pitches across the sensitive area of the detector. We demonstrate potential applications of the system by taking a clear transmission image of tiny metallic screws and nuts and a ladybug. The image was useful for characterizing the ladybug noninvasively. Detection efficiencies were low when the detector was operated at 4 K, so we plan to explore raising the operating temperature towards the critical temperature of the detector as a means to improve counting rates.
△ Less
Submitted 14 October, 2020;
originally announced October 2020.
-
Properties of DBD Plasma Jets using Powered Electrode With and Without Contact with the Plasma
Authors:
Fellype do Nascimento,
Konstantin Kostov,
Munemasa Machida,
Alexander Flacker
Abstract:
An experimental investigation comparing the properties of plasma jets in dielectric barrier discharge (DBD) configurations using a powered electrode with and without a dielectric barrier, while keeping a second dielectric barrier over the grounded electrode, is reported in this work. For this purpose, two different power sources were used to produce the plasma jets, with one of them producing a pu…
▽ More
An experimental investigation comparing the properties of plasma jets in dielectric barrier discharge (DBD) configurations using a powered electrode with and without a dielectric barrier, while keeping a second dielectric barrier over the grounded electrode, is reported in this work. For this purpose, two different power sources were used to produce the plasma jets, with one of them producing a pulsed high-voltage (HV) output and the other one producing a pulse-like HV output, which consists of a damped sine HV waveform. Measurements of plasma parameters were performed for both configurations using argon and helium as working gases. As a result, if the pulsed power source is used, significant differences were found in discharge power ($P_{plasma}$), rotational and vibrational temperatures ($T_r$ and $T_v$, respectively) when switching from one configuration to the other. On the other hand, using the pulse-like HV only the $P_{plasma}$ parameter presented significant differences when switching the electrode's configuration. For the pulsed source it has been observed that despite the remarkable increase in $P_{plasma}$ when changing from the double barrier configuration to the single barrier one, the values obtained for $T_r$ and $T_v$ also increased, but not in the same proportion as the increase in $P_{plasma}$, which suggests a non-linear dependency between temperatures and discharge power in the plasma jet. As an example for application of plasmas in both configurations, tests in an attempt to remove copper films deposited on alumina substrates were performed and, as a result, there was significant material removal only when the powered electrode was in contact with the plasma. As a general conclusion, if higher power is really required for some application it is better to use this configuration.
△ Less
Submitted 14 April, 2021; v1 submitted 21 August, 2020;
originally announced August 2020.
-
Linear transport in porous media
Authors:
Kenji Amagai,
Yuko Hatano,
Manabu Machida
Abstract:
The linear transport theory is developed to describe the time dependence of the number density of tracer particles in porous media. The advection is taken into account. The transport equation is numerically solved by the analytical discrete ordinates method. For the inverse Laplace transform, the double-exponential formula is employed.
The linear transport theory is developed to describe the time dependence of the number density of tracer particles in porous media. The advection is taken into account. The transport equation is numerically solved by the analytical discrete ordinates method. For the inverse Laplace transform, the double-exponential formula is employed.
△ Less
Submitted 25 May, 2020;
originally announced May 2020.
-
Decay behavior and optical parameter identification for spatial-frequency domain imaging by the radiative transport equation
Authors:
Manabu Machida,
Yoko Hoshi,
Keiichiro Kagawa,
Kazuki Takada
Abstract:
The decay behavior of the specific intensity is studied for the spatial-frequency domain imaging (SFDI). It is shown using the radiative transport equation that the decay is given by a superposition of different decay modes, and the decay rates of these modes are determined by spatial frequencies and Case's eigenvalues. This explains why SFDI can focus on shallow regions. The fact that light with…
▽ More
The decay behavior of the specific intensity is studied for the spatial-frequency domain imaging (SFDI). It is shown using the radiative transport equation that the decay is given by a superposition of different decay modes, and the decay rates of these modes are determined by spatial frequencies and Case's eigenvalues. This explains why SFDI can focus on shallow regions. The fact that light with nonzero spatial frequency rapidly decays makes it possible to exclusively extract optical properties of the top layer of a layered medium. We determine optical properties of the top layer of a solid phantom. This measurement is verified with different layered media of numerical phantoms.
△ Less
Submitted 17 August, 2022; v1 submitted 20 May, 2020;
originally announced May 2020.
-
Diffuse optical tomography by simulated annealing via a spin Hamiltonian
Authors:
Yu Jiang,
Manabu Machida,
Norikazu Todoroki
Abstract:
Diffuse optical tomography (DOT) is an imaging modality which uses near-infrared light. Although iterative numerical schemes are commonly used for its inverse problem, correct solutions are not obtained unless good initial guesses are chosen. We propose a numerical scheme of DOT which works even when good initial guesses of optical parameters are not available. We use simulated annealing (SA) whic…
▽ More
Diffuse optical tomography (DOT) is an imaging modality which uses near-infrared light. Although iterative numerical schemes are commonly used for its inverse problem, correct solutions are not obtained unless good initial guesses are chosen. We propose a numerical scheme of DOT which works even when good initial guesses of optical parameters are not available. We use simulated annealing (SA) which is a method of the Markov-chain Monte Carlo. To implement SA for DOT, a spin Hamiltonian is introduced in the cost function, and the Metropolis algorithm or single-component Metropolis-Hastings algorithm is used. By numerical experiments, it is shown that an initial random spin configuration is brought to a converged configuration by SA and targets in the medium are reconstructed. The proposed numerical method solves the inverse problem for DOT by finding the ground state of a spin Hamiltonian with SA.
△ Less
Submitted 3 June, 2021; v1 submitted 24 December, 2019;
originally announced December 2019.
-
Monte Carlo Radiation Transport Modelling of the Current-Biased Kinetic Inductance Detector
Authors:
Alex Malins,
Masahiko Machida,
The Dang Vu,
Kazuya Aizawa,
Takekazu Ishida
Abstract:
Radiation transport simulations were used to analyse neutron imaging with the current-biased kinetic inductance detector (CB-KID). The PHITS Monte Carlo code was applied for simulating neutron, $^{4}$He, $^{7}$Li, photon and electron transport, $^{10}$B(n,$α$)$^{7}$Li reactions, and energy deposition by particles within CB-KID. Slight blurring in simulated CB-KID images originated $^{4}$He and…
▽ More
Radiation transport simulations were used to analyse neutron imaging with the current-biased kinetic inductance detector (CB-KID). The PHITS Monte Carlo code was applied for simulating neutron, $^{4}$He, $^{7}$Li, photon and electron transport, $^{10}$B(n,$α$)$^{7}$Li reactions, and energy deposition by particles within CB-KID. Slight blurring in simulated CB-KID images originated $^{4}$He and $^{7}$Li ions spreading out in random directions from the $^{10}$B conversion layer in the detector prior to causing signals in the $X$ and $Y$ superconducting Nb nanowire meander lines. 478 keV prompt gamma rays emitted by $^{7}$Li nuclei from neutron-$^{10}$B reactions had negligible contribution to the simulated CB-KID images. Simulated neutron images of $^{10}$B dot arrays indicate that sub 10 $μ$m resolution imaging should be feasible with the current CB-KID design. The effect of the geometrical structure of CB-KID on the intrinsic detection efficiency was calculated from the simulations. An analytical equation was then developed to approximate this contribution to the detection efficiency. Detection efficiencies calculated in this study are upper bounds for the reality as the effects of detector temperature, the bias current, signal processing and dead-time losses were not taken into account. The modelling strategies employed in this study could be used to evaluate modifications to the CB-KID design prior to actual fabrication and testing, conveying a time and cost saving.
△ Less
Submitted 26 November, 2019; v1 submitted 26 November, 2019;
originally announced November 2019.
-
The linear Boltzmann equation in column experiments of porous media
Authors:
Kenji Amagai,
Motoko Yamakawa,
Manabu Machida,
Yuko Hatano
Abstract:
The use of the linear Boltzmann equation is proposed for transport in porous media in a column. By column experiments, we show that the breakthrough curve is reproduced by the linear Boltzmann equation. The advection-diffusion equation is derived from the linear Boltzmann equation in the asymptotic limit of large propagation distance and long time.
The use of the linear Boltzmann equation is proposed for transport in porous media in a column. By column experiments, we show that the breakthrough curve is reproduced by the linear Boltzmann equation. The advection-diffusion equation is derived from the linear Boltzmann equation in the asymptotic limit of large propagation distance and long time.
△ Less
Submitted 2 August, 2019;
originally announced August 2019.
-
Four-electrodes DBD plasma jet device with additional floating electrode
Authors:
Fellype do Nascimento,
Munemasa Machida,
Konstantin G. Kostov,
Stanislav Moshkalev
Abstract:
A Dielectric Barrier Discharge (DBD) plasma jet in a four electrodes configuration was investigated in order to improve the discharge parameters, such as, plasma power and rotational and vibrational temperatures of molecular species in the plasma plume. The improvement attempts were made by introducing an auxiliary floating electrode in a form of a metallic pin inside the DBD device. That piece wa…
▽ More
A Dielectric Barrier Discharge (DBD) plasma jet in a four electrodes configuration was investigated in order to improve the discharge parameters, such as, plasma power and rotational and vibrational temperatures of molecular species in the plasma plume. The improvement attempts were made by introducing an auxiliary floating electrode in a form of a metallic pin inside the DBD device. That piece was placed near the bottom of the main device, centered in relation to the four powered electrodes, which were covered with a dielectric material. By using metallic pins with different lengths, it was observed that there were considerable variations of the plasma parameters as a function of the pin length. Two carrier gases were tested: argon and helium. With helium as the working gas, it was found that there is an optimal pin length that maximizes the plasma power and its vibrational temperature. In addition, it was verified that for the pin of optimum length the relative intensity of light emissions from OH and NO species achieved higher values than in other conditions studied.
△ Less
Submitted 27 January, 2020; v1 submitted 19 July, 2019;
originally announced July 2019.
-
Fast and robust reconstruction algorithm for fluorescence diffuse optical tomography assuming a cuboid target
Authors:
Chunlong Sun,
Yu Jiang,
Jijun Liu,
Manabu Machida,
Gen Nakamura,
Goro Nishimura
Abstract:
A fast and robust algorithm for fluorescence diffuse optical tomography is proposed. We identify the location of a fluorescence target by assuming a cuboid. The proposed numerical method is verified by a numerical experiment and an in vivo experiment.
A fast and robust algorithm for fluorescence diffuse optical tomography is proposed. We identify the location of a fluorescence target by assuming a cuboid. The proposed numerical method is verified by a numerical experiment and an in vivo experiment.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
The role of vibrational temperature variations in a pulsed dielectric barrier discharge plasma device
Authors:
Fellype do Nascimento,
Stanislav Moshkalev,
Munemasa Machida
Abstract:
The use of dielectric barrier discharge (DBD) plasmas has become a practical way to carry out surface treatment and one seeks to do it in a more efficient way, which requires to have control of the plasma parameters like rotational and vibrational temperatures (T_rot and T_vib). Since the T_vib of an atmospheric pressure plasma jet is an important parameter related to improvement of surface treatm…
▽ More
The use of dielectric barrier discharge (DBD) plasmas has become a practical way to carry out surface treatment and one seeks to do it in a more efficient way, which requires to have control of the plasma parameters like rotational and vibrational temperatures (T_rot and T_vib). Since the T_vib of an atmospheric pressure plasma jet is an important parameter related to improvement of surface treatments, in this work it was analyzed two methods to increase the values of T_vib in the DBD plasma. One of the methods is to reduce the exit size of the DBD reactor, and the other is to change the gas flow rate. Explanations for the reasons that lead to the increment of the vibrational temperature are proposed in both cases.
△ Less
Submitted 30 November, 2018;
originally announced November 2018.
-
A hybrid inversion scheme combining Markov chain Monte Carlo and iterative methods for determining optical properties of random media
Authors:
Yu Jiang,
Yoko Hoshi,
Manabu Machida,
Gen Nakamura
Abstract:
Near-infrared spectroscopy (NIRS) including diffuse optical tomography is an imaging modality which makes use of diffuse light propagation in random media. When optical properties of a random medium is investigated from boundary measurements of reflected or transmitted light, iterative inversion schemes such as the Levenberg-Marquardt algorithm are known to fail when initial guesses are not close…
▽ More
Near-infrared spectroscopy (NIRS) including diffuse optical tomography is an imaging modality which makes use of diffuse light propagation in random media. When optical properties of a random medium is investigated from boundary measurements of reflected or transmitted light, iterative inversion schemes such as the Levenberg-Marquardt algorithm are known to fail when initial guesses are not close to the true value of the coefficient to be reconstructed. In this paper, we investigate how this weakness of iterative schemes is overcome by the use of Markov chain Monte Carlo. Using time-resolved measurements performed against a polyurethane-based phantom, we present a case that the Levenberg-Marquardt algorithm fails to work but the proposed hybrid method works well. Then with a toy model of diffuse optical tomography we illustrate that the evenberg-Marquardt method fails when it is trapped by a local minimum but the hybrid method can escape from local minima by using the Metropolis-Hastings Markov chain Monte Carlo algorithm until it reaches the valley of the global minimum. The proposed hybrid scheme can be applied to different inverse problems in NIRS which are solved iteratively. We find that for both numerical and phantom experiments optical properties such as the absorption and reduced scattering coefficients can be retrieved without being trapped by a local minimum when Monte Carlo simulation is run only about $100$ steps before switching to an iterative method. The hybrid method is compared with simulated annealing. Although the Metropolis-Hastings MCMC arrives at the steady state at about $10000$ Monte Carlo steps, in the hybrid method the Monte Carlo simulation can be stopped way before the burn-in time.
△ Less
Submitted 24 August, 2019; v1 submitted 2 February, 2018;
originally announced February 2018.
-
Transferred plasma jet from a dielectric barrier discharge for processing of poly(dimethylsiloxane) surfaces
Authors:
Fellype do Nascimento,
Munemasa Machida,
Mara A. Canesqui,
Stanislav Moshkalev
Abstract:
In this work we studied processing of poly(dimethylsiloxane) (PDMS) surfaces using dielectric barrier discharge (DBD) plasma in two different assemblies, one using the primary plasma jet obtained from a conventional DBD and the other using a DBD plasma jet transfer. The evolution of water contact angle (WCA) in function of plasma processing time and in function of aging time as well as the changes…
▽ More
In this work we studied processing of poly(dimethylsiloxane) (PDMS) surfaces using dielectric barrier discharge (DBD) plasma in two different assemblies, one using the primary plasma jet obtained from a conventional DBD and the other using a DBD plasma jet transfer. The evolution of water contact angle (WCA) in function of plasma processing time and in function of aging time as well as the changes in the surface roughness of PDMS samples for both plasma treatments have been studied. We also compared vibrational and rotational temperatures for both plasmas and for the first time the vibrational temperature (T_vib) for the transferred plasma jet has been shown to be higher as compared with the primary jet. The increment in the T_vib value seems to be the main reason for the improvements in adhesion properties and surface wettability for the transferred plasma jet. Possible explanations for the increase in the vibrational temperature are presented.
△ Less
Submitted 5 August, 2019; v1 submitted 24 June, 2016;
originally announced June 2016.
-
A screened automated structural search with semiempirical methods
Authors:
Yukihiro Ota,
Sergi Ruiz-Barragan,
Masahiko Machida,
Motoyuki Shiga
Abstract:
We developed an interface program between a program suite for an automated search of chemical reaction pathways, GRRM, and a program package of semiempirical methods, MOPAC. A two-step structural search is proposed as an application of this interface program. A screening test is first performed by semiempirical calculations. Subsequently, a reoptimization procedure is done by ab initio or density…
▽ More
We developed an interface program between a program suite for an automated search of chemical reaction pathways, GRRM, and a program package of semiempirical methods, MOPAC. A two-step structural search is proposed as an application of this interface program. A screening test is first performed by semiempirical calculations. Subsequently, a reoptimization procedure is done by ab initio or density functional calculations. We apply this approach to ion adsorption on cellulose. The computational efficiency is also shown for a GRRM search. The interface program is suitable for the structural search of large molecular systems for which semiempirical methods are applicable.
△ Less
Submitted 9 February, 2016;
originally announced February 2016.
-
Experimental analysis of DBD plasma jet properties using different gases and two kinds of transfer plate
Authors:
Fellype do Nascimento,
Stanislav Moshkalev,
Munemasa Machida
Abstract:
Dielectric Barrier Discharge (DBD) plasma jets has been studied extensively in recent years because of its wide range of applications. DBD plasmas can be produced using many different gases and can be applied to a broad variety of surfaces and substrates. In this work, we provide a comparison of DBD plasmas generated using argon (Ar), helium (He) and nitrogen (N2), as well as their mixtures with w…
▽ More
Dielectric Barrier Discharge (DBD) plasma jets has been studied extensively in recent years because of its wide range of applications. DBD plasmas can be produced using many different gases and can be applied to a broad variety of surfaces and substrates. In this work, we provide a comparison of DBD plasmas generated using argon (Ar), helium (He) and nitrogen (N2), as well as their mixtures with water vapor in order to know how some plasma properties are affected by the use of different gases. All plasmas were studied in two different conditions, using a transfer plate made of a conductive material and using a transfer plate made of an insulating one. We observed that the processes of excitation and ionization of nitrogen molecules by direct collisions with Ar or He are more evident and significant in He plasmas, which means that He atoms in metastable states have greater ability to transfer energy to molecules of nitrogen in the plasma. The collisions of He atoms in metastable states with N2 molecules determine the vibrational temperature (Tvib) values in He plasmas, while in Ar and N2 plasmas the Tvib values are determined mainly by collisions of electrons with N2 molecules.
△ Less
Submitted 5 August, 2019; v1 submitted 11 December, 2015;
originally announced December 2015.
-
Fields of View for Environmental Radioactivity
Authors:
Alex Malins,
Masahiko Okumura,
Masahiko Machida,
Hiroshi Takemiya,
Kimiaki Saito
Abstract:
The gamma component of air radiation dose rates is a function of the amount and spread of radioactive nuclides in the environment. These radionuclides can be natural or anthropogenic in origin. The field of view describes the area of radionuclides on, or below, the ground that is responsible for determining the air dose rate, and hence correspondingly the external radiation exposure. This work des…
▽ More
The gamma component of air radiation dose rates is a function of the amount and spread of radioactive nuclides in the environment. These radionuclides can be natural or anthropogenic in origin. The field of view describes the area of radionuclides on, or below, the ground that is responsible for determining the air dose rate, and hence correspondingly the external radiation exposure. This work describes Monte Carlo radiation transport calculations for the field of view under a variety of situations. Presented first are results for natural 40K and thorium and uranium series radionuclides distributed homogeneously within the ground. Results are then described for atmospheric radioactive caesium fallout, such as from the Fukushima Daiichi Nuclear Power Plant accident. Various stages of fallout evolution are considered through the depth distribution of 134Cs and 137Cs in soil. The fields of view for the natural radionuclides and radiocaesium are different. This can affect the responses of radiation monitors to these nuclides if the detector is partially shielded from the ground within its field of view. The field of view also sets the maximum reduction in air dose rates that can be achieved through local decontamination or remediation measures. This maximum efficiency can be determined quickly from the data presented here for the air dose rate versus the spatial extent of radioactive source on the ground.
△ Less
Submitted 2 November, 2015; v1 submitted 30 September, 2015;
originally announced September 2015.
-
Evaluation of ambient dose equivalent rates influenced by vertical and horizontal distribution of radioactive cesium in soil in Fukushima Prefecture
Authors:
Alex Malins,
Hiroshi Kurikami,
Shigeo Nakama,
Tatsuo Saito,
Masahiko Okumura,
Masahiko Machida,
Akihiro Kitamura
Abstract:
The air dose rate in an environment contaminated with 134Cs and 137Cs depends on the amount, depth profile and horizontal distribution of these contaminants within the ground. This paper introduces and verifies a tool that models these variables and calculates ambient dose equivalent rates at 1 m above the ground. Good correlation is found between predicted dose rates and dose rates measured with…
▽ More
The air dose rate in an environment contaminated with 134Cs and 137Cs depends on the amount, depth profile and horizontal distribution of these contaminants within the ground. This paper introduces and verifies a tool that models these variables and calculates ambient dose equivalent rates at 1 m above the ground. Good correlation is found between predicted dose rates and dose rates measured with survey meters in Fukushima Prefecture in areas contaminated with radiocesium from the Fukushima Dai-ichi Nuclear Power Plant accident. This finding is insensitive to the choice for modelling the activity depth distribution in the ground using activity measurements of collected soil layers, or by using exponential and hyperbolic secant fits to the measurement data. Better predictions are obtained by modelling the horizontal distribution of radioactive cesium across an area if multiple soil samples are available, as opposed to assuming a spatially homogeneous contamination distribution. Reductions seen in air dose rates above flat, undisturbed fields in Fukushima Prefecture are consistent with decrement by radioactive decay and downward migration of cesium into soil. Analysis of remediation strategies for farmland soils confirmed that topsoil removal and interchanging a topsoil layer with a subsoil layer result in similar reductions in the air dose rate. These two strategies are more effective than reverse tillage to invert and mix the topsoil.
△ Less
Submitted 14 September, 2015; v1 submitted 14 September, 2015;
originally announced September 2015.
-
Numerical algorithms of the radiative transport equation using rotated reference frames for optical tomography with structured illumination
Authors:
Manabu Machida
Abstract:
We consider optical tomography with structured illumination in spatial-frequency domain using the three-dimensional radiative transport equation. Without the diffusion approximation, the radiative transport equation is solved by the technique of rotated reference frames. In addition to the method of rotated reference frames (spherical-harmonic expansion), the three dimensional FN method is applied…
▽ More
We consider optical tomography with structured illumination in spatial-frequency domain using the three-dimensional radiative transport equation. Without the diffusion approximation, the radiative transport equation is solved by the technique of rotated reference frames. In addition to the method of rotated reference frames (spherical-harmonic expansion), the three dimensional FN method is applied to this optical tomography.
△ Less
Submitted 31 March, 2020; v1 submitted 26 July, 2015;
originally announced July 2015.
-
Polarization oscillations of near-field thermal emission
Authors:
Manabu Machida,
Evgenii Narimanov,
John C. Schotland
Abstract:
We consider the polarization of thermal emission in the near-field of various materials including dielectrics and metallic systems with resonant surface modes. We find that at thermal equilibrium, the degree of polarization exhibits spatial oscillations with a period of approximately half the optical wavelength, independent of material composition. This result contrasts with that of Setala, Kaivol…
▽ More
We consider the polarization of thermal emission in the near-field of various materials including dielectrics and metallic systems with resonant surface modes. We find that at thermal equilibrium, the degree of polarization exhibits spatial oscillations with a period of approximately half the optical wavelength, independent of material composition. This result contrasts with that of Setala, Kaivola and Friberg [Phys. Rev. Lett. 88, 123902 (2002)], who find monotonic decay of the degree of polarization for systems in local thermal equilibrium.
△ Less
Submitted 1 June, 2015;
originally announced June 2015.
-
Treatment of PDMS surfaces using pulsed DBD plasmas: comparing the use of different gases and its influence on adhesion
Authors:
Fellype do Nascimento,
Stanislav Moshkalev,
Munemasa Machida,
Sergio Parada
Abstract:
In this work we present some results of the treatment of polydimethylsiloxane (PDMS) surfaces using pulsed dielectric barrier discharge plasmas. The results of plasma treatment using different gases and mixtures, argon, argon plus water vapor, helium, helium plus water vapor, nitrogen and nitrogen plus water vapor, were compared testing the adhesion between two PDMS samples for each kind of plasma…
▽ More
In this work we present some results of the treatment of polydimethylsiloxane (PDMS) surfaces using pulsed dielectric barrier discharge plasmas. The results of plasma treatment using different gases and mixtures, argon, argon plus water vapor, helium, helium plus water vapor, nitrogen and nitrogen plus water vapor, were compared testing the adhesion between two PDMS samples for each kind of plasma. We also studied the water contact angle in function of plasma process time of PDMS surfaces with each kind of plasma treatment. The plasma was characterized by optical emission spectroscopy (OES) to identify the emitting species and determine the plasma temperatures. The plasma temperature for each process was estimated comparing the spectrum obtained by OES with the spectrum generated by SpecAir simulation code. Measurements of power delivered to the plasmas were also performed. As the results, all the process using different gases show good adhesion efficacy between PDMS samples when long exposure time (larger than 150 seconds) is applied. However, when only a few discharges are applied to PDMS samples the helium plasma process presented best results. Atomic Force Microscopy (AFM) analysis of PDMS samples treated with helium plasma showed reduction in the surface roughness, which increase the surface contact area and improves the adhesion.
△ Less
Submitted 19 May, 2015; v1 submitted 14 April, 2015;
originally announced April 2015.
-
Topographic Effects on Ambient Dose Equivalent Rates from Radiocesium Fallout
Authors:
Alex Malins,
Masahiko Okumura,
Masahiko Machida,
Kimiaki Saito
Abstract:
Land topography can affect air radiation dose rates by locating radiation sources closer to, or further from, detector locations when compared to perfectly flat terrain. Hills and slopes can also shield against the propagation of gamma rays. To understand the possible magnitude of topographic effects on air dose rates, this study presents calculations for ambient dose equivalent rates at a range o…
▽ More
Land topography can affect air radiation dose rates by locating radiation sources closer to, or further from, detector locations when compared to perfectly flat terrain. Hills and slopes can also shield against the propagation of gamma rays. To understand the possible magnitude of topographic effects on air dose rates, this study presents calculations for ambient dose equivalent rates at a range of heights above the ground for varying land topographies. The geometries considered were angled ground at the intersection of two planar surfaces, which is a model for slopes neighboring flat land, and a simple conical geometry, representing settings from hilltops to valley bottoms. In each case the radiation source was radioactive cesium fallout, and the slope angle was varied systematically to determine the effect of topography on the air dose rate. Under the assumption of homogeneous fallout across the land surface, and for these geometries and detector locations, the dose rates at high altitudes are more strongly affected by the underlying land topography than those close to ground level. At a height of 300m, uneven topographies can lead to a 50% change in air dose rates compared to if the ground were uniformly flat. However, in practice the effect will more often than not be smaller than this, and heterogeneity in the source distribution is likely to be a more significant factor in determining local air dose rates.
△ Less
Submitted 15 April, 2015; v1 submitted 13 February, 2015;
originally announced February 2015.
-
Comment on `Update of 40K and 226Ra and 232Th series $γ$-to-dose conversion factors for soil'
Authors:
Alex Malins,
Masahiko Machida,
Kimiaki Saito
Abstract:
A letter to the editor of the Journal of Environmental Radioactivity on the article: E. Gasser, A. Nachab, A. Nourreddine, Ch. Roy, and A. Sellam, `Update of 40K and 226Ra and 232Th series $γ$-to-dose conversion factors for soil', J. Environ. Radioactiv. 138, 68-71 (2014), DOI: 10.1016/j.jenvrad.2014.08.002.
A letter to the editor of the Journal of Environmental Radioactivity on the article: E. Gasser, A. Nachab, A. Nourreddine, Ch. Roy, and A. Sellam, `Update of 40K and 226Ra and 232Th series $γ$-to-dose conversion factors for soil', J. Environ. Radioactiv. 138, 68-71 (2014), DOI: 10.1016/j.jenvrad.2014.08.002.
△ Less
Submitted 19 January, 2015; v1 submitted 19 January, 2015;
originally announced January 2015.
-
Generation of Cold Argon Plasma Jet at the End of Flexible Plastic Tube
Authors:
Konstantin G. Kostov,
Munemasa Machida,
Vadym Prysiazhnyi
Abstract:
This brief communication reports a new method for generation of cold atmospheric pressure plasma jet at the downstream end of a flexible plastic tube. The device consists of a small chamber where dielectric barrier discharge (DBD) is ignited in Argon. The discharge is driven by a conventional low frequency AC power supply. The exit of DBD reactor is connected to a commercial flexible plastic tube…
▽ More
This brief communication reports a new method for generation of cold atmospheric pressure plasma jet at the downstream end of a flexible plastic tube. The device consists of a small chamber where dielectric barrier discharge (DBD) is ignited in Argon. The discharge is driven by a conventional low frequency AC power supply. The exit of DBD reactor is connected to a commercial flexible plastic tube (up to 4 meters long) with a thin floating Cu wire inside. Under certain conditions an Ar plasma jet can be extracted from the downstream tube end and there is no discharge inside the plastic tube. The jet obtained by this method is cold enough to be put in direct contact with human skin without electric shock and can be used for medical treatment and decontamination.
△ Less
Submitted 14 November, 2014;
originally announced November 2014.
-
Chandrasekhar polynomials -- A brief review
Authors:
Manabu Machida
Abstract:
A review on the Chandrasekhar polynomials is given. The polynomials often appear in transport theory. The relation to the method of rotated reference frames for the three-dimensional radiative transport equation is clarified.
A review on the Chandrasekhar polynomials is given. The polynomials often appear in transport theory. The relation to the method of rotated reference frames for the three-dimensional radiative transport equation is clarified.
△ Less
Submitted 24 May, 2020; v1 submitted 12 April, 2013;
originally announced April 2013.
-
Singular eigenfunctions for the three-dimensional radiative transport equation
Authors:
Manabu Machida
Abstract:
Case's method obtains solutions to the radiative transport equation as superpositions of elementary solutions when the specific intensity depends on one spatial variable. In this paper, we find elementary solutions when the specific intensity depends on three spatial variables in three-dimensional space. By using the reference frame whose z-axis lies in the direction of the wave vector, the angula…
▽ More
Case's method obtains solutions to the radiative transport equation as superpositions of elementary solutions when the specific intensity depends on one spatial variable. In this paper, we find elementary solutions when the specific intensity depends on three spatial variables in three-dimensional space. By using the reference frame whose z-axis lies in the direction of the wave vector, the angular part of each elementary solution becomes the singular eigenfunction for the one-dimensional radiative transport equation. Thus Case's method is generalized.
△ Less
Submitted 7 September, 2013; v1 submitted 6 April, 2013;
originally announced April 2013.
-
Transport of Entanglement
Authors:
Manabu Machida,
Vadim A Markel,
John C Schotland
Abstract:
We consider the propagation of two-photon light in a random medium. We show that the Wigner distribution of the two-photon wave function obeys an equation that is analogous to the radiative transport equation for classical light. Using this result, we predict that the entanglement of a photon pair is destroyed with propagation.
We consider the propagation of two-photon light in a random medium. We show that the Wigner distribution of the two-photon wave function obeys an equation that is analogous to the radiative transport equation for classical light. Using this result, we predict that the entanglement of a photon pair is destroyed with propagation.
△ Less
Submitted 17 October, 2012;
originally announced October 2012.
-
Energy Spectra of Quantum Turbulence: Large-scale Simulation and Modeling
Authors:
Narimasa Sasa,
Takuma Kano,
Masahiko Machida,
Victor S. L'vov,
Oleksii Rudenko,
Makoto Tsubota
Abstract:
In $2048^3$ simulation of quantum turbulence within the Gross-Pitaevskii equation we demonstrate that the large scale motions have a classical Kolmogorov-1941 energy spectrum E(k) ~ k^{-5/3}, followed by an energy accumulation with E(k) ~ const at k about the reciprocal mean intervortex distance. This behavior was predicted by the L'vov-Nazarenko-Rudenko bottleneck model of gradual eddy-wave cross…
▽ More
In $2048^3$ simulation of quantum turbulence within the Gross-Pitaevskii equation we demonstrate that the large scale motions have a classical Kolmogorov-1941 energy spectrum E(k) ~ k^{-5/3}, followed by an energy accumulation with E(k) ~ const at k about the reciprocal mean intervortex distance. This behavior was predicted by the L'vov-Nazarenko-Rudenko bottleneck model of gradual eddy-wave crossover [J. Low Temp. Phys. 153, 140-161 (2008)], further developed in the paper.
△ Less
Submitted 13 June, 2011; v1 submitted 18 August, 2010;
originally announced August 2010.