Researcher profile

Q. Wu

Q. Wu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
17topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2023arXiv

Flat optical conductivity in the topological kagome magnet TbMn$_6$Sn$_6$

Kagome magnet TbMn$_6$Sn$_6$ is a new type of topological material that is known to support exotic quantum magnetic states. Experimental work has identified that TbMn$_6$Sn$_6$ hosts Dirac electronic states that could lead to topological and Chern quantum phases, but the optical response of the Dirac fermions of TbMn$_6$Sn$_6$ and its properties remain to be explored. Here, we perform optical spectroscopy measurement combined with first-principles calculations on single-crystal sample of TbMn$_6$Sn$_6$ to investigate the associated exotic phenomena. TbMn$_6$Sn$_6$ exhibits frequency-independent optical conductivity spectra in a broad range from 1800 to 3000 cm$^{-1}$ (220-370 meV) in experiments. The theoretical band structures and optical conductivity spectra are calculated with several shifted Fermi energy to compare with the experiment. The theoretical spectra with 0.56 eV shift for Fermi energy are well consistent with our experimental results. Besides, the massive quasi-two-dimensional (quasi-2D) Dirac bands, which have linear band dispersion in $k_x$-$k_y$ plane and no band dispersion along the $k_z$ direction, exist close to the shifted Fermi energy. According to tight-binding model analysis, the quasi-2D Dirac bands give rise to a flat optical conductivity, while its value is smaller than, about one tenth of, that from the calculations and experiments. It indicates that the other trivial bands also contribute to the flat optical conductivity.

preprint2022arXiv

An 8\% Determination of the Hubble Constant from localized Fast Radio Bursts

The $Λ$CDM model successfully explains the majority of cosmological observations. However, the $Λ$CDM model is challenged by Hubble tension, a remarkable difference of Hubble constant $H_0$ between measurements from local probe and the prediction from Planck cosmic microwave background observations under $ Λ$CDM model. So one urgently needs new distance indicators to test the Hubble tension. Fast radio bursts (FRBs) are millisecond-duration pulses occurring at cosmological distances, which are attractive cosmological probes. However, there is a thorny problem that the dispersion measures (DMs) contributed by host galaxy and the inhomogeneities of intergalactic medium cannot be exactly determined from observations. Previous works assuming fixed values for them bring uncontrolled systematic error in analysis. A reasonable approach is to handle them as probability distributions extracted from cosmological simulations. Here we report a measurement of ${H_0} = 64.67^{+5.62}_{-4.66} {\rm \ km \ s^{-1} \ Mpc^{-1}}$ using fourteen localized FRBs, with an uncertainty of 8.7\% at 68.3 per cent confidence. Thanks to the high event rate of FRBs and localization capability of radio telescopes (i.e., Australian Square Kilometre Array Pathfinder and Very Large Array), future observations of a reasonably sized sample ($\sim$100 localized FRBs) will provide a new way of measuring $H_0$ with a high precision ($\sim$2.6\%) to test the Hubble tension.

preprint2022arXiv

Evidence for mechanical softening-hardening dual anomaly in transition metals from shock compressed vanadium

Solid usually becomes harder and tougher under compression, and turns softer at elevated temperature. Recently, compression-induced softening and heating-induced hardening (CISHIH) dual anomaly was predicted in group VB elements such as vanadium. Here, the evidence for this counterintuitive phenomenon is reported. By using accurate high-temperature high-pressure sound velocities measured at Hugoniot states generated by shock-waves, together with first-principles calculations, we observe not only the prominent compression-induced sound velocity reduction, but also strong heating-induced sound velocity enhancement, in shocked vanadium. The former corresponds to the softening in shear modulus by compression, whereas the latter reflects the reverse hardening by heat. These experiments also unveil another anomaly in Young's modulus that wasn't reported before. Based on the experimental and theoretical data, we infer that vanadium might transition from BCC into two different rhombohedral (RH1 and RH2) phases at about 79GPa and 116GPa along the Hugoniot, respectively, which implies a dramatic difference in static and dynamic loading, as well as the significance of deviatoric stress and rate-relevant effects in high-pressure phase transition dynamics.

preprint2022arXiv

First measurement of high-energy reactor antineutrinos at Daya Bay

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12~MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10~MeV is rejected with a significance of 6.2 standard deviations. A 29\% antineutrino flux deficit in the prompt energy region of 8-11~MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-$Q_β$ isotopes in commercial reactors.

preprint2022arXiv

Joint Determination of Reactor Antineutrino Spectra from $^{235}$U and $^{239}$Pu Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of $^{235}$U and $^{239}$Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of $^{235}$U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant $^{235}$U and $^{239}$Pu isotopes and improves the uncertainty of the $^{235}$U spectral shape to about 3\%. The ${}^{235}$U and $^{239}$Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the $^{235}$U and $^{239}$Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

preprint2022arXiv

Linear-in-Frequency Optical Conductivity over a broad range in the three-dimensional Dirac semimetal candidate Ir$_2$In$_8$Se

The optical conductivity of the new Dirac semimetal candidate Ir$_2$In$_8$Se is measured in a frequency range from 40 to 30000 cm$^{-1}$ at temperatures from 300 K down to 10 K. The measurement reveals that the compound is a low carrier density metal. We find that the real part of the conductivity $σ_1(ω)$ is linear in frequency over a broad range from 500 to 4000 cm$^{-1}$ at 300 K and varies slightly with cooling. This linearity strongly suggests the presence of three-dimensional linear electronic bands with band crossings near the Fermi level. Band structure calculations indicate the presence of type-II Dirac points. By comparing our data with the optical conductivity computed from the band structure, we conclude that the observed linear dependence mainly originates from the Dirac cones and the transition between the Dirac cones and the next lower bands. In addition, a weak energy gap feature is resolved below the charge density wave phase transition temperature in reflectivity spectra. An enhanced structure arising from the imperfect Fermi surface nesting is identified in the electronic susceptibility function, suggesting a Fermi surface nesting driven instability.

preprint2022arXiv

Optical spectroscopy and ultrafast pump-probe study of a quasi-one-dimensional charge density wave in CuTe

CuTe is a two-dimensional (2D) layered material, yet forming a quasi-one-dimensional (quasi-1D) charge-density-wave (CDW) along the a-axis in the ab-plane at high transition temperature $T_{CDW}=335$ K. However, the anisotropic properties of CuTe remain to be explored. Here we performed combined transport, polarized infrared reflectivity, and ultrafast pump-probe spectroscopy to investigate the underlying CDW physics of CuTe. Polarized optical measurement clearly revealed that an energy gap gradually forms along the a-axis upon cooling, while optical evidence of gap signature is absent along the b-axis, suggesting pronounced electronic anisotropy in this quasi-2D material. Time-resolved optical reflectivity measurement revealed that the amplitude and relaxation time of photo-excited quasiparticles change dramatically across the CDW phase transition. Taking fast Fourier transformation of the oscillation signals arising from collective excitations, we identify the 1.65-THz mode as the CDW amplitude mode, whose energy softens gradually at elevated temperatures. Consequently, we provide further evidence for the formation of completely anisotropic CDW order in CuTe, which is quite rare in quasi-2D materials.

preprint2022arXiv

Spin excitations in the quantum dipolar magnet Yb(BaBO$_3$)$_3$

We report results of magnetization, specific-heat and muon-spin relaxation measurements on single crystals of disorder-free Yb$^{3+}$ triangular lattice Yb(BaBO$_3$)$_3$. The magnetization experiments show anisotropic magnetic properties with Curie-Weiss temperatures $θ_{\perp}=-1.40$~K ($H \perp c$) and $θ_{\parallel}=-1.16$~K ($H \parallel c$) determined from low temperature data. The absence of both long-range antiferromagnetic order and spin freezing is confirmed down to 0.27 K at zero field. A two-level Schottky anomaly due to the opening of the ground-state Kramers doublet is observed from the low-temperature specific-heat measurements when the applied magnetic fields $μ_0H >0.7$~T. At zero field, the increase of both $C_{\rm mag}/T$ and the muon spin relaxation rate $λ$ below 1~K is due to the electronic spin excitations, which often exist in quantum magnets where dipole-dipole interaction creates an anisotropy of magnetic properties. The spin excitation is also supported by the unusual maximum of field dependence of $λ$ due to the field-induced increase of the density of excitations. We argue that dipolar interaction is dominant and induces the spin dynamics in the quantum magnet Yb(BaBO$_3$)$_3$.

preprint2021arXiv

A parsec-scale faint jet in the nearby changing-look Seyfert galaxy Mrk 590

Broad Balmer emission lines in active galactic nuclei (AGN) may display dramatic changes in amplitude, even disappearance and re-appearance in some sources. As a nearby galaxy at a redshift of z = 0.0264, Mrk 590 suffered such a cycle of Seyfert type changes between 2006 and 2017. Over the last fifty years, Mrk 590 also underwent a powerful continuum outburst and a slow fading from X-rays to radio wavelengths with a peak bolometric luminosity reaching about ten per cent of the Eddington luminosity. To track its past accretion and ejection activity, we performed very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.6 GHz in 2015. The EVN observations reveal a faint (~1.7 mJy) radio jet extending up to ~2.8 mas (projected scale ~1.4 pc) toward north, and probably resulting from the very intensive AGN activity. To date, such a parsec-scale jet is rarely seen in the known changing-look AGN. The finding of the faint jet provides further strong support for variable accretion as the origin of the type changes in Mrk 590.

preprint2021arXiv

Optical spectroscopy and ultrafast pump-probe study of structural phase transition in 1T'-TaTe2

1T'-TaTe2 exhibits an intriguing first-order structural phase transition at around 170 K. Understanding the electronic structural properties is a crucial way to comprehend the origin of the structural phase transition. We performed a combined optical and ultrafast pump-probe study on the compound across the transition temperature. The phase transition leads to abrupt changes of both optical spectra and ultrafast electronic relaxation dynamics. The measurements revealed a sudden reconstruction of band structure. We elaborate that the phase transition is of the first order and can not be attributed to the conventional density-wave type instability. Our work is illuminating for understanding the origin of the structural phase transition.

preprint2020arXiv

A New Method to Measure Hubble Parameter $H(z)$ using Fast Radio Bursts

The Hubble parameter $H(z)$ is directly related to the expansion of our Universe. It can be used to study dark energy and constrain cosmology models. In this paper, we propose that $H(z)$ can be measured using fast radio bursts (FRBs) with redshift measurements. We use dispersion measures contributed by the intergalactic medium, which is related to $H(z)$, to measure Hubble parameter. We find that 500 mocked FRBs with dispersion measures and redshift information can accurately measure Hubble parameters using Monte Carlo simulation. The maximum deviation of $H(z)$ from standard $Λ$CDM model is about 6\% at redshift $z= 2.4$. We also test our method using Monte Carlo simulation. Kolmogorov-Smirnov (K-S) test is used to check the simulation. The $p$-value of K-S test is 0.23, which confirms internal consistency of the simulation. In future, more localizations of FRBs make it as an attractive cosmological probe.

preprint2020arXiv

An {\it ab-initio} Gamow shell model approach with a core

Gamow shell model (GSM) is usually performed within the Woods-Saxon (WS) basis in which the WS parameters need to be determined by fitting experimental single-particle energies including their resonance widths. In the multi-shell case, such a fit is difficult due to the lack of experimental data of cross-shell single-particle energies and widths. In this paper, we develop an {\it ab-initio} GSM by introducing the Gamow Hartree-Fock (GHF) basis that is obtained using the same interaction as the one used in the construction of the shell-model Hamiltonian. GSM makes use of the complex-momentum Berggren representation, then including resonance and continuum components. Hence, GSM gives a good description of weakly bound and unbound nuclei. Starting from chiral effective field theory and employing many-body perturbation theory (MBPT) (called nondegenerate $\hat Q$-box folded-diagram renormalization) in the GHF basis, a multi-shell Hamiltonian ({\it sd-pf} shells in this work) can be constructed. The single-particle energies and their resonance widths can also been obtained using MBPT. We investigated $^{23-28}$O and $^{23-31}$F isotopes, for which multi-shell calculations are necessary. Calculations show that continuum effects and the inclusion of the {\it pf} shell are important elements to understand the structure of nuclei close to and beyond driplines.

preprint2020arXiv

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the $θ_{μe}$ mixing angle are derived that constitute the most stringent limits to date over five orders of magnitude in the sterile mass-squared splitting $Δm^2_{41}$, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL$_s$ for $Δm^2_{41}<5\,$eV$^2$.Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL$_s$ for $Δm^2_{41}$ $<$ 1.2 eV$^2$.

preprint2020arXiv

Neutron-rich calcium isotopes within realistic Gamow shell model calculations with continuum coupling

Based on the realistic nuclear force of the high-precision CD-Bonn potential, we have performed comprehensive calculations for neutron-rich calcium isotopes using the Gamow shell model (GSM) which includes resonance and continuum. The realistic GSM calculations produce well binding energies, one- and two-neutron separation energies, predicting that $^{57}$Ca is the heaviest bound odd isotope and $^{70}$Ca is the dripline nucleus. Resonant states are predicted, which provides useful information for future experiments on particle emissions in neutron-rich calcium isotopes. Shell evolutions in the calcium chain around neutron numbers \textit{N} = 32, 34 and 40 are understood by calculating effective single-particle energies, the excitation energies of the first $2^+$ states and two-neutron separation energies. The calculations support shell closures at $^{52}$Ca (\textit{N} = 32) and $^{54}$Ca (\textit{N} = 34) but show a weakening of shell closure at $^{60}$Ca (\textit{N} = 40). The possible shell closure at $^{70}$Ca (\textit{N} = 50) is predicted.

preprint2020arXiv

Nuclear multipole responses from chiral effective field theory interaction

We probe nuclear multipole resonances in the framework of the random-phase approximation by using the interaction obtained from the chiral effective field theory. The three-nucleon force is included in a form of the in-medium two-nucleon interaction which was derived from the chiral three-nucleon force. The isoscalar monopole, isoscalar dipole, isovector dipole and isoscalar quadrupole resonances of the closed-shell $^{56,68,78}$Ni have been investigated. The calculations reasonably reproduce the experimental multipole resonances of $^{56,68}$Ni, and well describe the pygmy dipole resonance and dipole polarizability measured in $^{68}$Ni. The multipole resonances of $^{78}$Ni, including pygmy dipole resonance and dipole polarizability, are predicted. The detailed effects of the tensor force and three-body force are analyzed by dissecting the chiral interaction. We find that in general the tensor force effect on electric giant resonances is not as significant as the effect from the three-body force, although the tensor force provides more than half of the binding energy. The effect from three-body force is strong in light nuclei. Particularly, three-body force is crucial for the formation of the pygmy resonance in calculations.

preprint2020arXiv

NWChem: Past, Present, and Future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principledriven methodologies to model complex chemical and materials processes. Over the last few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach and outlook.

preprint2020arXiv

Optimization of the JUNO liquid scintillator composition using a Daya Bay antineutrino detector

To maximize the light yield of the liquid scintillator (LS) for the Jiangmen Underground Neutrino Observatory (JUNO), a 20 t LS sample was produced in a pilot plant at Daya Bay. The optical properties of the new LS in various compositions were studied by replacing the gadolinium-loaded LS in one antineutrino detector. The concentrations of the fluor, PPO, and the wavelength shifter, bis-MSB, were increased in 12 steps from 0.5 g/L and <0.01 mg/L to 4 g/L and 13 mg/L, respectively. The numbers of total detected photoelectrons suggest that, with the optically purified solvent, the bis-MSB concentration does not need to be more than 4 mg/L. To bridge the one order of magnitude in the detector size difference between Daya Bay and JUNO, the Daya Bay data were used to tune the parameters of a newly developed optical model. Then, the model and tuned parameters were used in the JUNO simulation. This enabled to determine the optimal composition for the JUNO LS: purified solvent LAB with 2.5 g/L PPO, and 1 to 4 mg/L bis-MSB.

preprint2020arXiv

Orbital localization error of density functional theory in shear properties of vanadium and niobium

It is believed that the density functional theory (DFT) describes most elements with s, p and d orbitals very well, except some materials that having strongly localized and correlated valence electrons. In this work, we find that the widely employed exchange-correlation (xc) functionals, including LDA, GGA and meta-GGA, underestimate the shear modulus and phase stability of V and Nb greatly. The advanced hybrid functional that is usually better for correlated system, on the other hand, completely fails in these two simple metals. This striking failure is revealed due to the orbital localization error in GGA, which is further deteriorated by hybrid functionals. This observation is corroborated by a similar failure of DFT+U and van der Waals functionals when applied to V and Nb. To remedy this problem, an semi-empirical approach of DFT+J is proposed which can delocalize electrons by facilitating the on-site exchange. Furthermore, it is observed that including density derivatives slightly improves the performance of the semi-local functionals, with meta-GGA outperforms GGA, and the latter is better than LDA. This discovery indicates the possibility and necessity to include higher-order density derivatives beyond the Laplacian level for the purpose to remove the orbital localization error (mainly from d orbitals) and delocalization error (mainly from s and p orbitals) completely in V and Nb, so that to achieve a better description of their electronic structures. The same strategy can be applied to other d electron system and f electron system.

preprint2020arXiv

Thermodynamic anomalies and three distinct liquid-liquid transitions in warm dense liquid hydrogen

The properties of hydrogen at high pressure have wide implications in astrophysics and high-pressure physics. Its phase change in the liquid is variously described as a metallization, H2-dissociation, density discontinuity or plasma phase transition. It has been tacitly assumed that these phenomena coincide at a first-order liquid-liquid transition (LLT). In this work, the relevant pressure-temperature conditions are thoroughly explored with first-principles molecular dynamics. We show there is a large dependency on exchange-correlation functional and significant finite size effects. We use hysteresis in a number of measurable quantities to demonstrate a first-order transition up to a critical point, above which molecular and atomic liquids are indistinguishable. At higher temperature beyond the critical point, H2-dissociation becomes a smooth cross-over in the supercritical region that can be modelled by a pseudo-transition, where the H2-2H transformation is localized and does not cause a density discontinuity at metallization. Thermodynamic anomalies and counter-intuitive transport behavior of protons are also discovered even far beyond the critical point, making this dissociative transition highly relevant to the interior dynamics of Jovian planets. Below the critical point, simulation also reveals a dynamic H2-2H chemical equilibrium with rapid interconversion, showing that H2 and H are miscible. The predicted critical temperature lies well below the ionization temperature. Our calculations unequivocally demonstrate that there are three distinct regimes in the liquid-liquid transition of warm dense hydrogen.

preprint2020arXiv

Understanding FRB 200428 in the synchrotron maser shock model: consistency and possible challenge

Recently, the discovery of Galactic FRB 200428 associated with a X-ray burst (XRB) of SGR 1935+2154 has built a bridge between FRBs and magnetar activities. In this paper, we assume that the XRB occurs in the magnetar magnetosphere. We show that the observational properties of FRB 200428 and the associated XRB are consistent with the predictions of synchrotron maser emission at ultrarelativistic magnetized shocks, including radiation efficiency, similar energy occurrence frequency distributions, and simultaneous arrive times. It requires that the upstream medium is a mildly relativistic baryonic shell ejected by a previous flare. The energy injection by flares responsible for the radio bursts will produce a magnetar wind nebula, which has been used to explain the persistent radio source associated FRB 121102. We find that the radio continuum around SGR 1935+2154 can be well understood in the magnetar wind nebula model, by assuming the same energy injection rate $\dot{E} \propto t^{-1.37}$ as FRB 121102. The required baryonic mass is also estimated form the observations of FRB 121102 by GBT and FAST. By assuming the same radiation efficiency $η\sim 10^{-5}$, the total baryonic mass ejected from the central magnetar is about 0.005 solar mass. This value is much larger than the typical mass of a magnetar outer crust, but is comparable to the total mass of a magnetar crust.

preprint2019arXiv

Ultrafast Quasiparticle Dynamics and Electron-Phonon Coupling in (Li0.84Fe0.16)OHFe0.98Se

Distinctive superconducting behaviors between bulk and monolayer FeSe make it challenging to obtain a unified picture of all FeSe-based superconductors. We investigate the ultrafast quasiparticle (QP) dynamics of an intercalated superconductor (Li1-xFex)OHFe1-ySe, which is a bulk crystal but shares a similar electronic structure with single-layer FeSe on SrTiO3. We obtain the electron-phonon coupling (EPC) constant λA1g (0.22 +/- 0.04), which well bridges that of bulk FeSe crystal and single-layer FeSe on SrTiO3. Moreover, we find that such a positive correlation between λA1g and superconducting Tc holds among all known FeSe-based superconductors, even in line with reported FeAs-based superconductors. Our observation indicates possible universal role of EPC in the superconductivity of all known categories of iron-based superconductors, which is a critical step towards achieving a unified superconducting mechanism for all iron-based superconductors.

preprint2018arXiv

Resonances of unbound quantum many-body systems of nuclei

Resonance is a general phenomenon which can happen in classic or quantum systems. An unbound many-body quantum system can undergo a self-resonant process. It has long been a challenge how to describe unbound many-body quantum systems in resonances. In this paper, we develop a novel first-principles method that is capable of describing resonant quantum systems. We exploit, for the first time, the advanced in-medium similarity renormalization group (IMSRG) in the complex-energy Gamow-Berggren representation with resonance and non-resonant continuum. The ab initio Gamow IMSRG has broad applications, such as, to the electromagnetic-interaction systems of atoms, molecules or quantum dots, and strong-interaction atomic nuclei. In the present work, we apply the method to loosely bound or unbound nuclear systems. Carbon and oxygen isotopes have been investigated with an optimized chiral effective field theory potential. The resonant states observed in neutron-rich 22O and 24O are well reproduced. The loose halo structure of the Borromean nucleus 22C is clearly seen by the density calculation, in which the continuum s-waves play a crucial role. Further, we predict low-lying resonant excited states in 22C. This method provides rigorous and tractable theoretical calculations for weakly-bound or unbound open quantum systems.