This year, Liu textit{et al}. [Phys. Rev. B textbf{104}, L041405 (2021)] proposed a new class of topological phonons (TPs; i.e., one-nodal surface (NS) phonons), which provides an effective route for realizing one-NSs in phonon systems. In this work,
based on first-principles calculations and symmetry analysis, we extended the types of NS phonons from one- to three-NS phonons. The existence of three-NS phonons (with NS states on the $k_{i}$ = $pi$ ($i$ = $x$, $y$, $z$) planes in the three-dimensional Brillouin zone (BZ)) is enforced by the combination of two-fold screw symmetry and time reversal symmetry. We screened all 230 space groups (SGs) and found nine candidate groups (with the SG numbers (Nos.) 19, 61, 62, 92, 96, 198, 205, 212, and 213) hosting three-NS phonons. Interestingly, with the help of first-principles calculations, we identified $P2_{1}$2$_{1}$2$_{1}$-type YCuS$_{2}$ (SG No. 19), $Pbca$-type NiAs$_{2}$ (SG No. 61), $Pnma$-type SrZrO$_{2}$ (SG No. 62), $P4_{1}$2$_{1}$2-type LiAlO$_{2}$ (SG No. 92), $P4_{3}$2$_{1}$2-type ZnP$_{2}$ (SG No. 96), $P2_{1}$3-type NiSbSe (SG No. 198), $Pabar{3}$-type As$_{2}$Pt (SG No. 205), $P4_{3}$32-type BaSi$_{2}$ (SG No. 212), and $P4_{1}$32-type CsBe$_{2}$F$_{5}$ (SG No. 213) as realistic materials hosting three-NS phonons. The results of our presented study enrich the class of NS states in phonon systems and provide concrete guidance for searching for three-NS phonons and singular Weyl point phonons in realistic materials.
Quantum spin systems may offer the first opportunities for beyond-classical quantum computations of scientific interest. While general quantum simulation algorithms likely require error-corrected qubits, there may be applications of scientific intere
st prior to the practical implementation of quantum error correction. The variational quantum eigensolver (VQE) is a promising approach to find energy eigenvalues on noisy quantum computers. Lattice models are of broad interest for use on near-term quantum hardware due to the sparsity of the number of Hamiltonian terms and the possibility of matching the lattice geometry to the hardware geometry. Here, we consider the Kitaev spin model on a hardware-native square-octagon qubit connectivity map, and examine the possibility of efficiently probing its rich phase diagram with VQE approaches. By benchmarking different choices of variational ansatz states and classical optimizers, we illustrate the advantage of a mixed optimization approach using the Hamiltonian variational ansatz (HVA). We further demonstrate the implementation of an HVA circuit on Rigettis Aspen-9 chip with error mitigation.
We discuss the implementation of quantum algorithms for lattice $Phi^4$ theory on circuit quantum electrodynamics (cQED) system. The field is represented on qudits in a discretized field amplitude basis. The main advantage of qudit systems is that it
s multi-level characteristic allows the field interaction to be implemented only with diagonal single-qudit gates. Considering the set of universal gates formed by the single-qudit phase gate and the displacement gate, we address initial state preparation and single-qudit gate synthesis with variational methods.
Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods to accurately predict error scaling
as a function of the resolution and parameters of the model that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the Finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wavefunction sampling and the finite representation size. Numerical simulations of small size $Phi^4$ problems demonstrate that the boson mass optimizing the sampling of the ground state wavefunction is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wavefunctions does not necessarily result in accurate representation. We develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.
In this ISSI-supported series of studies on magnetic helicity in the Sun, we systematically implement different magnetic helicity calculation methods on high-quality solar magnetogram observations. We apply finite-volume, discrete flux tube (in parti
cular, connectivity-based) and flux-integration methods to data from Hinodes Solar Optical Telescope. The target is NOAA active region 10930 during a ~1.5 day interval in December 2006 that included a major eruptive flare (SOL2006-12-13T02:14X3.4). Finite-volume and connectivity-based methods yield instantaneous budgets of the coronal magnetic helicity, while the flux-integration methods allow an estimate of the accumulated helicity injected through the photosphere. The objectives of our work are twofold: A cross-validation of methods, as well as an interpretation of the complex events leading to the eruption. To the first objective, we find (i) strong agreement among the finite-volume methods, (ii) a moderate agreement between the connectivity-based and finite-volume methods, (iii) an excellent agreement between the flux-integration methods, and (iv) an overall agreement between finite-volume and flux-integration based estimates regarding the predominant sign and magnitude of the helicity. To the second objective, we are confident that the photospheric helicity flux significantly contributed to the coronal helicity budget, and that a right-handed structure erupted from a predominantly left-handed corona during the X-class flare. Overall, we find that the use of different methods to estimate the (accumulated) coronal helicity may be necessary in order to draw a complete picture of an active-region corona, given the careful handling of identified data (preparation) issues, which otherwise would mislead the event analysis and interpretation.
Porosity variations in tubular scaffolds are critical to reproducible, sophisticated applications of electrospun fibers in biomedicine. Established laser micrometry techniques produced ~14,000 datapoints enabling thickness and porosity plots versus b
oth the azimuthal (Phi) and axial (Z) directions following cylindrical mandrel deposition. These 3D datasets could then be unrolled into maps revealing variations in thickness and porosity versus 0, -5, and -15 kV collector bias. As bias increases, thinner, more focused depositions occur. Simultaneous decreases in net porosity versus bias (91.1% 0kV > 83.4% -5kV > 80.2% -15 kV) are sensible, but significant changes in the distribution were unexpected. Surprisingly, at 0 kV, extensive mesoscale surface roughness is evident. Optical profilometry revealed unique features ~1600-420 mum in size, standing ~210 mum above the surrounding surface. These shrink to only ~440-150 mum in size and ~30 mum higher at -5 kV bias and disappear entirely at -15 kV. Scanning electron microscopy (SEM) resolved these into novel, localized domains containing tightly aligned fibers oriented parallel to the mandrel axis. Unexpectedly, we also observed substantial orientation along the mandrel axis. By modifying classical bending instability models to incorporate cylindrical electric fields, simulation revealed that horizontal components in the modified electric field alter bending loop shape, causing the observed alignment. This provides a new, easily utilized tool enabling facile, efficient tuning of orientation.
Generating and manipulating Dirac points in artificial atomic crystals has received attention especially in photonic systems due to their ease of implementation. In this paper, we propose a two-dimensional photonic crystal made of a Kekule lattice of
pure dielectrics, where the internal rotation of cylindrical pillars induces optical Dirac-degeneracy breaking. Our calculated dispersion reveals that the synchronized rotation reverses bands and switches parity as well so as to induce a topological phase transition. Our simulation demonstrates that such topologically protected edge states can achieve robust transmission in defect waveguides under deformation, and therefore provides a pragmatically tunable scheme to achieve reconfigurable topological phases.
Using a data sample of 980~fb$^{-1}$ collected with the Belle detector operating at the KEKB asymmetric-energy $e^+e^-$ collider, we present evidence for the $Omega(2012)^-$ in the resonant substructure of $Omega_{c}^{0} to pi^+ (bar{K}Xi)^{-}$ ($(ba
r{K}Xi)^{-}$ = $K^-Xi^0$ + $bar{K}^0 Xi^-$) decays. The significance of the $Omega(2012)^-$ signal is 4.2$sigma$ after considering the systematic uncertainties. The ratio of the branching fraction of $Omega_{c}^{0} to pi^{+} Omega(2012)^- to pi^+ (bar{K}Xi)^{-}$ relative to that of $Omega_{c}^{0} to pi^{+} Omega^-$ is calculated to be 0.220 $pm$ 0.059(stat.) $pm$ 0.035(syst.). The individual ratios of the branching fractions of the two isospin modes are also determined, and found to be ${cal B}(Omega_{c}^0 to pi^+ Omega(2012)^-) times {cal B}(Omega(2012)^- to K^-Xi^0)/{cal B}(Omega_{c}^0 to pi^+ K^- Xi^0)$ = (9.6 $pm$ 3.2(stat.) $pm$ 1.8(syst.))% and ${cal B}(Omega_{c}^0 to pi^+ Omega(2012)^-) times {cal B}(Omega(2012)^- to bar{K}^0 Xi^-)/{cal B}(Omega_{c}^0 to pi^+ bar{K}^0 Xi^-)$ = (5.5 $pm$ 2.8(stat.) $pm$ 0.7(syst.))%.
We present a statistical analysis for the characteristics and spatial evolution of the interplanetary discontinuities (IDs) in the solar wind, from 0.13 to 0.9 au, by using the Parker Solar Probe measurements on Orbits 4 and 5. 3948 IDs have been col
lected, including 2511 rotational discontinuities (RDs) and 557 tangential discontinuities (TDs), with the remnant unidentified. The statistical results show that (1) the ID occurrence rate decreases from 200 events/day at 0.13 au to 1 events/day at 0.9 au, following a spatial scaling r-2.00, (2) the RD to TD ratio decreases quickly with the heliocentric distance, from 8 at r<0.3 au to 1 at r>0.4 au, (3) the magnetic field tends to rotate across the IDs, 45{deg} for TDs and 30{deg} for RDs in the pristine solar wind within 0.3 au, (4) a special subgroup of RDs exist within 0.3 au, characterized by small field rotation angles and parallel or antiparallel propagations to the background magnetic fields, (5) the TD thicknesses normalized by local ion inertial lengths (di) show no clear spatial scaling and generally range from 5 to 35 di, and the normalized RD thicknesses follow r-1.09 spatial scaling, (6) the outward (anti-sunward) propagating RDs predominate in all RDs, with the propagation speeds in the plasma rest frame proportional to r-1.03. This work could improve our understandings for the ID characteristics and evolutions and shed light on the study of the turbulent environment in the pristine solar wind.
Readout errors on near-term quantum computers can introduce significant error to the empirical probability distribution sampled from the output of a quantum circuit. These errors can be mitigated by classical postprocessing given the access of an exp
erimental emph{response matrix} that describes the error associated with measurement of each computational basis state. However, the resources required to characterize a complete response matrix and to compute the corrected probability distribution scale exponentially in the number of qubits $n$. In this work, we modify standard matrix inversion techniques using two perturbative approximations with significantly reduced complexity and bounded error when the likelihood of high order bitflip events is strongly suppressed. Given a characteristic error rate $q$, our first method recovers the probability of the all-zeros bitstring $p_0$ by sampling only a small subspace of the response matrix before inverting readout error resulting in a relative speedup of $text{poly}left(2^{n} / big(begin{smallmatrix} n w end{smallmatrix}big)right)$, which we motivate using a simplified error model for which the approximation incurs only $O(q^w)$ error for some integer $w$. We then provide a generalized technique to efficiently recover full output distributions with $O(q^w)$ error in the perturbative limit. These approximate techniques for readout error correction may greatly accelerate near term quantum computing applications.
