Cross-shaped metal–semiconductor–metal plasmonic crystal for terahertz modulator

The transmission and tuning properties of a cross-shaped plasmonic crystal based on periodic metal–semiconductor–metal (MSM) structures have been investigated in the terahertz (THz) regime. According to the mode analysis, we find that the different resonance modes in the plasmonic crystal show the different changes when this device is actively controlled by the carrier injection of the MSM structures. The longitudinal modes disappear, while the horizontal mode moves to a higher frequency. The former leads to an intensity modulation at 0.5 THz and 1.1 THz when the groove depth h = 60 μm, and the later leads to a band blue-shift from 1.325 THz to 1.38 THz. These results will be applied to THz modulation and tunable filtering.     

Fast directional multilevel summation for oscillatory kernels based on Chebyshev interpolation

Many applications lead to large systems of linear equations with dense matrices. Direct matrix-vector products become prohibitive, since the computational cost increases quadratically with the size of the problem. By exploiting specific kernel properties fast algorithms can be constructed.A directional multilevel algorithm for translation-invariant oscillatory kernels of the type K(x, y) = G(x ? y)e^?k∣x?y∣, with G(x ? y) being any smooth kernel, will be presented. We will first present a general approach to build fast multipole methods (FMMs) based on Chebyshev interpolation and the adaptive cross approximation (ACA) for smooth kernels. The Chebyshev interpolation is used to transfer information up and down the levels of the FMM. The scheme is further accelerated by compressing the information stored at Chebyshev interpolation points using ACA and QR decompositions. This leads to a nearly optimal computational cost with a small pre-processing time due to the low computational cost of ACA. This approach is in particular faster than performing singular value decompositions.This does not address the difficulties associated with the oscillatory nature of K. For that purpose, we consider the following modification of the kernel K^u = K(x, y)e^??ku·(x?y), where u is a unit vector (see Brandt [1]). We proved that the kernel K^u can be interpolated efficiently when x ? y lies in a cone of direction u. This result is used to construct an FMM for the kernel K.Theoretical error bounds will be presented to control the error in the computation as well as the computational cost of the method. The paper ends with the presentation of 2D and 3D numerical convergence studies, and computational cost benchmarks.     

Aqueous ionic liquid based ultrasonic assisted extraction of four acetophenones from the Chinese medicinal plant Cynanchum bungei Decne

In this study, an aqueous ionic liquid based ultrasonic assisted extraction (ILUAE) method for the extraction of the four acetophenones, namely 4-hydroxyacetophenone (1), 2,5-dihydroxyacetophenone (2), baishouwubenzophenone (3) and 2,4-dihydroxyacetophenone (4) from the Chinese medicinal plant Cynanchum bungei was developed. Three kinds of aqueous l-alkyl-3-methylimidazolium ionic liquids with different anion and alkyl chain were investigated. The results indicated that ionic liquids (ILs) showed remarkable effects on the extraction efficiency of acetophenones. In addition, the ILUAE, including several ultrasonic parameters, such as the ILs concentration, solvent to solid ratio, power, particle size, temperature, and extraction time have been optimized. Under these optimal conditions (e.g., with 0.6 M [C₄MIM]BF₄, solvent to solid ratio of 35:1, power of 175 W, particle size of 60–80 mesh, temperature of 25 °C and time of 50 min), this approach gained the highest extraction yields of four acetophenones 286.15, 21.65, 632.58 and 205.38 μg/g, respectively. The proposed approach has been evaluated by comparison with the conventional heat-reflux extraction (HRE) and regular UAE. The results indicated that ILUAE is an alternative method for extracting acetophenones from C. bungei.     

Design and VNA-measurement of coplanar waveguide (CPW) on benzocyclobutene (BCB) at THz frequencies

The low permittivity and the low loss tangent of the benzocyclobutene polymer (BCB) offers to coplanar waveguides (CPW) a low dispersive propagation properties at THz frequency. These transmission lines have been designed, modeled with a three dimensional (3D) solver of Maxwell equations based on finite element method (FEM) from 20 to 1000 GHz at various characteristic impedances (Z_c). Their dispersion and losses (radiation, conduction and dielectric) have been investigated separately versus the waveguide size, the nature of the substrate (dielectric or semiconductor) to optimize the THz signal propagation. Monomode CPW on BCB numerically designed for various Z_c were realized and measured with vector network analyzer (VNA). S-parameters of CPW are de-embedded by optimization of the accesses’ model. A good agreement is found between experimental and numerical results with low attenuation constants of 2.7 dB/mm and 3.5 dB/mm at 400 GHz and 500 GHz, respectively.     

A waveguide high-pass filter system for measuring the spectrum of pulsed terahertz sources

We propose a system for measuring spectra of terahertz (THz) pulses, including single pulses, which is based on high-pass filters (HPFs). The system consists of channels for measuring amplitudes of pulses (initial pulses and those transmitted via HPFs with different cutoff frequencies) and an algorithm for processing of the obtained data. The pulse spectrum is restored by using the iteration method or the amplitude–frequency method. The iteration method of spectrum restoration is applicable in the range of THz pulse durations from 10⁻⁹ s to 10⁻⁷ s. The amplitude–frequency method is applicable to THz pulses with durations exceeding 10⁻⁸ s. The system for measuring of THz pulse spectra was simulated by using the characteristics of specially developed waveguide HPFs. The relative simulation error of determining the central frequency by the amplitude–frequency method is equal to 2 · 10⁻⁶ for THz pulse durations of 10⁻⁵ s and longer.     

Ion conduction in proton- and related defect (super) ionic conductors: Mechanical,electronic and structure parameters

We find the relationships among optical dielectric constant ε_∞, activation energy E_ac, averaged atomic mass per a formula unit, ∑_jm_j / N, volume V and transition temperature T_c for various type ion conductors with forms of E_ac = α / (ε_∞ ? β), E_ac = A₀ + δ / [(∑_jm_j / N) ? d], E_ac = A_v / V^2/3, and ln(T_c) = g ? hln(∑_jm_j), where α, β, δ, A₀, d, A_v, g and h are constants depending on the kinds of conduction elements. We derive those proportional forms from a simple equation of motion under the assumption of ion hopping assisted by enhanced vibration displacement of host lattice. The enhancement is induced from the large fourth-order term of the host lattice potential originating from the electronic shielding effect of Coulomb force, heavy atomic mass of constituent ion, and volume expansion under the long-range periodicity of crystal structure. This mechanism is ascertained from characteristic phenomena of various kinds of conduction elements. For impurity-type H⁺-ion or defect conductor, the proportional form is shifted from that of superionic conductor because of weakened effect of host lattice vibration mode on H⁺-ion or O-ion defect. Photo-induced spectra of mobile ion in AgCl are understood, and a small quantum effect of H⁺ -ion is suggested.     

Characterization of Young’s modulus of silicon versus temperature using a “beam deflection” method with a four-point bending fixture

Young’s modulus (E) and Poisson’s ratio (ν) are dependent upon the direction on the silicon surface. In this work, E and ν of silicon have been calculated analytically for any crystallographic direction of silicon by using compliance coefficients (s₁₁, s₁₂, and s₄₄), and the values of E are confirmed experimentally by using a “beam deflection” method with a four-point bending fixture. Experimental results for E as a function of temperature from −150 °C to +150 °C are presented for (0 0 1) and (1 1 1) silicon wafers.     

High resolution emission spectroscopy of the E2Π–X2Σ+ transition of SrH and SrD

Emission spectra of SrH and SrD have been studied at high resolution using a Fourier transform spectrometer. The molecules have been produced in a high temperature furnace from the reaction of strontium metal vapor with H₂/D₂ in the presence of a slow flow of Ar gas. The spectra observed in the 18 000–19 500 cm^?1 region consist of the 0–0 and 1–1 bands of the E²Π–X²Σ⁺ transition of the two isotopologues. A rotational analysis of these bands has been obtained by combining the present measurements with previously available pure rotation and vibration–rotation measurements for the ground state, and improved spectroscopic constants have been obtained for the E²Π state. The present analysis provides spectroscopic constants for the E²Π state as ΔG(½) = 1166.1011(15) cm^?1, B_e = 3.805503(32) cm^?1, α_e = 0.098880(47) cm^?1, r_e = 2.1083727(89) Å for SrH, and ΔG(½) = 839.1283(23) cm^?1, B_e = 1.918564(15) cm^?1, α_e = 0.034719(23) cm^?1, r_e = 2.1121943(83) Å for SrD.     

Resonances in the final Rydberg state population of multiply charged ions ArVIII,KrVIII, and XeVIII escaping solid surfaces

The appearance of resonances (pronounced maxima at n_A = n_res) in the probability distributions for the population of the Rydberg state (n_A, l_A, m_A) of multiply charged ions (Z ? 1) escaping solid surfaces at intermediate velocities (v ≈ 1 a.u.) is discussed. Within the framework of the time-symmetrized two-state vector model, in which the state of a single active electron is described by two wave functions Ψ₁ and Ψ₂, the resonances are explained by means of an electron tunneling in the very vicinity of the ion–surface potential barrier top. To include this specific feature of electron transitions into the model, the appropriate etalon equation method is used in the calculation of the function Ψ₁. We consider the ions ArVIII, KrVIII, and XeVIII with the same core charges Z = 8 a.u., but with different core polarizations. The effect of the ionic core polarization is associated with the function Ψ₂. The population probabilities for n_A ≈ n_res are complemental to those obtained recently for n_A < n_res, and in sufficiently good agreement with available beam-foil experimental data. The pronounced resonances in the final population distributions are recognized only in the case of ArVIII ion and for the lower values of the solid work function (argon anomaly).     

Cytogenetics of four Omophoita species (Coleoptera,Chrysomelidae, Alticinae): A comparative analysis using mitotic and meiotic cells submitted to the standard staining and C-banding technique

Omophoita belongs to the tribe Oedionychini and is endemic from Neotropical region. The species of the tribe Oedionychini have revealed certain singular chromosomal features, such as sex chromosomes with extremely large size, asynapsis, and synthelic or amphithelic orientation during meiosis. Additionally, some species also showed post-reductional segregation of the gigantic sex chromosomes in meiotic division. The purpose of this work was to characterize cytogenetically four Omophoita species (O. magniguttis, O. octoguttata, O. personata, and O. sexnotata) in relation to their diploid number, chromosomal morphology, type of sex chromosome system, and constitutive heterochromatin pattern in mitotic and meiotic cells, and compare the obtained data with those of related species to establish the mechanism involved in the chromosomal differentiation of these species during the evolutionary process. The diploid number, 2n = 22 = 20 + X + y, and meiotic formulae, 10II + X + y, observed in these species were similar to those of the same genus and other species related. The autosomal morphology was acrocentric in O. magniguttis and O. octoguttata, metacentric in O. personata, and predominantly metacentric in O. sexnotata. In all these species, the sex chromosomes were metacentric. The secondary constriction occurred in pair 6 and X chromosome of O. personata, and in pair 6 and y chromosome of O. sexnotata. The constitutive heterochromatin was pericentromeric in O. magniguttis and centromeric in O. sexnotata, with the exception of the mitotic sex chromosomes of O. sexnotata, in which centromeric C band was lacking. Additional C bands in the sex chromosomes of O. magniguttis and certain autosomes and sex chromosomes of O. sexnotata were observed. Collochores were indirectly identified in the spermatocytes of O. octoguttata, O. personata, and O. sexnotata. The main mechanisms involved in the karyotype evolution of these species were discussed.     

Experimental study of minimum ignition energy of lean H2-N2O mixtures

Ignition energies for short duration (<50 ns) spark discharges were measured for undiluted and nitrogen-diluted H₂-N₂O mixtures of equivalence ratios ? = 0.15 and 0.2, dilution of 0% and 20% N₂, and initial pressures of 15–25 kPa. The ignition events were analyzed using statistical tools and the probability of ignition versus spark energy density (spark energy divided by the spark length) was obtained. The simple cylindrical ignition kernel model was compared against the results from the present study. Initial pressure has a significant effect on the width of the probability distribution, ranging from a broad (P = 15 kPa) to a narrow (P = 25 kPa) probability distribution indicating that the statistical variation of median spark energy density increases as initial pressure of the mixture decreases. A change in the equivalence ratio from 0.15 to 0.2 had a small effect on the median spark energy density. The addition of 20% N₂ dilution caused a significant increase in the median spark energy density when compared to no dilution. The extrapolation of the present results to atmospheric pressure, stoichiometric H₂-N₂O indicates that the electrostatic discharge ignition hazards are comparable to or greater than H₂-O₂ mixtures.     

Terahertz time-domain spectroscopy and optical properties of AM2O8 (A=Zr,Hf and M=W,Mo)

Terahertz time domain spectroscopy and optical properties of cubic ZrW₂O₈ and HfW₂O₈, and trigonal ZrMo₂O₈ and HfMo₂O₈ have been studied for the first time. It is shown that α-ZrW₂O₈ and α-HfW₂O₈ exhibit absorption peaks around 1.23 THz (41 cm^?1) and 1.18 THz (39.5 cm^?1), respectively, whereas no corresponding absorptions are observed for trigonal ZrMo₂O₈ and HfMo₂O₈. These external modes are suggested to originate from the librational and/or translational motions of the linked ZrO₆/HfO₆ and WO₄ polyhedra and to contribute to the negative thermal expansion of α-ZrW₂O₈ and α-HfW₂O₈. The absorption coefficients and refractive indices of these materials are derived from the transmission amplitudes of the THz field passing through the samples. The results suggest that cubic ZrW₂O₈ and HfW₂O₈ have the potential for THz filters or attenuators.     

Synthesis,growth and characterization of a semiorganic nonlinear optical single crystal: L-Valine cadmium bromide

A semiorganic nonlinear optical material L-valine cadmium bromide was synthesized and single crystal grown from an aqueous solution by the method of slow evaporation technique at room temperature. The grown crystal was characterized by Powder X-ray diffraction and FT IR studies. The LVCB crystallizes in monoclinic system with cell parameters a = 10.144(2) Å, b = 5.54(1) Å, c = 12.07(2) Å, β = 109.115(2)° with space group P₂₁. Thermal behavior and stability of crystal were studied using thermogravimetric analyses (TGA) and differential thermal analysis (DTA) techniques. The suitability of this material for NLO application was studied by optical absorption studies and second harmonic generation (SHG) efficiency measurement by Kurtz–Perry powder method.     

First high-resolution analysis of the ν15, ν12, ν5, ν10 and ν2 bands of oxirane

Fourier transform spectra of oxirane (ethylene oxide, c-C₂H₄O) have been recorded in the 730–1560 cm^?1 (6.4–13.7 μm) spectral region using a Bruker IFS125HR spectrometer at a resolution of 0.0019 cm^?1. A total of six vibration bands, ν₁₅, ν₁₂, ν₅, ν₃, ν₁₀ and ν₂, have been observed and analyzed. The corresponding upper state ro-vibrational levels were fit using Hamiltonian matrices accounting for various interactions. Satisfactory fits were obtained using the following polyads {15¹, 12¹, 5¹} and {10¹, 2¹} of interacting states. As a result, an accurate and extended set of Hamiltonian constants were obtained. The following band centers were derived: ν₀ (ν₁₅) = 808.13518(60) cm^?1, ν₀ (ν₁₂) = 822.27955(37) cm^?1, ν₀ (ν₅) = 876.72592(15), ν₀ (ν₃) = 1270.37032(10) cm^?1, ν₀ (ν₁₀) = 1471.35580(50) cm^?1 and ν₀ (ν₂) = 1497.83309(15) cm^?1 where the uncertainties are one standard deviation.     

Synthesis,crystal growth and characterization of a semiorganic material: Calcium dibromide bis(glycine) tetrahydrate

Single crystals of semiorganic material calcium dibromide bis(glycine) tetrahydrate were grown from aqueous solution. The crystal belongs to monoclinic system, with a = 13.261(5) Å, b = 6.792(2) Å, c = 15.671(9) Å and β = 91.68(4)°. The presence of the elements in the title compound was confirmed by energy dispersive X-ray analysis. The solubility and metastable zone width were found. The grown crystals were tested by powder XRD, FTIR, Thermo Gravimetric and Differential Thermal Analysis, UV–vis–NIR analysis, dielectrical and mechanical studies. The transmittance of calcium dibromide bis(glycine) tetrahydrate crystal has been used to calculate the refractive index n, the extinction coefficient K and both the real ɛ_r and imaginary ɛ_i components of the dielectric constant as functions of wavelength. The optical band gap of calcium dibromide bis(glycine) tetrahydrate is 3.23 eV.     

Numerical and perturbative computations of solitary waves of the Benjamin–Ono equation with higher order nonlinearity using Christov rational basis functions

Computation of solitons of the cubically-nonlinear Benjamin–Ono equation is challenging. First, the equation contains the Hilbert transform, a nonlocal integral operator. Second, its solitary waves decay only as O(1/∣x∣²). To solve the integro-differential equation for waves traveling at a phase speed c, we introduced the artificial homotopy H(u_XX) ? c u + (1 ? δ)u² + δu³ = 0, δ ∈ [0, 1] and solved it in two ways. The first was continuation in the homotopy parameter δ, marching from the known Benjamin–Ono soliton for δ = 0 to the cubically-nonlinear soliton at δ = 1. The second strategy was to bypass continuation by numerically computing perturbation series in δ and forming Padé approximants to obtain a very accurate approximation at δ = 1. To further minimize computations, we derived an elementary theorem to reduce the two-parameter soliton family to a parameter-free function, the soliton symmetric about the origin with unit phase speed. Solitons for higher order Benjamin–Ono equations are also computed and compared to their Korteweg–deVries counterparts. All computations applied the pseudospectral method with a basis of rational orthogonal functions invented by Christov, which are eigenfunctions of the Hilbert transform.     

Diffusion-promoted-desorption mechanism for D2 desorption from Si(100) surfaces

Temperature-programmed-desorption (TPD) spectra and isothermal desorption rates of D₂ molecules from a Si(100) surface have been calculated to reproduce experimental β_1, A-TPD spectra and isothermal desorption rate curves. In the diffusion-promoted-desorption (DPD) mechanism, hydrogen desorption from the Si(100) (2 × 1) surfaces takes place via D atom diffusion from doubly-occupied Si dimers (DODs) to their adjacent unoccupied Si dimers (UODs). Taking a clustering interaction among DODs into consideration, coverages θ_DU of desorption sites consisting of a pair of a DOD and UOD are evaluated by a Monte Carlo (MC) method. The TPD spectra for the β_1, A peak are obtained by numerically integrating the desorption rate equation R = ν_A exp(? E_d, A / k_BT)θ_DU, where ν_A is the pre-exponential factor and E_d, A is the desorption barrier. The TPD spectra calculated for E_d, A = 1. 6 eV and ν_A = 2.7 × 10⁹ /s are found to be in good agreement with the experimental TPD data for a wide coverage range from 0.01 to 0.74 ML. Namely, the deviation from first-order kinetics observed in the coverage dependent TPD spectra as well as in the isothermal desorption rate curves can be reproduced by the model simulations. This success in reproducing both the experimental TPD data and the very low desorption barrier validates the proposed DPD mechanism.     

Microfluidic acoustic trapping force and stiffness measurement using viscous drag effect

It has recently been demonstrated that it was possible to individually trap 70 μm droplets flowing within a 500 μm wide microfluidic channel by a 24 MHz single element piezo-composite focused transducer. In order to further develop this non-invasive approach as a microfluidic particle manipulation tool of high precision, the trapping force needs to be calibrated to a known force, i.e., viscous drag force arising from the fluid flow in the channel. However, few calibration studies based on fluid viscosity have been carried out with focused acoustic beams for moving objects in microfluidic environments.In this paper, the acoustic trapping force (F_trapping) and the trap stiffness (or compliance k) are experimentally determined for a streaming droplet in a microfluidic channel. F_trapping is calibrated to viscous drag force produced from syringe pumps. Chebyshev-windowed chirp coded excitation sequences sweeping the frequency range from 18 MHz to 30 MHz is utilized to drive the transducer, enabling the beam transmission through the channel/fluid interface for interrogating the droplets inside the channel. The minimum force (F_min,trapping) required for initially immobilizing drifting droplets is determined as a function of pulse repetition frequency (PRF), duty factor (DTF), and input voltage amplitude (V_in) to the transducer. At PRF = 0.1 kHz and DTF = 30%, F_min,trapping is increased from 2.2 nN for V_in = 22 V_pp to 3.8 nN for V_in = 54 V_pp. With a fixed V_in = 54 V_pp and DTF = 30%, F_min,trapping can be varied from 3.8 nN at PRF = 0.1 kHz to 6.7 nN at PRF = 0.5 kHz. These findings indicate that both higher driving voltage and more frequent beam transmission yield stronger traps for holding droplets in motion.The stiffness k can be estimated through linear regression by measuring the trapping force (F_trapping) corresponding to the displacement (x) of a droplet from the trap center. By plotting F_trapping – x curves for certain values of V_in (22/38/54 V_pp) at DTF = 10% and PRF = 0.1 kHz, k is measured to be 0.09, 0.14, and 0.20 nN/μm, respectively. With variable PRF from 0.1 to 0.5 kHz at V_in = 54 V_pp, k is increased from 0.20 to 0.42 nN/μm. It is shown that a higher PRF leads to a more compliant trap formation (or a stronger F_trapping) for a given displacement x. Hence the results suggest that this acoustic trapping method has the potential as a noninvasive manipulation tool for individual moving targets in microfluidics by adjusting the transducer’s excitation parameters.     

Volumetric properties and viscosities of the 2-pyrrolidone + 1,2-propanediol + water ternary system and its binary constituents at T = 313.15 K

Densities and viscosities of ternary mixtures of 2-pyrrolidone + 1,2-propanediol + water and corresponding binary mixtures of 1,2-propanediol + water, 2-pyrrolidone + water and 2-pyrrolidone + 1,2-propanediol have been measured over the whole composition range at 313.15 K. From the obtained data, the excess molar volumes (V^E), the deviations in viscosity (Δη) and the excess Gibbs free energy of activation (ΔG^?E) have been calculated. The V^E, Δη and ΔG^?E results were correlated and fitted by the Redlich–Kister equation for binary mixtures and by the Cibulka equation for ternary mixtures, as a function of mole fraction. Several predictive empirical relations were applied to predict the excess molar volumes of ternary mixtures from the binary mixing data.     

Interfacial reactions of Ba2YCu3O6+z with coated conductor buffer layer,LaMnO3

Chemical interactions between the Ba₂YCu₃O_6+x superconductor and the LaMnO₃ buffer layers employed in coated conductors have been investigated experimentally by determining the phases formed in the Ba₂YCu₃O_6+x–LaMnO₃ system. The Ba₂YCu₃O_6+x–LaMnO₃ join within the BaO–(Y₂O₃–La₂O₃)–MnO₂–CuO_x multi-component system is non-binary. At 810 °C (p_O2 = 100 Pa) and at 950 °C in purified air, four phases are consistently present along the join, namely, Ba_2?x(La_1+x?yY_y)Cu₃O_6+z, Ba(Y_2?xLa_x)CuO₅, (La_1?xY_x)MnO₃, (La,Y)Mn₂O₅. The crystal chemistry and crystallography of Ba(Y_2?xLa_x)CuO₅ and (La_1?xY_x)Mn₂O₅ were studied using the X-ray Rietveld refinement technique. The Y-rich and La-rich solid solution limits for Ba(Y_2?xLa_x)CuO₅ are Ba(Y_1.8La_0.2)CuO₅ and Ba(Y_0.1La_1.9)CuO₅, respectively. The structure of Ba(Y_1.8La_0.2)CuO₅ is Pnma (No. 62), a = 12.2161(5) Å, b = 5.6690(2) Å, c = 7.1468(3) Å, V = 494.94(4) Å³, and D_x = 6.29 g cm^?3. YMn₂O₅ and LaMn₂O₅ do not form solid solution at 810 °C (p_O2 = 100 Pa) or at 950 °C (in air). The structure of YMn₂O₅ was confirmed to be Pbam (No. 55), a = 7.27832(14) Å, b = 8.46707(14) Å, c = 5.66495(10) Å, and V = 349.108(14) Å³. A reference X-ray pattern was prepared for YMn₂O₅.