{"id":1021,"date":"2026-01-08T15:08:04","date_gmt":"2026-01-08T07:08:04","guid":{"rendered":"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/?page_id=1021"},"modified":"2026-01-08T15:08:04","modified_gmt":"2026-01-08T07:08:04","slug":"research","status":"publish","type":"page","link":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/research\/","title":{"rendered":"Research"},"content":{"rendered":"\n<ul class=\"link\"><li><a href=\"#GW\">Gravitational wave detection with Gaia astrometry<\/a><\/li><li><a href=\"#tianyu\">Tianyu project (\u5929\u8bed\u8ba1\u5212)<\/a><\/li><li><a href=\"#dark\">Detect black holes and compact objects and other dark companions using Gaia and other data<\/a><\/li><li><a href=\"#gaia\">Detect cold Jupiters through combined analyses of RV and Gaia-Hipparcos data<\/a><\/li><li><a href=\"#auto\">Highly automatical detection of exoplanets<\/a><\/li><li><a href=\"#pexo\">Precise exoplanetology<\/a><\/li><li><a href=\"#earth\">Develop techniques to detect Earth analogs<\/a><\/li><li><a href=\"#binary\">Stability of wide binaries<\/a><\/li><li><a href=\"#OC\">Formation and evolution of the Oort cloud in light of Sedna-like objects<\/a><\/li><li><a href=\"#habitability\">Architecture habitability of exoplanets<\/a><\/li><li><a href=\"#exocomet\">Find exocomets in the Solar System<\/a><\/li><li><a href=\"#encounter\">Find the imprint of culprit stars in long period comets<\/a><\/li><li><a href=\"#climate\">Quantum analog of terrestrial climate<a><\/li><li><a href=\"#impact\">The history of comet impact rate<\/a><\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"GW\">Gravitational wave detection with Gaia astrometry<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"624\" height=\"340\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72471.png\" alt=\"\" class=\"wp-image-922\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72471.png 624w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72471-300x163.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72471-500x272.png 500w\" sizes=\"auto, (max-width: 624px) 100vw, 624px\" \/><figcaption>\u8be5\u56fe\u6539\u81eaEllis et al, 2024, Physical Review D, Volume 109, Issue 2<\/figcaption><\/figure>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"936\" height=\"204\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72472.png\" alt=\"\" class=\"wp-image-923\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72472.png 936w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72472-300x65.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72472-768x167.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2025\/09\/\u56fe\u72472-500x109.png 500w\" sizes=\"auto, (max-width: 936px) 100vw, 936px\" \/><\/figure>\n\n\n\n<p>Since LIGO\/Virgo\u2019s first detection, <strong>hundreds<\/strong> of stellar-mass BH mergers have been observed, but supermassive BH binaries (SMBBHs) remain elusive. Pulsar timing arrays (PTAs) recently reported \u223c4\u03c3 evidence for a nano-hertz GW background via the Hellings\u2013Downs correlation, with strain \u223c10\u207b\u00b9\u2075. We propose a complementary <strong>astrometric<\/strong> search using <strong>Gaia<\/strong>: compare DR2 vs. DR3 to detect stellar accelerations induced by SMBBH GWs. With ~10\u2078 stars at ~10 \u03bcas precision, Gaia could sense ~1 nas deflections\u2014again corresponding to strain ~10\u207b\u00b9\u2075, comparable to PTA sensitivity. Unlike PTAs, Gaia offers uniform all-sky sampling and sensitivity that remains effective at tens of nHz. This project is funded by China\u2019s National Key R&amp;D Program (MOST) in collaboration with Shanghai Astronomical Observatory and Shenzhen Technology University.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"tianyu\">Tianyu project (\u5929\u8bed\u8ba1\u5212)<\/h2>\n\n\n\n<figure class=\"wp-block-gallery columns-3 is-cropped wp-block-gallery-1 is-layout-flex wp-block-gallery-is-layout-flex\"><ul class=\"blocks-gallery-grid\"><li class=\"blocks-gallery-item\"><figure><img loading=\"lazy\" decoding=\"async\" width=\"680\" height=\"886\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/tianyu.png\" alt=\"\" data-id=\"608\" data-full-url=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/tianyu.png\" data-link=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/?attachment_id=608\" class=\"wp-image-608\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/tianyu.png 680w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/tianyu-230x300.png 230w\" sizes=\"auto, (max-width: 680px) 100vw, 680px\" \/><figcaption class=\"blocks-gallery-item__caption\">Tianyu telecope <\/figcaption><\/figure><\/li><li class=\"blocks-gallery-item\"><figure><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"570\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map-1024x570.png\" alt=\"\" data-id=\"606\" data-full-url=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map.png\" data-link=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/?attachment_id=606\" class=\"wp-image-606\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map-1024x570.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map-300x167.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map-768x427.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map-500x278.png 500w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/survey_map.png 1312w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"blocks-gallery-item__caption\">Tianyu survey footprint<\/figcaption><\/figure><\/li><li class=\"blocks-gallery-item\"><figure><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"633\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit-1024x633.png\" alt=\"\" data-id=\"603\" data-full-url=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit.png\" data-link=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/?attachment_id=603\" class=\"wp-image-603\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit-1024x633.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit-300x185.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit-768x474.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit-486x300.png 486w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/detection_limit.png 1036w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"blocks-gallery-item__caption\">Tianyu detection limit<\/figcaption><\/figure><\/li><li class=\"blocks-gallery-item\"><figure><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"692\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion-1024x692.png\" alt=\"\" data-id=\"601\" data-full-url=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion.png\" data-link=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/?attachment_id=601\" class=\"wp-image-601\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion-1024x692.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion-300x203.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion-768x519.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion-444x300.png 444w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/comparion.png 1060w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"blocks-gallery-item__caption\">Tianyu compared with other surveys<\/figcaption><\/figure><\/li><\/ul><figcaption class=\"blocks-gallery-caption\">-Tianyu is aimed at detecting the second solar system including cold Jupiters and Earth-like planets as well as searching for time-domain signals such as supernovae, TDEs, GWEMs, FRB\/GRB optical counterparts, variable stars, solar system minor bodies, and Be stars. <\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\">Detection of Dark Companions by periodogram analyses of Gaia catalog data<\/h3>\n\n\n\n<p>Tianyu project is funded by the JiaoDa 2030 B program and the Tianyu telescope will start its scientific observation at the Lenghu Saishiteng Mountain (&gt;4000m altitude) in 2026. It will produce about 3PB raw data per year and will produce real-time alerts and data reduction with a time delay less than 1min. While Tianyu is mainly to detect long period transiting planets, it will also conduct time-domain science. <\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li><strong>F. Feng<\/strong>*, Y. Rui*, et al. 2024,&nbsp;<a href=\"https:\/\/ui.adsabs.harvard.edu\/abs\/2024arXiv240407149F\/abstract\"><em>Tianyu: search for the second solar system and explore the dynamic universe<\/em>,<\/a> accepted by Acta Astronomica Sinica<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"dark\">Detect black holes and compact objects and other dark companions using Gaia and other data<\/h2>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter size-large\"><img decoding=\"async\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/dark-1-1024x836.png\" alt=\"\" class=\"wp-image-597\" \/><figcaption>Top left: simulate the Gaia epoch data and refit using linear astrometric model to model Gaia DRs; top right: the principle of modeling Gaia catalog data using GOST; bottom left: black holes in binaries (El-Badry et al. 2023); bottom right: Gaia BH2 recovered by periodogram anlaysis of Gaia DRs. <\/figcaption><\/figure><\/div>\n\n\n\n<h3 class=\"wp-block-heading\">Detection of Dark Companions by periodogram analyses of Gaia catalog data<\/h3>\n\n\n\n<p>Although Gaia has released three data releases, it will not released its epoch data (or time series data) until Gaia DR4. To use Gaia DRs to detect dark companions such as black holes, neutron stars, wide dwarfs, exoplanets and other compact objects, we use GOST to simulate the Gaia epoch data and refit the synthetic epoch data using the 5-parameter linear astrometric model. By modeling the Gaia DRs in this way and calculating the chi2 or Bayes Factor for a sample of periods, we calculate the so-called astrometric Bayes Factor Periodogram (ABFP) to find periodic signals in Gaia DRs. As shown in the above figure, we are able to cover the astrometric signal for Gaia BH2, demonstrating the potential of ABFP in finding millions of dark companions and binaries. <\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li><strong>F. Feng<\/strong>, Y. Rui, Y. Xuan, H. Jones 2024,&nbsp;<em><a rel=\"noreferrer noopener\" href=\"https:\/\/ui.adsabs.harvard.edu\/abs\/2024ApJS..271...50F\/abstract\" target=\"_blank\">Modeling and Calibration of Gaia, Hipparcos, and Tycho-2 Astrometric Data for the Detection of Dark Companions<\/a><\/em>, ApJS, 271(2), p.50.<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"gaia\">Detect cold Jupiters through combined analyses of RV and Gaia-Hipparcos data<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"769\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3-1024x769.png\" alt=\"\" class=\"wp-image-208\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3-1024x769.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3-300x225.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3-768x577.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3-400x300.png 400w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-3.png 1524w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>One-planet fit to the radial velocity (top left), the proper motion (bottom left) and position (bottom right) difference between Gaia DR2 and Hipparcos catalog data for eps Indi A.  <\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\">Synergy between RV and Gaia astrometry<\/h3>\n\n\n\n<p>Since Gaia and Hipparcos measure precise positions of stars over a 24 year time span, the difference in their positions and proper motions for a given star can be used to constrain the stellar acceleration due to large planets or brown dwarfs. Our modeling of both RV and Gaia-Hipparcos offsets enable us to discover the nearest Jupiter analog, eps Indi A b. This approach will be applied to thousands of nearby stars to find a large sample of cold Jupiters and brown dwarfs. <\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li><strong>F. Feng<\/strong>, G. Anglada-Escude, M. Tuomi, H. R. A. Jones, J. Chaname, P. R. Butler, M. Janson 2019,&nbsp;<a rel=\"noreferrer noopener\" href=\"https:\/\/www.google.com\/url?q=https%3A%2F%2Facademic.oup.com%2Fmnras%2Farticle-abstract%2F490%2F4%2F5002%2F5601391%3FredirectedFrom%3Dfulltext&amp;sa=D&amp;sntz=1&amp;usg=AFQjCNFICIqdOB57i79EBG84vPnBBbq_XQ\" target=\"_blank\"><em>Detection of the nearest Jupiter analog in radial velocity and astrometry data<\/em><\/a>, MNRAS, 490, 5002<\/li><li><strong>Fabo Feng<\/strong>, R Paul Butler, et al., 2021.&nbsp;<a href=\"https:\/\/academic.oup.com\/mnras\/article-abstract\/507\/2\/2856\/6338119?redirectedFrom=fulltext\">Optimized modelling of Gaia\u2013Hipparcos astrometry for the detection of the smallest cold Jupiter and confirmation of seven low-mass companions<\/a>.&nbsp;MNRAS, Issue 2, pp.2856-2868<\/li><li><strong>F. Feng<\/strong>, et al. 2023, <a href=\"https:\/\/ui.adsabs.harvard.edu\/abs\/2023MNRAS.525..607F\/exportcitation\">Revised orbits of the two nearest Jupiters<\/a>, MNRAS, 525, 607<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"auto\">Highly automatical detection of exoplanets<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"637\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-1024x637.png\" alt=\"\" class=\"wp-image-207\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-1024x637.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-300x186.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-768x477.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-1536x955.png 1536w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2-483x300.png 483w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2.png 1908w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>Agatha is highly automized to detect planets efficiently and reliably.<\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\">The automated Agatha includes the following functions.<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>Tempered MCMC explore the global maxima of posterior <\/li><li>Cold chains explore the local maxima and constrain the signals<\/li><li>Compare models with different number of planet signals in the Bayesian framework<\/li><li>Diagnose signals using the Bayes factor periodograms (BFPs) for activity indices, different noise models and different instruments to check the signal consistency<\/li><li>Diagnose signals using the moving periodogram to check consistency over time<\/li><\/ul>\n\n\n\n<h3 class=\"wp-block-heading\">RELEVANT PUBLICATIONS:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li><strong>Feng, F.<\/strong>, Shectman, S.A., Clement, M.S., Vogt, S.S., Tuomi, M., Teske, J.K., Burt, J., Crane, J.D., Holden, B., Wang, S.X. and Thompson, I.B., 2020. Search for Nearby Earth Analogs. III. Detection of 10 New Planets, 3 Planet Candidates, and Confirmation of 3 Planets around 11 Nearby M Dwarfs.&nbsp;<em>The Astrophysical Journal Supplement Series<\/em>,&nbsp;<em>250<\/em>(2), p.29.<\/li><li><strong>Feng, F.<\/strong>, Butler, R.P., Shectman, S.A., Crane, J.D., Vogt, S., Chambers, J., Jones, H.R., Wang, S.X., Teske, J.K., Burt, J. and Diaz, M.R., 2020. Search for Nearby Earth Analogs. II. detection of five new planets, eight planet candidates, and confirmation of three planets around nine nearby M dwarfs.&nbsp;<em>The Astrophysical Journal Supplement Series<\/em>,&nbsp;<em>246<\/em>(1), p.11.<\/li><li><strong>F. Feng<\/strong>, J. D. Crane, S. X. Wang, J. K. Teske, S. A. Shectman, M. R. D\u00edaz, I. B. Thompson, H. R. A. Jones, R. Paul Butler 2019,&nbsp;<a rel=\"noreferrer noopener\" href=\"https:\/\/www.google.com\/url?q=https%3A%2F%2Fiopscience.iop.org%2Farticle%2F10.3847%2F1538-4365%2Fab1b16&amp;sa=D&amp;sntz=1&amp;usg=AFQjCNFBbH2hov5jhrulqcTTM0ubiFZYRA\" target=\"_blank\">S<em>earch for Nearby Earth Analogs. I. 15 Planet Candidates Found in PFS Data<\/em><\/a><em>,<\/em>&nbsp;ApJS, 242, 2<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"pexo\">Precise exoplanetology<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"655\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-1024x655.png\" alt=\"\" class=\"wp-image-204\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-1024x655.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-300x192.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-768x491.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-1536x982.png 1536w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-2048x1310.png 2048w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2021\/01\/image-469x300.png 469w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>PEXO is designed to implement as most high order effects as possible. It can model timing to nanosecond, astrometry to microarcsecond, and radial velocity to cm\/s precision. <\/figcaption><\/figure>\n\n\n\n<figure class=\"wp-block-image\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"571\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-1024x571.png\" alt=\"\" class=\"wp-image-630\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-1024x571.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-300x167.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-768x428.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-1536x857.png 1536w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-2048x1142.png 2048w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2024\/06\/pexo_github-500x279.png 500w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>PEXO is a packaged designed to model astrometry to microarcsecond, RV to cm\/s, and timing to nanosecond precision. The relevant paper was published on ApJS. The R and python interface are all available at GitHub: <a href=\"https:\/\/github.com\/phillippro\/pexo\">https:\/\/github.com\/phillippro\/pexo<\/a> and binder: <a href=\"https:\/\/mybinder.org\/v2\/gh\/phillippro\/pexo\/HEAD?filepath=demos%2FIntroduction.ipynb\">https:\/\/mybinder.org\/v2\/gh\/phillippro\/pexo\/HEAD?filepath=demos%2FIntroduction.ipynb<\/a><\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\">RELEVANT PUBLICATIONS:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li><strong>F. Feng<\/strong>, M. Lisogorskyi, H. R. A. Jones, S. M. Kopeikin, R. P. Butler, G. Anglada-Escude, A. P. Boss 2019, PEXO: a global modeling framework for nanosecond timing, microsecond astrometry, and \u03bcm\/s radial velocities, ApJS, 244, 39<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"earth\">Develop techniques to detect Earth analogs<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"858\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti-1024x858.png\" alt=\"\" class=\"wp-image-98\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti-1024x858.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti-300x251.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti-768x643.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti-358x300.png 358w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/tau_ceti.png 1280w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>Four super-Earths orbiting around tau Ceti corresponding to the weakest planet-induced radial velocity signals ever been detected. <\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>Goldilocks principle in noise modeling<\/strong><\/h3>\n\n\n\n<p>During postdoctoral research at the University of Hertfordshire, I have applied the Bayesian reference in the analysis of radial velocity data to detect exoplanets. I have introduced the so-called \u00e2\u20ac\u0153Goldilocks principle\u00e2\u20ac\u009d in the noise modeling of radial velocity data to avoid false positives and false negatives (Feng et al. 2016). According to my research, the moving average model is optimal in modeling correlated noise.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>Disentangling periodic signals from correlated noise in the Agatha framework<\/strong><\/h3>\n\n\n\n<p>I have developed the red noise periodograms called Bayes factor periodogram (BFP) and marginalized likelihood periodogram (MLP) based on likelihood optimization and marginal- ization. These periodograms are combined to form a framework, which is implemented in the Agatha web application (http:\/\/www.agatha.herts.ac.uk). Agatha has been found to out- perform previous periodograms in terms of removing correlated noise, selecting noise models, and testing consistency of signals (Feng et al. 2017b). It is also one of the data analysis tools for the Red Dots Campaign.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>Detection of the weakest signals in radial velocity data<\/strong><\/h3>\n\n\n\n<p>I have also pioneered a new noise proxy called \u00e2\u20ac\u0153differential radial velocities\u00e2\u20ac\u009d to model the wavelength-dependent noise in the radial velocity data. This proxy is able to weight a posteriori the radial velocity data measured within different wavelength ranges. The application of this proxy and the other statistical methods I have developed in the analysis of the radial velocity data of HD 20794 and \u00cf\u201e Ceti leads to detections of a few Earth-size planets corresponding to radial velocity variation as low as 0.3 m\/s, the weakest signal ever detected using the radial velocity method (Feng et al. 2017c,d). The work on \u00cf\u201e Ceti has been covered by media such as CNN, BBC, Daily Mail, Sky News, Gizmodo, Sky &amp; Telescope, etc.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng, M. Tuomi, H. R. A. Jones, J. Barnes, G. Anglada-Escude, S. Vogt &amp; R. P. Butler, 2017, Color difference makes a difference: four planet candidates around tau Ceti, AJ, 154, 135<\/li><li>F. Feng, M. Tuomi &amp; H. R. A. Jones, 2017, Evidence for at least three planet candidates orbiting HD20794, A&amp;A, 605, A103<\/li><li>F. Feng, M. Tuomi &amp; H. R. A. Jones, 2017, Agatha: Disentangling periodic signal from correlated noise in a periodogram framework, MNRAS, 470, 4794, the corresponding online application is available at http:\/\/www.agatha.herts.ac.uk<\/li><li>F. Feng, M. Tuomi, H. R. A. Jones, R. P. Butler &amp; S. Vogt, 2016, A Goldilocks principle for modeling radial velocity noise, MNRAS, 461, 2440<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"binary\">Stability of wide binaries<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"448\" height=\"492\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/proxima_capture.png\" alt=\"\" class=\"wp-image-104\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/proxima_capture.png 448w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/proxima_capture-273x300.png 273w\" sizes=\"auto, (max-width: 448px) 100vw, 448px\" \/><figcaption><em>Instability of Proxima clones under the pertubation from the Galactic tide and stellar encounters based on simulations from -5 Gyr to +5 Gyr. The black and blue lines denote the orbital elements of clones which were\/will be ejected through combined perturbations of stellar encounters and the Galactic tide. The red and green lines denote the orbits of stable clones.<\/em><\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>Capture scenario of the Proxima-alpha Centauri system<\/strong><\/h3>\n\n\n\n<p>I am currently investigating the influence of stellar encounters on the dynamical evolution of wide binaries. I have found that Proxima was probably captured by alpha Centauri based on simulations of the Proxima\u00e2\u20ac\u2122s orbits under the perturbations from the Galactic tide and stellar encounters over \u00c2\u00b15 Gyr. This hypothesis is also supported by the metallicity difference between Proxima and alpha Centauri A and B (Feng &amp; Jones 2017b). If the capture scenario is true, Proxima b (Anglada-Escude et al. 2016) is more likely to be in the habitable zone for a long time due to a shorter time scale of the interaction between Proxima b and alpha Centauri. This work has been covered by IFLScience, New Scientist, Gizmodo, etc.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><strong>A new classification of wide binaries and multiples<\/strong><\/h3>\n\n\n\n<p>In this paper, I argue that Fomalhaut C is a so-called gravitational pair\u009d of Foma-haut A and B. In other words, Fomalhaut C is in a meta-stable state such that it becomes bound and unbound repeatedly with respect to Fomalhaut A and B. Like the Cooper pair mechanism in superconductors, this phenomena only appears once the orbital energy of a com- ponent becomes comparable with the energy fluctuations caused by the environment. Gaia will provide high precision astrometry for millions of stars which provide six dimensional initial conditions of stars with velocity errors less than 1 km\/s. Numerical integration of the orbits of these stars along with abundance analyses can enable discoveries of coeval and non-coeval gravitational pairs in bound and unbound states. This project may also provide candidates of gravitational pairs of the solar system and test the Nemesis hypothesis (Whitmire &amp; Jackson 1984).<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; H. R. A. Jones, 2018, Understanding Fomalhaut as a Cooper pair, MNRAS, 474, 4412<\/li><li>F. Feng &amp; H. R. A. Jones, 2018, Was Proxima captured by alpha Centauri A and B?, MNRAS, 473, 3185<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"OC\">Formation and evolution of the Oort cloud in light of Sedna-like objects<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"643\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide-1024x643.png\" alt=\"\" class=\"wp-image-107\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide-1024x643.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide-300x188.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide-768x482.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide-478x300.png 478w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oort_cloud_tide.png 1280w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption>Asymmetric tidal perturbation on a planetary system<\/figcaption><\/figure>\n\n\n\n<p>The orbits of Sedna and 2012 VP112 have generated new interest in the formation and evolution of the Oort Cloud. The perihelia of Sedna and 2012 VP113 are 76 and 81 AU, and their eccentricities are 0.86 and 0.70, respectively. Considering their high eccentricities, they could not form in their current locations, otherwise, they would be disrupted by planetesimals with large relative velocities. Thus they may have migrated to their current eccentric orbits through interactions with one or a few massive objects, which could be an unseen planet or a solar companion beyond the Kuiper belt, a single star, or stars in the open cluster in which the Sun formed. All these scenarios more or less assume the existence of an unseen solar companion or stellar encounters. However, the recent result of WISE survey seems to exclude most of the combinations of the mass and orbital elements of an unseen solar companion. In addition, the stochastic stellar encounters in this scenario may not properly explain the concentration of the arguments of perihelia of Sedna-like objects at omega=0.<\/p>\n\n\n\n<p>Instead of invoking an unknown planet or stellar encounter, I will model the formation of the Sedna-like bodies based on the solar apex motion, which could induce a preferred plane where the tidal force imposed on the Oort cloud is stronger than that far from the plane. This project aims to model the formation of the inner Oort cloud and the extended scattered disc, and test whether the perihelion concentration can be caused by the solar apex motion.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publication<\/h3>\n\n\n\n<p>F. Feng &amp; C.A.L. Bailer-Jones, 2014, Exploring the role of the Sun&#8217;s motion in terrestrial comet impacts, MNRAS, 442, 3653<\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"habitability\">Architecture habitability of exoplanets<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"563\" height=\"583\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/habitability.png\" alt=\"\" class=\"wp-image-110\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/habitability.png 563w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/habitability-290x300.png 290w\" sizes=\"auto, (max-width: 563px) 100vw, 563px\" \/><figcaption>Hierarchical habitability assessment<\/figcaption><\/figure>\n\n\n\n<p>The concept of the habitable zone (HZ) is useful in selecting candidates of habitable planets for follow up observations and studies. Exoplanets in the HZ may not be habitable due to many factors such as a lack of water or unstable orbits. These factors are more or less related to the architecture of planetary systems that host terrestrial planets in their HZ. Thus I intend to build a model to predict which architectures are habitable in the sense of allowing the formation and dynamical stability of terrestrial planets in the HZ. Based on simulations of the dynamics of the terrestrial planets under the perturbations from the whole system, I will also build models of the orbital variations, such as eccentricity, obliquity and precession, and asteroid\/comet impacts, and assess the influence of these factors on the climate of terrestrial planets. Finally, I will piece together these models to find the main characteristics of habitable architectures, and combine my architecture habitability models with the HZ to better predict planetary habitability.<\/p>\n\n\n\n<p>Furthermore, I would like to assess the habitability of exoplanets in using a hierarchical model. The calculation of the HZ assumes that the incoming stellar radiation is the main factor that determines whether a rocky planet can support liquid water on its surface. Although refined by accounting for many other factors, such as the atmosphere, cloud cover, the albedo and the orbital eccentricity of the rocky planet, the HZ only provides a basic selection of habitable planets. A more self-consistent and comprehensive way to find habitable worlds is to model planetary habitability hierarchically. Although many physical and architectural properties of rocky planets and their host systems are unknown, these uncertainties can be included in the priors of a Bayesian hierarchical model. Specifically, the hierarchical structure of this model is: (1) the galaxy evolution determines whether the chemical composition in a galaxy can support life, which defines the concept of &#8220;cosmic habitable age&#8221; or CHA; (2) the habitability of a planet may also be influenced by its interstellar environment which has components such as star formation, stellar encounters and cosmic rays (see Bailer-Jones 2009 and Feng &amp; Bailer-Jones 2013 for reviews). This defines the &#8220;Galactic habitable zone&#8221; or GHZ; (3) if a planet is located in the GHZ, its habitability is further determined by the architecture of its host planetary system; (4) in addition to these external factors, the internal structure and composition of the planet finally determines whether it is really habitable.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; C.A.L. Bailer-Jones, 2014, Exploring the role of the Sun&#8217;s motion in terrestrial comet impacts, MNRAS, 442, 3653<\/li><li>F. Feng &amp; C.A.L. Bailer-Jones, 2015, Pleistocene deglaciations paced by variations of the Earth\u00e2\u20ac\u2122s orbit, Quaternary Science Review, 122, 166<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"exocomet\">Find exocomets in the Solar System<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/oumuamua.png\" alt=\"\" class=\"wp-image-111\" width=\"517\" height=\"289\" \/><figcaption>&#8216;Oumuamua, the first interstellar visitor to the Solar System<\/figcaption><\/figure>\n\n\n\n<p>I have studied the first insterstellar object &#8216;Omuamua. I simulate the trajectories of this object and nearby stars using the TGAS catalog and find that it is likely a young objects ejected from a star in the so-called &#8220;Local Association&#8221;. I have found the variation of the rotation period of this object which lead to the proposal that this object is tumbling and thus further support that it is probably a young object. I have also written an article on Conversation to introduce my discoveries.<\/p>\n\n\n\n<p>There are possibly many extra-solar comets in the solar system which were captured by the Sun from the interstellar space (e.g. Machholz 1 and Hyakutake C\/1996 B2). I am interested in finding how many comets were captured by the Sun from stellar encounters. I will first simulate the interaction between the Sun (with its Oort cloud) and fly-bys (with its OC). By adopting different parameters of flybys, I can build a semi-analytical model to estimate how many exocomets can be captured by the Sun given a certain type of flybys. Then I am able to estimate the total population of exocomets in the solar system over its history. Finally, I may find whether the Sun have captured comets from know stellar encounters. If there are such exocomets, I will further assess which comets in the solar system could be the best candidates.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; C.A.L. Bailer-Jones, 2014, Exploring the role of the Sun&#8217;s motion in terrestrial comet impacts, MNRAS,442, 3653<\/li><li>F. Feng &amp; H. R. A. Jones, 2018, `Oumuamua as a messenger from the Local Association, ApJL, 852, L27<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"encounter\">Find the imprint of culprit stars in long period comets<\/h2>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-flow wp-block-group-is-layout-flow\">\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-flow wp-block-group-is-layout-flow\">\n<figure class=\"wp-block-image size-large is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/stellar_encounter.png\" alt=\"\" class=\"wp-image-112\" width=\"433\" height=\"238\" \/><\/figure>\n\n\n\n<figure class=\"wp-block-image size-large is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/encounter2.png\" alt=\"\" class=\"wp-image-113\" width=\"437\" height=\"236\" \/><figcaption><em>An Aitoff projection in Galactic coordinates of the aphelia\/perihelia of LPCs\/encounters for the reliable sequence of past encounters (left panel) and HIP 103738 (Gamma Microscopii, right panel). In the left panel, the blue crosses show the perihelion of the encounters. The size of the cross is proportional to the number of LPCs the encounter. The red crosses show the aphelia of dynamically new LPCs in the simulation and the black points show the aphelia of observed dynamically new LPCs. The black star symbol indicates the solar apex, and the dashed line shows the preferred tidal plane, the normal to which is the solar apex. The right panel is similar, but is only for HIP 103738, and does not show the observed LPCs and additionally shows the aphelia of LC LPCs in the simulation (yellow open circles). The blue curve is actually a sequence of points showing the encounter trajectory when its pericentric distance is less than 4 pc. Each point gives the star\u2019s position at a equally-separated time interval of 0.001 Myr.<\/em><\/figcaption><\/figure>\n<\/div><\/div>\n<\/div><\/div>\n\n\n\n<p>Long period comets (LPCs), which have high eccentric orbits, are delivered from the Oort cloud by the tidal forces from the Galatic tide (smooth potential) and stellar encounters (discrete\/stochastic potential). These tidal forces can change the perihelia of Oort cloud comets and make them discernable when their perihelia are reduced to less than 5 AU. The studies of LPCs are important in understanding the bombardment, geological, biological and climatic history on terrestrial planets\/satellites. In addition, the investigation of the origin of LPCs will advance our understanding of the birth environment of the solar system by constraining the population of the inner Oort cloud. I plan to build a model to assess the connection between LPCs and stellar encounters (or culprit stars) based on simulations of the Oort cloud under the perturbations of the Galactic tide and the sample of stellar encounters found by Bailer-jones 2015. This work will provide a method to find the imprint of a culprit star in the currently observed LPCs, and may first identify one or a few culprits from these LPCs. It will also be a prototype of finding culprit stars using future Gaia sample of stellar encounters.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; C.A.L. Bailer-Jones, 2014, Exploring the role of the Sun&#8217;s motion in terrestrial comet impacts, MNRAS,442, 3653<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"climate\">Quantum analog of terrestrial climate<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"576\" height=\"261\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/milankovitch.png\" alt=\"\" class=\"wp-image-114\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/milankovitch.png 576w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/milankovitch-300x136.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/milankovitch-500x227.png 500w\" sizes=\"auto, (max-width: 576px) 100vw, 576px\" \/><figcaption><em>Climate variations over the Pleistocene. The present day is at time zero on the right. The \u03b4<sup>18<\/sup>O record (lower solid line) stacked by Lisiecki and Raymo (2005) is compared with the daily-averaged insolation at the summer solstice at 65<sup>\u25e6<\/sup>N, Q<sup>\u0304day<\/sup>65<sup>\u25e6<\/sup>N(upper solid line), the obliquity (dashed line), and the eccentricity (dotted line) calculated by Laskar et al. (2004). The latter two have been scaled to have a common amplitude. The grey region around \u22121000 kyr represents the MPT extending from \u22121250 kyr to \u2212770 kyr (Clark et al., 2006). The grey bar extending from \u2212423 to \u2212362 kyr represents marine isotope stage (MIS) 11. The\u03b4<sup>18 <\/sup>O variations are dominated by 41 kyr and 100 kyr cycles before and after the MPT respectively.<\/em><\/figcaption><\/figure>\n\n\n\n<p>The climate on the Earth over the last 2.5 Myr is characterized by 41 and 100 kyr glacial-interglacial cycles. Milankovitch proposed that the orbital variations (eccentricity, obliquity and precession) of the Earth can modulate the insolation (incoming solar radiation), and thus produce variations with corresponding frequencies in paleoclimatic proxies such as delta18O (a measure of the ratio between 18O and 16O. The role of different orbital elements in triggering glacial terminations has been addressed in my recent work. However, linear models of the climate response to the orbital variations (or orbital forcings) cannot explain the transition from the 41 kyr to the 100 kyr cycles occurred about 800 kyr ago (known as mid-Pleistocene transition). In addition, the sawtooth-like glacial-interglacial cycles over the last 1 Myr are found to be caused by non-linear climate responses. To solve this so-called &#8220;100-kyr problem&#8221;, climate models with various feedback mechanisms have been proposed, but only a few are successful in predicting most features in climate proxies. Common factors in these favored models are phase transitions, thresholds and bifurcations. These concepts more or less reflect a quantum-like climate system. Modeling the terrestrial climate using a perceptual analogy is not unusual in modeling climate. For example, Broecker 1991 depicted the thermohaline circulation as a \u00e2\u20ac\u02dcconveyor belt\u00e2\u20ac\u2122. With this kind of analogy approach, I plan to first identify the &#8220;transitions&#8221;, such as the transition from a fully frozen state (so-called &#8220;snowball Earth&#8221;) to a warm state, the mid-Pleistocene transition, and transitions from glaciations to deglaciations, in the paleoclimatic proxies. Assuming that the Earth can store and release certain discrete amount of energy, I intend to use &#8220;energy levels&#8221; of the Earth system to model these transitions, and to explain these energy levels based on energy balance models. For example, the ice sheets and the oceans are like an energy reservoir, which can store a certain amount of energy. On each energy level, the oceanic and ice-sheet dynamics may stabilize the climate. However, the transition from one level to another would occur if the insolation and the latent heat absorbed by the energy reservoir exceed its capacity. This quantum analogy model treats the Earth as a whole system, and thus can be easily applied to study the climate of exoplanets.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; C.A.L. Bailer-Jones, 2015, Pleistocene deglaciations paced by variations of the Earth&#8217;s orbit, Quaternary Science Review, 122, 166<\/li><\/ul>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"impact\">The history of comet impact rate<\/h2>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"771\" src=\"http:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact-1024x771.png\" alt=\"\" class=\"wp-image-115\" srcset=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact-1024x771.png 1024w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact-300x226.png 300w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact-768x578.png 768w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact-398x300.png 398w, https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-content\/uploads\/sites\/12\/2020\/12\/comet_impact.png 1182w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<p>The fossil record indicates that the biodiversity has varied a lot over the Phanerozoic eon (past 550 Myr). We try to find the extraterrestrial causes for the biodiversity variation through building astronomical models for different geological data sets. The craters on the Earth, Moon and other planets are caused by the impacts of either asteroids or comets. Considering the relatively isolated environment of asteroids, the most intensive variation of impact rate comes from the varying flux of new comets originated from the Oort cloud. The Oort cloud can be perturbed by Galactic tide and stellar encounters which are probably modulated by the solar motion. To build an end-to-end model for the comet impacting flux, we have first studied the comet flux caused by tidal force and stellar encounters separately. Including only galactic tide, 1e6 comet orbits are simulated over 1 Gyr and the new comets entering into the &#8220;loss cone&#8221; are modulated obviously by solar motion mainly through the vertical tidal force. Only close encounters can trigger comet shower which is comparable with the background flux injected by Galactic tide. In the first Gyr, we find the synergy effect between Galactic tide and encounters is not obvious. This is because the phase space corresponding to the trajectories entering into the &#8220;loss cone&#8221; are not depleted yet. We also find that the flux caused by the tide is modulated by the solar motion mainly through modulated vertical tidal force. This work will be the first to combine modulated galactic tide and constant encounter flux. It would be interesting to find the sensitivity of final comet flux to model parameters such as initial condition of the Sun.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Relevant publications:<\/h3>\n\n\n\n<ul class=\"wp-block-list\"><li>F. Feng &amp; C.A.L. Bailer-Jones, 2014, Exploring the role of the Sun&#8217;s motion in terrestrial comet impacts, MNRAS, 442, 3653<\/li><\/ul>\n\n\n\n<p><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Gravitational wave detection with Gaia a &hellip; <a href=\"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/research\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":11,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-1021","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/pages\/1021","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/users\/11"}],"replies":[{"embeddable":true,"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/comments?post=1021"}],"version-history":[{"count":1,"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/pages\/1021\/revisions"}],"predecessor-version":[{"id":1022,"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/pages\/1021\/revisions\/1022"}],"wp:attachment":[{"href":"https:\/\/web.tdli.sjtu.edu.cn\/ffeng\/wp-json\/wp\/v2\/media?parent=1021"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}