Search

Strong spin-orbit coupling and relatively weak hyperfine interactions make germanium hole spin qubits a promising candidate for semiconductor quantum processors. The two-dimensional hole gas structure of strained Ge quantum wells serves as the primary material platform for spin hole qubits.A low disorder material environment is essential for this process. In this work, we fabricated a Ge/SiGe heterojunction with a 60 nm buried quantum well layer on a Si substrate using reduced pressure chemical vapor deposition technology. At a temperature of 16 mK, when the carrier density is 1.87*10^11/cm2, we obtained a mobility as high as 308.64*10^4cm2/Vs. Concurrently, double quantum dot and planar germanium coupling with microwave cavities were also successfully achieved.This fully demonstrates that this structure can be used for the preparation of higher-performance hole spin qubits.

Speech enhancement is designed to enhance the intelligibility and quality of speech across diverse noise conditions. Recently, diffusion model has gained lots of attention in speech enhancement area, achieving competitive results. Current diffusion-based methods blur the signal with isotropic Gaussian noise and recover clean speech from the prior. However, these methods often suffer from a substantial computational burden. We argue that the computational inefficiency partially stems from the oversight that speech enhancement is not purely a generative task; it primarily involves noise reduction and completion of missing information, while the clean clues in the original mixture do not need to be regenerated. In this paper, we propose a method that introduces noise with anisotropic guidance during the diffusion process, allowing the neural network to preserve clean clues within noisy recordings. This approach substantially reduces computational complexity while exhibiting robustness against various forms of noise and speech distortion. Experiments demonstrate that the proposed method achieves state-of-the-art results with only approximately 4.5 million parameters, a number significantly lower than that required by other diffusion methods. This effectively narrows the model size disparity between diffusion-based and predictive speech enhancement approaches. Additionally, the proposed method performs well in very noisy scenarios, demonstrating its potential for applications in highly challenging environments., Comment: We make reassessment and update the author list. All authors have approved this version of the manuscript

In this paper, we present our approach for the Track 1 of the Chinese Auditory Attention Decoding (Chinese AAD) Challenge at ISCSLP 2024. Most existing spatial auditory attention decoding (Sp-AAD) methods employ an isolated window architecture, focusing solely on global invariant features without considering relationships between different decision windows, which can lead to suboptimal performance. To address this issue, we propose a novel streaming decoding architecture, termed StreamAAD. In StreamAAD, decision windows are input to the network as a sequential stream and decoded in order, allowing for the modeling of inter-window relationships. Additionally, we employ a model ensemble strategy, achieving significant better performance than the baseline, ranking First in the challenge.

The spatial auditory attention decoding (Sp-AAD) technology aims to determine the direction of auditory attention in multi-talker scenarios via neural recordings. Despite the success of recent Sp-AAD algorithms, their performance is hindered by trial-specific features in EEG data. This study aims to improve decoding performance against these features. Studies in neuroscience indicate that spatial auditory attention can be reflected in the topological distribution of EEG energy across different frequency bands. This insight motivates us to propose Prototype Training, a neuroscience-inspired method for Sp-AAD. This method constructs prototypes with enhanced energy distribution representations and reduced trial-specific characteristics, enabling the model to better capture auditory attention features. To implement prototype training, an EEGWaveNet that employs the wavelet transform of EEG is further proposed. Detailed experiments indicate that the EEGWaveNet with prototype training outperforms other competitive models on various datasets, and the effectiveness of the proposed method is also validated. As a training method independent of model architecture, prototype training offers new insights into the field of Sp-AAD.

In this paper, we studied the space-time estimates for the solution to the Schr\"odinger equation. By polynomial partitioning, induction arguments, bilinear to linear arguments and broad norm estimates, we set up several maximal estimates for the Schr\"odinger equation with high-frequency input data. By these maximal estimates, we obtain the sharp global space-time estimate when $n=2$ and improve the known results in the critical cases when $n\geq 3$. The maximal estimate for $n=2$ is also used to extend the results of the local space-time estimates for the solution to the Schr\"odinger equation.

In this paper, we consider the $L_x^p(\mathbb{R}^2)\rightarrow L_{x,u}^q(\mathbb{R}^2\times [1,2])$ estimate for the operator $T$ along a dilated plane curve $(ut,u\gamma(t))$, where $$Tf(x,u):=\int_{0}^{1}f(x_1-ut,x_2-u \gamma(t))\,\textrm{d}t,$$ $x:=(x_1,x_2)$ and $\gamma$ is a general plane curve satisfying some suitable smoothness and curvature conditions. We show that $T$ is $L_x^p(\mathbb{R}^2)$ to $L_{x,u}^q(\mathbb{R}^2\times [1,2])$ bounded whenever $(\frac{1}{p},\frac{1}{q})\in \square \cup \{(0,0)\}\cup \{(\frac{2}{3},\frac{1}{3})\}$ and $1+(1 +\omega)(\frac{1}{q}-\frac{1}{p})>0$, where the trapezium $\square:=\{(\frac{1}{p},\frac{1}{q}):\ \frac{2}{p}-1\leq\frac{1}{q}\leq \frac{1}{p}, \frac{1}{q}>\frac{1}{3p}, \frac{1}{q}>\frac{1}{p}-\frac{1}{3}\}$ and $\omega:=\limsup_{t\rightarrow 0^{+}}\frac{\ln|\gamma(t)|}{\ln t}$. This result is sharp except for some borderline cases. On the other hand, in a smaller $(\frac{1}{p},\frac{1}{q})$ region, we also obtain the almost sharp estimate $T : L_x^p(\mathbb{R}^2)\rightarrow L_{x}^q(\mathbb{R}^2)$ uniformly for $u\in [1,2]$. These results imply that the operator $T$ has the so called local smoothing phenomenon, i.e., the $L^q$ integral about $u$ on $[1,2]$ extends the region of $(\frac{1}{p},\frac{1}{q})$ in uniform estimate $T : L_x^p(\mathbb{R}^2)\rightarrow L_{x}^q(\mathbb{R}^2)$.

As the latest advancements in natural language processing, large language models (LLMs) have achieved human-level language understanding and generation abilities in many real-world tasks, and even have been regarded as a potential path to the artificial general intelligence. To better facilitate research on LLMs, many open-source LLMs, such as Llama 2 and Falcon, have recently been proposed and gained comparable performances to proprietary models. However, these models are primarily designed for English scenarios and exhibit poor performances in Chinese contexts. In this technical report, we propose YAYI 2, including both base and chat models, with 30 billion parameters. YAYI 2 is pre-trained from scratch on a multilingual corpus which contains 2.65 trillion tokens filtered by our pre-training data processing pipeline. The base model is aligned with human values through supervised fine-tuning with millions of instructions and reinforcement learning from human feedback. Extensive experiments on multiple benchmarks, such as MMLU and CMMLU, consistently demonstrate that the proposed YAYI 2 outperforms other similar sized open-source models.

A traffic system is a random and complex large system, which is difficult to conduct repeated modelling and control research in a real traffic environment. With the development of automatic driving technology, the requirements for testing and evaluating the development of automatic driving technology are getting higher and higher, so the application of computer technology for traffic simulation has become a very effective technical means. Based on the micro-traffic flow modelling, this paper adopts the vehicle motion model based on cellular automata and the theory of bicycle intelligence to build the simulation environment of autonomous vehicle flow. The architecture of autonomous vehicles is generally divided into a perception system, decision system and control system. The perception system is generally divided into many subsystems, responsible for autonomous vehicle positioning, obstacle recognition, traffic signal detection and recognition and other tasks. Decision systems are typically divided into many subsystems that are responsible for tasks such as path planning, path planning, behavior selection, motion planning, and control. The control system is the basis of the selfdriving car, and each control system of the vehicle needs to be connected with the decision-making system through the bus, and can accurately control the acceleration degree, braking degree, steering amplitude, lighting control and other driving actions according to the bus instructions issued by the decision-making system, so as to achieve the autonomous driving of the vehicle.

The low gain avalanche detectors (LGADs) are thin sensors with fast charge collection which in combination with internal gain deliver an outstanding time resolution of about 30 ps. High collision rates and consequent large particle rates crossing the detectors at the upgraded Large Hadron Collider (LHC) in 2028 will lead to radiation damage and deteriorated performance of the LGADs. The main consequence of radiation damage is loss of gain layer doping (acceptor removal) which requires an increase of bias voltage to compensate for the loss of charge collection efficiency and consequently time resolution. The Institute of High Energy Physics (IHEP), Chinese Academy of Sciences (CAS) has developed a process based on the Institute of Microelectronics (IME), CAS capability to enrich the gain layer with carbon to reduce the acceptor removal effect by radiation. After 1 MeV neutron equivalent fluence of 2.5$\times$10$^{15}$ n$_{eq}$/cm$^{2}$, which is the maximum fluence to which sensors will be exposed at ATLAS High Granularity Timing Detector (HGTD), the IHEP-IME second version (IHEP-IMEv2) 50 $\mu$m LGAD sensors already deliver adequate charge collection > 4 fC and time resolution < 50 ps at voltages < 400 V. The operation voltages of these 50 $\mu$m devices are well below those at which single event burnout may occur.

We fabricate an undoped Ge quantum well under 30 nm Ge0.8Si0.2 shallow barrier with reverse grading technology. The under barrier is deposited by Ge0.8Si0.2 followed by Ge0.9Si0.1 so that the variation of Ge content forms a sharp interface which can suppress the threading dislocation density penetrating into undoped Ge quantum well. And the Ge0.8Si0.2 barrier introduces enough in-plane parallel strain -0.41% in the Ge quantum well. The heterostructure field-effect transistors with a shallow buried channel get a high two-dimensional hole gas (2DHG) mobility over 2E6 cm2/Vs at a low percolation density of 2.51 E-11 cm2. We also discover a tunable fractional quantum Hall effect at high densities and high magnetic fields. This approach defines strained germanium as providing the material basis for tuning the spin-orbit coupling strength for fast and coherent quantum computation., Comment: 11 pages, 5 figures

Preliminary mission design requires an efficient and accurate approximation to the low-thrust rendezvous trajectories, which might be generally three-dimensional and involve multiple revolutions. In this paper, a new shaping method using cubic spline functions is developed for the analytical approximation, which shows advantages in the optimality and computational efficiency. The rendezvous constraints on the boundary states and transfer time are all satisfied analytically, under the assumption that the boundary conditions and segment numbers of cubic spline functions are designated in advance. Two specific shapes are then formulated according to whether they have free optimization parameters. The shape without free parameters provides an efficient and robust estimation, while the other one allows a subsequent optimization for the satisfaction of additional constraints such as the constraint on the thrust magnitude. Applications of the proposed method in combination with the particle swarm optimization algorithm are discussed through two typical interplanetary rendezvous missions, that is, an inclined multi-revolution trajectory from the Earth to asteroid Dionysus and a multi-rendezvous trajectory of sample return. Simulation examples show that the proposed method is superior to existing methods in terms of providing good estimation for the global search and generating suitable initial guess for the subsequent trajectory optimization.

Low Gain Avalanche Diode (LGAD) is applied for the High-Granularity Timing Detector (HGTD), and it will be used to upgrade the ATLAS experiment. The first batch IHEP-IME LGAD sensors were designed by the Institute of High Energy Physics (IHEP) and fabricated by the Institute of Microelectronics (IME). Three IHEP-IME sensors (W1, W7 and W8) were irradiated by the neutrons up to the fluence of 2.5 x 10^15 n_eq/cm^2 to study the effect of the shallow carbon and deep N++ layer on the irradiation hardness. Taking W7 as a reference, W1 has an extra shallow carbon applied, and W8 has a deeper N++ layer., Comment: This work has been submitted to the IEEE for possible publication

This letter proposes a new blind source separation (BSS) framework termed minimum variance independent component analysis (MVICA), which can potentially achieve the maximum output signal-to-interference ratio (SIR) while also allowing more flexibility in real implementations. The statistical independence assumption has been the foundation of the most dominant BSS techniques in recent decades. However, this assumption does not always hold true and the accurate probabilistic modeling of source is inherently difficult. To overcome these limitations and improve the separation performance, the MVICA framework is rigorously derived by optimizing the design of these independence-based BSS algorithms with the maximum SIR criterion. A deep neural network-supported implementation of MVICA is subsequently described. Experimental results under various conditions show the superiority of MVICA over the state-of-the-art BSS algorithms, in terms of not only SIR but also signal-to-distortion ratio and automatic speech recognition rate., Comment: 5 pages, 2 figures

This paper shows the simulation and test results of 50um thick Low Gain Avalanche Detectors (LGAD) sensors designed by the Institute of High Energy Physics (IHEP) and fabricated by the Institute of Microelectronics of the Chinese Academy of Sciences (IME). Three wafers have been produced with four different gain layer implant doses each. Different production processes, including variation in the n++ layer implant energy and carbon co-implantation were used. Test results show that the IHEP-IME sensors with the higher dose of gain layer have lower breakdown voltages and higher gain layer voltages from capacitance-voltage properties, which are consistent with the TCAD simulation. Beta test results show that the time resolution of IHEP-IME sensors is better than 35ps when operated at high voltage and the collected charges of IHEP-IME sensors are larger than 15fC before irradiation, which fulfill the required specifications of sensors before irradiations for the ATLAS HGTD project., Comment: 11 pages,8 figures

Serverless computing is increasingly popular because of its lower cost and easier deployment. Several cloud service providers (CSPs) offer serverless computing on their public clouds, but it may bring the vendor lock-in risk. To avoid this limitation, many open-source serverless platforms come out to allow developers to freely deploy and manage functions on self-hosted clouds. However, building effective functions requires much expertise and thorough comprehension of platform frameworks and features that affect performance. It is a challenge for a service developer to differentiate and select the appropriate serverless platform for different demands and scenarios. Thus, we elaborate the frameworks and event processing models of four popular open-source serverless platforms and identify their salient idiosyncrasies. We analyze the root causes of performance differences between different service exporting and auto-scaling modes on those platforms. Further, we provide several insights for future work, such as auto-scaling and metric collection.

Quick network address translation (NAT) is proposed to improve the network performance of the NAT system on the commodity server by three ways. First, the quick NAT search algorithm is designed to use the Hash search instead of the sequential search to reduce latency when looking up the NAT rule table. Second, to leverage the power of the multi-core central processing unit (CPU) and the multi-queue network interface card, Quick NAT enables multiple CPU cores to process in parallel. The localized connection tracking table and the compare-and-swap based lock-free NAT Hash tables are designed to eliminate the lock overhead. Third, Quick NAT uses the polling and zero-copy delivery to reduce the cost of interrupt and packet copies. The evaluation results show that Quick NAT obtains high scalability and line-rate throughput on the commodity server., Comment: in Chinese language

NAT gateway is an important network system in today's IPv4 network when translating a private IPv4 address to a public address. However, traditional NAT system based on Linux Netfilter cannot achieve high network throughput to meet modern requirements such as data centers. To address this challenge, we improve the network performance of NAT system by three ways. First, we leverage DPDK to enable polling and zero-copy delivery, so as to reduce the cost of interrupt and packet copies. Second, we enable multiple CPU cores to process in parallel and use lock-free hash table to minimize the contention between CPU cores. Third, we use hash search instead of sequential search when looking up the NAT rule table. Evaluation shows that our Quick NAT system significantly improves the performance of NAT on commodity platforms.

Long Short Term Memory(LSTM) models are the building blocks of many state-of-the-art natural language processing(NLP) and speech enhancement(SE) algorithms. However, there are a large number of parameters in an LSTM model. This usually consumes a large number of resources to train the LSTM model. Also, LSTM models suffer from computational inefficiency in the inference phase. Existing model compression methods (e.g., model pruning) can only discriminate based on the magnitude of model parameters, ignoring the issue of importance distribution based on the model information. Here we introduce the MPO decomposition, which describes the local correlation of quantum states in quantum many-body physics and is used to represent the large model parameter matrix in a neural network, which can compress the neural network by truncating the unimportant information in the weight matrix. In this paper, we propose a matrix product operator(MPO) based neural network architecture to replace the LSTM model. The effective representation of neural networks by MPO can effectively reduce the computational consumption of training LSTM models on the one hand, and speed up the computation in the inference phase of the model on the other hand. We compare the MPO-LSTM model-based compression model with the traditional LSTM model with pruning methods on sequence classification, sequence prediction, and speech enhancement tasks in our experiments. The experimental results show that our proposed neural network architecture based on the MPO approach significantly outperforms the pruning approach., Comment: 2 figures, 5 tables

The deep neural network (DNN) based speech enhancement approaches have achieved promising performance. However, the number of parameters involved in these methods is usually enormous for the real applications of speech enhancement on the device with the limited resources. This seriously restricts the applications. To deal with this issue, model compression techniques are being widely studied. In this paper, we propose a model compression method based on matrix product operators (MPO) to substantially reduce the number of parameters in DNN models for speech enhancement. In this method, the weight matrices in the linear transformations of neural network model are replaced by the MPO decomposition format before training. In experiment, this process is applied to the causal neural network models, such as the feedforward multilayer perceptron (MLP) and long short-term memory (LSTM) models. Both MLP and LSTM models with/without compression are then utilized to estimate the ideal ratio mask for monaural speech enhancement. The experimental results show that our proposed MPO-based method outperforms the widely-used pruning method for speech enhancement under various compression rates, and further improvement can be achieved with respect to low compression rates. Our proposal provides an effective model compression method for speech enhancement, especially in cloud-free application., Comment: 11 pages, 6 figures, 7 tables

In this paper, we determine the $L^p(\mathbb{R})\times L^q(\mathbb{R})\rightarrow L^r(\mathbb{R})$ boundedness of the bilinear Hilbert transform $H_{\gamma}(f,g)$ along a convex curve $\gamma$ $$H_{\gamma}(f,g)(x):=\mathrm{p.\,v.}\int_{-\infty}^{\infty}f(x-t)g(x-\gamma(t)) \,\frac{\textrm{d}t}{t},$$ where $p$, $q$, and $r$ satisfy $\frac{1}{p}+\frac{1}{q}=\frac{1}{r}$, and $r>\frac{1}{2}$, $p>1$, and $q>1$. Moreover, the same $L^p(\mathbb{R})\times L^q(\mathbb{R})\rightarrow L^r(\mathbb{R})$ boundedness property holds for the corresponding (sub)bilinear maximal function $M_{\gamma}(f,g)$ along a convex curve $\gamma$ $$M_{\gamma}(f,g)(x):=\sup_{\varepsilon>0}\frac{1}{2\varepsilon}\int_{-\varepsilon}^{\varepsilon}|f(x-t)g(x-\gamma(t))| \,\textrm{d}t.$$

In this paper we consider the initial boundary value problem (IBVP) for the nonlinear biharmonic Schr\"odinger equation posed on a bounded interval $(0,L)$ with non-homogeneous Navier or Dirichlet boundary conditions, respectively. For Navier boundary IBVP, we set up its local well-posedness if the initial data lies in $H^s(0, L)$ with $s\geq 0$ and $s\neq n+1/2, n\in \mathbb{N}$, and the boundary data are selected from the appropriate spaces with optimal regularities, i.e., the $j$-th order data are chosen in $H_{loc}^{(s+3-j)/4}(\mathbb {R}^+)$, for $j=0,2$. For Dirichlet boundary IBVP the corresponding local well-posedness is obtained when $s>10/7$ and $s\neq n+1/2, n\in \mathbb{N}$, and the boundary data are selected from the appropriate spaces with optimal regularities, i.e., the $j$-th order data are chosen in $H_{loc}^{(s+3-j)/4}(\mathbb {R}^+)$, for $j=0,1$.

This paper shows $$ \sup_{f\in H^s(\mathbb{R}^n)}\dim _H\left\{x\in\mathbb{R}^n:\ \lim_{t\rightarrow0}e^{it(-\Delta)^\alpha}f(x)\neq f(x)\right\}\leq n+1-\frac{2(n+1)s}{n}\ \ \text{under}\ \ \begin{cases} n\geq2;\\ \alpha>\frac12; \frac{n}{2(n+1)}

In this paper, Theorems 1.1- 1.2 show that the Boussinesq operator $\mathcal{B}_tf$ converges pointwise to its initial data $f\in H^s(\mathbb{R})$ as $t\to 0$ provided $s\geq\frac{1}{4}$ -- more precisely -- on the one hand, by constructing a counterexample in $\mathbb{R}$ we discover that the optimal convergence index $s_{c,1}=\frac14$; on the other hand, we find that the Hausdorff dimension of the disconvergence set for $\mathcal{B}_tf$ is \begin{align*} \alpha_{1,\mathcal{B}}(s)&=\begin{cases} 1-2s&\ \ \text{as}\ \ \frac{1}{4}\leq s\leq\frac{1}{2};\\ 1 &\ \ \text{as}\ \ 0

Serverless computing is increasingly popular because of the promise of lower cost and the convenience it provides to users who do not need to focus on server management. This has resulted in the availability of a number of proprietary and open-source serverless solutions. We seek to understand how the performance of serverless computing depends on a number of design issues using several popular open-source serverless platforms. We identify the idiosyncrasies affecting performance (throughput and latency) for different open-source serverless platforms. Further, we observe that just having either resource-based (CPU and memory) or workload-based (request per second (RPS) or concurrent requests) auto-scaling is inadequate to address the needs of the serverless platforms.

We find an interesting fact that fictitious retrograde co-orbitals of Saturn, or small bodies inside the retrograde 1:1 resonance with Saturn, are highly unstable in our numerical simulations. It is shown that in the presence of Jupiter, the retrograde co-orbitals will get ejected from Saturn's co-orbital space within a timescale of 10 Myr. This scenario reminds us of the instability of Saturn Trojans caused by both the Great Inequality and the secular resonances. Therefore, we carry out in-depth inspections on both mechanisms and prove that the retrograde resonance overlap, raised by Great Inequality, cannot serve as an explanation for the instability of retrograde co-orbitals, due to the weakness of the retrograde 2:5 resonance with Jupiter at a low eccentricity. However, we discover that both $\nu_5$ and $\nu_6$ secular resonances contribute to the slow growth of the eccentricity, therefore, are possibly the primary causes of the instability inside Saturn's retrograde co-orbital space., Comment: 6 pages, 6 figures, accepted for publication in MNRAS

As the discoveries of more minor bodies in retrograde resonances with giant planets, such as 2015 BZ509 and 2006 RJ2, our curiosity about the Kozai-Lidov dynamics inside the retrograde resonance has been sparked. In this study, we focus on the 3D retrograde resonance problem and investigate how the resonant dynamics of a minor body impacts on its own Kozai-Lidov cycle. Firstly we deduce the action-angle variables and canonical transformations that deal with the retrograde orbit specifically. After obtaining the dominant Hamiltonian of this problem, we then carry out the numerical averaging process in closed form to generate phase-space portraits on a $e-\omega$ space. The retrograde 1:1 resonance is particularly scrutinized in detail, and numerical results from a CRTBP model shows a great agreement with the our semi-analytical portraits. On this basis, we inspect two real minor bodies currently trapped in retrograde 1:1 mean motion resonance. It is shown that they have different Kozai-Lidov states, which can be used to analyze the stability of their unique resonances. In the end, we further inspect the Kozai-Lidov dynamics inside the 2:1 and 2:5 retrograde resonance, and find distinct dynamical bifurcations of equilibrium points on phase-space portraits., Comment: 10 pages, 8 figures. Accepted for publication in MNRAS

Background --The objective of this study was to examine the association of routine blood test results with coronary heart disease (CHD) risk, to incorporate them into coronary prediction models and to compare the discrimination properties of this approach with other prediction functions. Methods and Results --This work was designed as a retrospective, single-center study of a hospital-based cohort. The 5060 CHD patients (2365 men and 2695 women) were 1 to 97 years old at baseline with 8 years (2009-2017) of medical records, 5051 health check-ups and 5075 cases of other diseases. We developed a two-layer Gradient Boosting Decision Tree(GBDT) model based on routine blood data to predict the risk of coronary heart disease, which could identify 86% of people with coronary heart disease. We built a dataset with 15,000 routine blood tests results. Using this dataset, we trained the two-layer GBDT model to classify healthy status, coronary heart disease and other diseases. As a result of the classification after machine learning, we found that the sensitivity of detecting the health data was approximately 93% for all data, and the sensitivity of detecting CHD was 93% for disease data that included coronary heart disease. On this basis, we further visualized the correlation between routine blood results and related data items, and there was an obvious pattern in health and coronary heart disease in all data presentations, which can be used for clinical reference. Finally, we briefly analyzed the results above from the perspective of pathophysiology. Conclusions --Routine blood data provides more information about CHD than what we already know through the correlation between test results and related data items. A simple coronary disease prediction model was developed using a GBDT algorithm, which will allow physicians to predict CHD risk in patients without overt CHD.

In this paper, the $L^2$ boundedness of the Hilbert transform along variable flat curve $(t,P(x_1)\gamma(t))$ $$H_{P,\gamma}f(x_1,x_2):=\mathrm{p.\,v.}\int_{-\infty}^{\infty}f(x_1-t,x_2-P(x_1)\gamma(t))\,\frac{\textrm{d}t}{t},\quad \forall\, (x_1,x_2)\in\mathbb{R}^2,$$ is studied, where $P$ is a real polynomial on $\mathbb{R}$. A new sufficient condition on the curve $\gamma$ is introduced.

Let $p\in (1,\infty)$. In this paper, for any given measurable function $u:\ \mathbb{R}\rightarrow \mathbb{R}$ and a generalized plane curve $\gamma$ satisfying some conditions, the $L^p(\mathbb{R}^2)$ boundedness of the Hilbert transform along the variable plane curve $u(x_1)\gamma$ $$H_{u,\gamma}f(x_1,x_2):=\mathrm{p.\,v.}\int_{-\infty}^{\infty}f(x_1-t,x_2-u(x_1)\gamma(t)) \,\frac{\textrm{d}t}{t}, \quad \forall\, (x_1,x_2)\in\mathbb{R}^2, $$ is obtained. At the same time, the $L^p(\mathbb{R})$ boundedness of the corresponding Carleson operator along the general curve $\gamma$ $$\mathcal{C}_{u,\gamma}f(x):=\mathrm{p.\,v.}\int_{-\infty}^{\infty}e^{iu(x)\gamma (t)}f(x-t)\,\frac{\textrm{d}t}{t}, \quad\forall\, x\in\mathbb{R}, $$ is also obtained. Moreover, all the bounds are independent of the measurable function $u$.

Asteroids in mean motion resonances with giant planets are common in the solar system, but it was not until recently that several asteroids in retrograde mean motion resonances with Jupiter and Saturn were discovered. A retrograde co-orbital asteroid of Jupiter, 2015 BZ509 is confirmed to be in a long-term stable retrograde 1:1 mean motion resonance with Jupiter, which gives rise to our interests in its unique resonant dynamics. In this paper, we investigate the phase-space structure of the retrograde 1:1 resonance in detail within the framework of the circular restricted three-body problem. We construct a simple integrable approximation for the planar retrograde resonance using canonical contact transformation and numerically employ the averaging procedure in closed form. The phase portrait of the retrograde 1:1 resonance is depicted with the level curves of the averaged Hamiltonian. We thoroughly analyze all possible librations in the co-orbital region and uncover a new apocentric libration for the retrograde 1:1 resonance inside the planet's orbit. We also observe the significant jumps in orbital elements for outer and inner apocentric librations, which are caused by close encounters with the perturber., Comment: 8 pages, 5 figures. Accepted for publication in AJ

Image deblurring, a.k.a. image deconvolution, recovers a clear image from pixel superposition caused by blur degradation. Few deep convolutional neural networks (CNN) succeed in addressing this task. In this paper, we first demonstrate that the minimum-mean-square-error (MMSE) solution to image deblurring can be interestingly unfolded into a series of residual components. Based on this analysis, we propose a novel iterative residual deconvolution (IRD) algorithm. Further, IRD motivates us to take one step forward to design an explicable and effective CNN architecture for image deconvolution. Specifically, a sequence of residual CNN units are deployed, whose intermediate outputs are then concatenated and integrated, resulting in concatenated residual convolutional network (CRCNet). The experimental results demonstrate that proposed CRCNet not only achieves better quantitative metrics but also recovers more visually plausible texture details compared with state-of-the-art methods., Comment: rejected by AAAI 2019

The stability and topological structure of equilibrium points in the potential field of the asteroid 101955 Bennu have been investigated with a variable density and rotation period. A dimensionless quantity is introduced for the nondimensionalization of the equations of motion, and this quantity can indicate the effect of both the rotation period and bulk density of the asteroid. Using the polyhedral model of the asteroid Bennu, the number and position of the equilibrium points are calculated and illustrated by a contour plot of the gravitational effective potential field. The topological structure and the stability of the equilibrium points are also investigated using the linearized method. The results show that there are nine equilibrium points in the potential field of the asteroid Bennu, eight in the exterior of the body and one in the interior of the body. Moreover, bifurcation will occur with a variation of the density and rotation period. Different equilibrium points will encounter each other and mix together. Thus, the number of equilibrium points will change. The stability and topological structure of the equilibrium points will also change because of the variation of the density and rotation period of the asteroid. When considering the error of the density of Bennu, the range of the dimensionless quantity covers the critical values that will lead to bifurcation. This means that the stability of the equilibrium points is uncertain, making the dynamical environment of Bennu much more complicated. These bifurcations can help better understand the dynamic environment of an irregular-shaped asteroid., Comment: 11 pages, 11 figures, 4 tables

A method for classifying orbits near asteroids under a polyhedral gravitational field is presented, and may serve as a valuable reference for spacecraft orbit design for asteroid exploration. The orbital dynamics near asteroids are very complex. According to the variation in orbit characteristics after being affected by gravitational perturbation during the periapsis passage, orbits near an asteroid can be classified into 9 categories: (1) surroundingto-surrounding, (2) surrounding-to-surface, (3) surroundingto-infinity, (4) infinity-to-infinity, (5) infinity-to-surface, (6) infinity-to-surrounding, (7) surface-to-surface, (8) surfaceto-surrounding, and (9) surface-to- infinity. Assume that the orbital elements are constant near the periapsis, the gravitation potential is expanded into a harmonic series. Then, the influence of the gravitational perturbation on the orbit is studied analytically. The styles of orbits are dependent on the argument of periapsis, the periapsis radius, and the periapsis velocity. Given the argument of periapsis, the orbital energy before and after perturbation can be derived according to the periapsis radius and the periapsis velocity. Simulations have been performed for orbits in the gravitational field of 216 Kleopatra. The numerical results are well consistent with analytic predictions., Comment: 26 pages, 10 figures

As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion., Comment: 66 pages, 24 figures, submitted to Icarus on 25/Aug/2016

We study global well-posedness for the Kadomtsev-Petviashvili II equation in three space dimensions with small initial data. The crucial points are new bilinear estimates and the definition of the function spaces. As by-product we obtain that all solutions to small initial data scatter., Comment: 33 pages, 1 figure, final version

This manuscript considers the main features of the nonlinear dynamics of multiple irregular celestial body systems. The gravitational potential, static electric potential, and magnetic potential are considered. Based on the three established potentials, we show that three conservative values exist for this system, including a Jacobi integral. The equilibrium conditions for the system are derived and their stability analyzed. The equilibrium conditions of a celestial system comprised of n irregular bodies are reduced to 12n minus 9 equations. The dynamical results are applied to simulate the motion of multiple-asteroid systems. The simulation is useful for the study of the stability of multiple irregular celestial body systems and for the design of spacecraft orbits to triple asteroid systems discovered in the solar system. The dynamical configurations of the five triple-asteroid systems 45 Eugenia, 87 Sylvia, 93 Minerva, 216 Kleopatra, and 136617 1994CC, and the six-body system 134340 Pluto are calculated and analyzed., Comment: 60 pages, 20 figures, The online version of this article (doi:10.1007/s10509-016-2884-z) contains supplementary material, which is available to authorized users. 6 videos, 5 rar files, 1 pdf file of figures

Recently, the spinning tethered system is regarded as a typical and fundamental space structure attracting great interest of the aerospace engineers, and has been discussed primarily for specific space missions in past decades, including on-orbit capture and propellantless orbit transfer etc. The present work studies the dynamical behaviours of a fast spinning tethered binary system under central gravitational field, and derives principles of the basic laws of orbital maneuver. Considering the characteristics of coupled librational and orbital motions, an averaging method is introduced to deal with the slow-fast system equation, thus a definite equivalent model is derived. The general orbit motion is completely determined analytically, including the orbit geometry, periodicity, conversations and moving region etc. Since the possibility of orbit control using tether reaction has been proved by previous studies, special attention is paid to the transportation mode of angular momentum and mechanical energy between the orbit and libration. The effect of tether length change on the orbit shape is verified both in the averaged model and original model. The results show the orbit angular momentum and mechanical energy can be controlled independently, and the operating principles of tether reactions are derived for special modification of orbit shape.

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equibria or the stable periodic orbits will not escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric froce, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we clacluated the strengths of different forces., Comment: 36 pages, 8 figures

We study the orbits and manifolds near the equilibrium points of a rotating asteroid. The linearised equations of motion relative to the equilibrium points in the gravitational field of a rotating asteroid, the characteristic equation and the stable conditions of the equilibrium points are derived and discussed. First, a new metric is presented to link the orbit and the geodesic of the smooth manifold. Then, using the eigenvalues of the characteristic equation, the equilibrium points are classified into 8 cases. A theorem is presented and proved to describe the structure of the submanifold as well as the stable and unstable behaviours of a massless particle near the equilibrium points. The linearly stable, the non-resonant unstable, and the resonant equilibrium points are discussed. There are three families of periodic orbits and four families of quasi-periodic orbits near the linearly stable equilibrium point. For the non-resonant unstable equilibrium points, there are four cases; for the periodic orbit and the quasi-periodic orbit, the structures of the submanifold and the subspace near the equilibrium points are studied for each case. For the resonant equilibrium points, the dimension of the resonant manifold is greater than four, and we find at least 1 family of periodic orbits near the resonant equilibrium points. Besides, this theory is applied to asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia, and 6489 Golevka.

Recently, Near Earth Objects (NEOs) have been attracting great attention, and thousands of NEOs have been found to date. This paper examines the NEOs' orbital dynamics using the framework of an accurate solar system model and a Sun-Earth-NEO three-body system when the NEOs are close to Earth to search for NEOs with low-energy orbits. It is possible for such an NEO to be temporarily captured by Earth; its orbit would thereby be changed and it would become an Earth-orbiting object after a small increase in its velocity. From the point of view of the Sun-Earth-NEO restricted three-body system, it is possible for an NEO whose Jacobian constant is slightly lower than C1 and higher than C3 to be temporarily captured by Earth. When such an NEO approaches Earth, it is possible to change its orbit energy to close up the zero velocity surface of the three-body system at point L1 and make the NEO become a small satellite of the Earth. Some such NEOs were found; the best example only required a 410m/s increase in velocity., Comment: 22 pages, 6 figures, accepted for publication in Research in Astronomy and Astrophysics (Chinese Journal of Astronomy and Astrophysics)

Near Earth Asteroids have a possibility of impacting with the Earth and always have a thread on the Earth. This paper proposes a way of changing the trajectory of the asteroid to avoid the impaction. Solar sail evolving in a H-reversal trajectory is utilized for asteroid deflection. Firstly, the dynamics of solar sail and the characteristics of the H-reversal trajectory are analyzed. Then, the attitude of the solar sail is optimized to guide the sail to impact with the object asteroid along a H-reversal trajectory. The impact velocity depends on two important parameters: the minimum solar distance along the trajectory and lightness number. A larger lightness number and a smaller solar distance lead to a higher impact velocity. Finally, the deflection capability of a solar sail impacting with the asteroid along the H-reversal is discussed. The results show that a 10 kg solar sail with a lead-time of one year can move Apophis out of a 600-m keyhole area in 2029 to eliminate the possibility of its resonant return in 2036.

Advanced solar sailing has been an increasingly attractive propulsion system for highly non-Keplerian orbits. Three new applications of the orbital angular momentum reversal (H-reversal) trajectories using solar sails are presented in this paper: space observation, heliocentric orbit transfer, and collision orbits with asteroids. A theoretical proof for the existence of double H-reversal trajectories (referred to as 'H2RTs') is given, and the characteristics of the H2RTs are introduced before the discussion of the mission applications. A new family of H2RTs was obtained using a 3D dynamic model of the two-body frame. In a time-optimal control model, the minimum period H2RTs both inside and outside the ecliptic plane were examined using an ideal solar sail. Due to the quasi-heliostationary property at its two symmetrical aphelia, the H2RTs were deemed suitable for space observation. For the second application, the heliocentric transfer orbit was able to function as the time-optimal H-reversal trajectory, as its perihelion velocity is circular or elliptic velocity. Such a transfer orbit can place the sailcraft into a clockwise orbit in the ecliptic plane, with a high inclination or displacement above or below the Sun. The third application of the H-reversal trajectory was simulated impacting an asteroid passing near Earth in a head-on collision. The collision point can be designed through selecting different perihelia or different launch windows. Sample orbits of each application were presented through numerical simulation. The results can serve as a reference for theoretical research and engineering design., Comment: 19 pages, 14 figures, 2 tables. RAA 2011

In this paper we establish the local and global well-posedness of the real valued fifth order Kadomstev-Petviashvili I equation in the anisotropic Sobolev spaces with nonnegative indices. In particular, our local well-posedness improves Saut-Tzvetkov's one and our global well-posedness gives an affirmative answer to Saut-Tzvetkov's $L^2$-data conjecture., Comment: 17pages

The Kawahara and modified Kawahara equations are fifth-order KdV type equations and have been derived to model many physical phenomena such as gravity-capillary waves and magneto-sound propagation in plasmas. This paper establishes the local well-posedness of the initial-value problem for Kawahara equation in $H^s({\mathbf R})$ with $s>-\frac74$ and the local well-posedness for the modified Kawahara equation in $H^s({\mathbf R})$ with $s\ge-\frac14$. To prove these results, we derive a fundamental estimate on dyadic blocks for the Kawahara equation through the $[k; Z]$ multiplier norm method of Tao \cite{Tao2001} and use this to obtain new bilinear and trilinear estimates in suitable Bourgain spaces., Comment: 17pages

We consider the fifth order Kadomtsev-Petviashvili I (KP-I) equation as $\partial_tu+\alpha\partial_x^3u+\partial^5_xu+\partial_x^{-1}\partial_y^2u+uu_x=0,$ while $\alpha\in \mathbb{R}$. We introduce an interpolated energy space $E_s$ to consider the well-posedeness of the initial value problem (IVP) of the fifth order KP-I equation. We obtain the local well-posedness of IVP of the fifth order KP-I equation in $E_s$ for $0

In this paper we consider some dissipative versions of the modified Korteweg de Vries equation $u_t+u_{xxx}+|D_x|^{\alpha}u+u^2u_x=0$ with $0<\alpha\leq 3$. We prove some well-posedness results on the associated Cauchy problem in the Sobolev spaces $H^s({\Bbb R})$ for $s>1/4-\alpha/4$ on the basis of the $[k; Z]-$multiplier norm estimate obtained by Tao in \cite{Tao} for KdV equation., Comment: 20 pages