CONFERENCE PROGRAM 1 Nov – 5 Nov, 2018 HUAQIAO UNIVERSITY

DDAP10, Huaqiao University, Xiamen China, 2018

The 10th Dynamics Days Asia Pacific (DDAP10)

CONFERENCE PROGRAM

1 Nov – 5 Nov, 2018 HUAQIAO UNIVERSITY, XIAMEN，CHINA

http://ddap10.hqu.edu.cn/

DDAP10, Huaqiao University, Xiamen China, 2018

Table of Contents

I. Sponsors, Organizers and Committees ...... 3 1.1 Introduction ...... 3 1.2 Sponsors ...... 3 1.3 Financial Support ...... 4 1.4 Local Organization Committee ...... 4 1.5 Conference Working Committee ...... 4 1.6 International Advisory Committee (IAC) ...... 4 II. General Information ...... 6 III. Maps ...... 8 3.1 Campus Map ...... 8 3.2 Map for Hotels and WYX Conference Center ...... 9 3.3 Map for Airport, North Train Station and Hotels ...... 10 IV. Programs ...... 11 4.1 An Overview ...... 11 4.2 Program ...... 12 V. Speakers and Abstracts ...... 24 5.1 List of Speakers...... 24 5.1.1 Plenary Speakers ...... 24 5.1.2 Session Speakers ...... 25 5.1.3 Posters ...... 27 5.2 Abstracts ...... 29 5.2.1 Plenary Talks (PT) ...... 29 5.2.2 Session Talks (ST) ...... 43 5.2.3 Posters ...... 100 VI. Contacts ...... 123 VII. Acknowledgement ...... 124

DDAP10, Huaqiao University, Xiamen China, 2018

I. Sponsors, Organizers and Committees

1.1 Introduction

The 10th Dynamics Days Asia-Pacific (DDAP10) will be held at Xiamen Campus of Huaqiao University (Xiamen, Fujian Province, CHINA) from 2018-10-27 to 2018-11-05. The acticity of DDAP10 contains the following two parts: ①Winter School on Complex and Nonlinear Systems: The winter school is scheduled from 2018.10-28 to 2018.11.1. The school aims at a systematic series of lectures on nonlinear dynamics and complexity science. The participants mainly covers young scientists, postdocs, PhD students and master students from China. The language will be Chinese. We will invites experts as the lecturers to give talks. The scheduled number of participants will be about 100 people. ②The DDAP10 Conference: The conference is scheduled from 2018.11.1 to 2018.11.5. The participants will cover colleagues from Asia Pacific Region (including North America Region) and Europe. The number of participants will be around 500 people. The DDAP10 Conference is the dominant part of all activities of DDAP10. There will be 20 plenary talks, 100 session talks (including session invited and oral talks) and posters (no number limit) during the Conference. Dynamics Days Asia-Pacific (DDAP) is one of the four influential and regular international series of conferences on dynamical systems, including Dynamics Days USA, Dynamics Days Europe and Dynamics Days South America. The first conference (DDAP1) was sponsored by Professor Bambi Hu, dean of Physics Department, Director of Center for Nonlinear Studies, at Hong Kong Baptist University, and was launched in 1999 at HKBU. The Previous DDAP's have been held respectively in Hong Kong (DDAP1, 1999), Hangzhou (DDAP2, 2002), Singapore (DDAP3, 2004), Pohang (DDAP4, 2006), Nara (DDAP5, 2008), Sydney(DDAP6, 2010), Taipei (DDAP7, 2012), Chennai (DDAP8, 2014) and Hong Kong (DDAP9, 2016). DDAP10 is financially supported by National Natural Science Foundation of China, Funds from Fujian Province and Huaqiao University.

1.2 Sponsors

Institute of Systems, Huaqiao University College of Information Science and Engineering, Huaqiao University

DDAP10, Huaqiao University, Xiamen China, 2018

1.3 Financial Support

National Natural Science Foundation of CHINA Fujian Provincial Government Huaqiao University

1.4 Local Organization Committee

 Zhigang Zheng (Chair, Huaqiao University)  Jixiong Pu (Vice Chair, Huaqiao University)  Sumiyoshi Abe (Vice Chair, Huaqiao University)  Leihan Tang (Vice Chair, Hong Kong Baptist University)  Changsong Zhou (Vice Chair, Hong Kong Baptist University)  Dahai He (Xiamen University)  Nianbei Li (Huaqiao University)  Bihong Lin (Huaqiao University)  Qing Lin (Huaqiao University)  Jiliang Luo (Huaqiao University)  Congjie Ou (Huaqiao University)  Hong Zhao (Xiamen University)  Weidong Zou (Longyan University)

1.5 Conference Working Committee

 Jixiong Pu (Group Leaders, Huaqiao University)  Rongmei Chen (Group Leaders, Huaqiao University)  Zhigang Zheng (Group Leaders, Huaqiao University)  Hongbin Chen (Huaqiao University)  Haixiang Fu (Huaqiao University)  Wenjing Jia (Huaqiao University)  Zhifu Huang (Huaqiao University)  Baoqing Lin (Huaqiao University)  Longkun Tang (Huaqiao University)  Can Xu (Huaqiao University)  Yun Zhai (Huaqiao University)

1.6 International Advisory Committee (IAC)

Sumiyoshi Abe (Huaqiao University, China) Daniel M. Abrams (Northwestern University, USA) Giulio Casati (University of Insubria, Italy) Jinde Cao (Southeast University, China) 4

DDAP10, Huaqiao University, Xiamen China, 2018

Guanrong Chen (City University of Hong Kong, HKSAR, China) Zengru Di (Beijing Normal University, China) Celso Grebogi (The University of Aberdeen, UK) Peter Hanggi (University of Augsburg, Germany) Chin Kun Hu (Academia Sinica, Taipei) Gang Hu (Beijing Normal University, China) Liang Huang (Lanzhou University, China) Pak Ming Hui (Chinese University of Hong Kong, HKSAR, China) Kunihiko Kaneko (The University of Tokyo, Japan) Doochul Kim (Institute for Basic Science, Korea) Seunghwan Kim (Pohang University of Science and Technology, Korea) Hyunggyu Park (Korea Institute for Advanced Study, Korea) Yuri S. Kivshar (Australian National University, Australia) Yoshiki Kuramoto (Kyoto University, Japan) Jurgen Kurths (Potsdam University, PIK, Potsdam, Germany) Choy-Heng Lai (National University of Singapore, Singapore) Yingcheng Lai (Arizona State University, USA) Carlo R. Laing (Massey University, New Zealand) Baowen Li (Tongji University/University of Colorado at Boulder, China/USA) Jibin Li (Huaqiao University, China) Jie Liu (IAPCM, China) Fabio Marchesoni (Tongji University, China) Adilson E. Motter (Northwest University, USA) Qi Ouyang (Peking University, China) Zhongcan Ouyang (Chinese Academy of Sciences, China) Akardy Pikovsky (Potsdam University, Germany) Ramakrishna Ramaswamy (Jawaharlal Nehru University, India) Jianwei Shuai (Xiamen University, China) Eugene Stanley (Boston University, USA) Chao Tang (Peking University, China) Leihan Tang (Hong Kong Baptist University, HKSAR, China) Mark Timme (MPI for Dynamics and Self-Organization, Germany) Penger Tong (HKUST, HKSAR, China) Tamas Vicsek (Eötvös Loránd University (ELTE), Hungary) Bing-Hong Wang (China Univ of Science & Technology, China) Wei Wang (Nanjing University, China) K. Y. Michael Wong (HKUST, HKSAR, China) Chenxu Wu (Xiamen University, China) Jinghua Xiao (Beijing University of Post and Telecommunications, China) Lu Yu (Chinese Academy of Sciences, China) Bo Zheng (Zhejiang University, China) Haijun Zhou (Chinese Academy of Sciences, China)

DDAP10, Huaqiao University, Xiamen China, 2018

II. General Information

1. Conference Venue: The DDAP10 will be held during Nov. 1-5, 2018 in Wang-Yuan-Xing International Conference Center （ WICC, 王源兴国际会议中心） at Huaqiao University (Xiamen Campus).

2. Arrival Registration:  The reception time is 8am-8pm on Nov. 1 and each participant can register at the reception desks in his/her booked hotel (the three hotels referenced by the conference: North Bay Wanda Hotel 北海湾惠龙万达酒店, Ling Ling Hotel 灵玲大酒店 and Mingzhu Harbor Hotel 明珠海湾大酒店).  There will be Reception Buffet Dinner on Nov. 1 at each hotel for registered participants (North Bay Wanda Hotel: 6pm-9pm; Ling Ling Hotel: 6pm-8pm; Mingzhu Harbor Hotel: 6pm-8:30pm).  During the conference time between Nov. 2 and Nov. 5, participants can make on-site registration at the conference center of WICC.

3. Shuttle Buses: The departure time for shuttle bus from the hotel to WICC is 7:40am for North Bay Wanda Hotel and Ling Ling Hotel, 7:30am for Mingzhu Harbor Hotel. The shuttle bus will be arranged outside the WICC and take guests back to their hotel after the last talk scheduled on that day.

4. Duration for Talks: The time for Plenary Talks is 40 minutes, for Invited Session Talks is 25/30 minutes, for regular Contributed Session Talks is 20 minutes, and for short session talks is 15 minutes. The last 3 minutes is reserved for questions and answers.

5. Posters: Posters will be displayed in the hall of the Wang-Yuan-Xing International Conference Center （WICC）.

6. Conference Banquet: The Conference Banquet will be held on Friday, Nov. 02 at 7pm in Ling Ling Hotel. The shuttle bus will be arranged to take guests back to their own hotels after the banquet.

7. Cafeteria and Food:  The Box Lunch will be provided for every participant at the WICC. The Buffet Dinner will be provided by each Hotel at 6:30pm.  There are many Chinese and western restaurants outside the west gate of 6

DDAP10, Huaqiao University, Xiamen China, 2018

university campus. 8. Tips for the Transport on your arrival:  There are only two terminals (T3 and T4) in Xiamen Gaoqi Airport.  The distance between the Airport and the North Bay Wanda Hotel （北海湾惠龙万达酒店）is about 11km and the cost for taxi is about 30-40 RMB, between the Airport and the Mingzhu Harbor Hotel （明珠海湾大酒店） is about 12km and the cost is about 40 RMB, between the Airport and the Ling Ling Hotel （灵玲大酒店） is about 14km and the cost is about 45 RMB.  The distance between the Airport and the Conference Center (Wang-Yuan-Xing Center at Huaqiao University 华侨大学王源兴国际会议中心) is about 15km and the cost is about 50 RMB.  For participants from the Xiamen Railway Station or North Station, please take the taxi to hotels.  If you arrive in Xiamen and wish go directly to Huaqiao University, please take the BRT-1 express bus （快 1 线） to Huaqiao-University station.

DDAP10, Huaqiao University, Xiamen China, 2018

III. Maps

3.1 Campus Map

Map for the Wang-Yuan-Xing Center （王源兴国际会议中心） and the campus

DDAP10, Huaqiao University, Xiamen China, 2018

3.2 Map for Hotels and WYX Conference Center

DDAP10, Huaqiao University, Xiamen China, 2018

3.3 Map for Airport, North Train Station and Hotels

The distance between the Airport and the North Bay Wanda Hotel （北海湾惠龙万达酒店） is about 11km and the cost for taxi is about 30-40 RMB, between the Airport and the Mingzhu Harbor Hotel （明珠海湾大酒店） is about 12km and the cost is about 40 RMB, between the Airport and the Ling Ling Hotel （灵玲大酒店） is about 14km and the cost is about 45 RMB, between the Airport and the Conference Center (Wang-Yuan-Xing Center at Huaqiao University 华侨大学王源兴国际会议中心) is about 15km and the cost is about 40 RMB.

DDAP10, Huaqiao University, Xiamen China, 2018

IV. Programs

4.1 An Overview

All academic activities will be arranged in Wang-Yuan-Xing International Conference Center (WICC), Huaqiao University, Xiamen. Plenary Session A Session B* Session C Session D WICC-G204 WICC-G201 WICC-G204 WICC-G101 WICC-G103 *Note: Session B will be temporarily changed to WICC-G203 on November 3rd morning.

11.1 Desk registration (8:00-22:00)

11.2 8:30-9:00 Opening Ceremony Session-1 A B C D Banquet

8:50-9:00 Photo Taking 13:30-14:00 ST-05 ST-37 ST-51 ST-70

Plenary-1 14:00-14:30 ST-03 ST-38 ST-67 ST-71

9:00-9:40 PT-11: Yingcheng Lai 14:30-15:00 ST-07 ST-42 ST-49 ST-72

9:40-10:20 PT-09: Kunihiko Kaneko 15:00-15:30 ST-09 ST-40 ST-50 ST-92

10:20-10:40 Coffee break 15:30-16:00 Coffee break

10:40-11:20 PT-01: Ravindra Amritkar 16:00-16:30 ST-06 ST-41 ST-47 ST-89

11:20-12:00 PT-03: Sergej Flach 16:30-17:00 ST-14 ST-43 ST-53 ST-91

17:00-17:20 ST-02 ST-44 ST-59 ST-93

17:20-17:40 ST-04 ST-45 ST-48 ST-90

17:40-18:00 ST-08 ST-46 ST-52 ST-88

11.3 Session-2 A B C D Plenary-2 TBA

8:30-9:00 ST-16 ST-39 ST-57 ST-73 13:30-14:10 PT-10: Jurgen Kurths

9:00-9:30 ST-19 ST-34 ST-63 ST-81 14:10-14:50 PT-13: Adilson E. Motter

9:30-10:00 ST-13 ST-35 ST-66 ST-80 14:50-15:30 PT-18: Masaki Sano

10:00-10:20 Coffee break 15:30-16:00 Coffee break

10:20-10:45 ST-21 ST-27 ST-60 ST-82 16:00-16:40 PT-19: Zhisong Wang

10:45-11:10 ST-17 ST-94 ST-29 ST-84 16:40-17:20 PT-07: Lin I

11:10-11:30 ST-10 ST-26 ST-55 ST-79 17:20-18:00 PT-02: Eberhard Bodenschatz

11:30-11:45 ST-18 ST-33 ST-61 ST-83

11:45-12:00 ST-11 ST-36 ST-56 ST-86

11.4 Plenary-3 Free TBA

8:30-9:10 PT-05: Peter Hanggi

9:10-9:50 PT-14: Qi Ouyang

9:50-10:10 Coffee break

10:10-10:50 PT-20: Changsong Zhou

10:50-11:30 PT-16: Peter Robinson

11:30-12:10 PT-08: Hawoong Jeong

DDAP10, Huaqiao University, Xiamen China, 2018

11.5 Session-3 A B C D Plenary-4

8:30-9:00 ST-22 ST-31 ST-58 ST-77 13:30-14:10 PT-17: Miguel A.F. Sanjuan

9:00-9:30 ST-23 ST-24 ST-68 ST-78 14:10-14:50 PT-12: Carlo R. Laing

9:30-10:00 ST-25 ST-32 ST-54 ST-75 14:50-15:30 PT-06: Hiroshi H. Hasegawa

10:00-10:20 Coffee break 15:30-15:50 Coffee break

10:20-10:45 ST-20 ST-28 ST-65 ST-76 15:50-16:30 PT-15: Zhilin Qu

10:45-11:10 ST-15 ST-95 ST-96 ST-74 16:30-17:10 PT-04: Celso Grebogi

11:10-11:30 ST-01 ST-30 ST-64 ST-85 17:10-17:20 Best Presentation Award

11:30-11:45 ST-12 ST-99 ST-69 ST-87 17:20-17:40 Closing Ceremony

11:45-12:00 ST-98 ST-97

11.6 Departure

4.2 Program

 11.2 Plenary-1 8:30-12:00 WICC G204

Time Speaker No. Title

Chairman: Zhigang Zheng 8:30-9:00 Opening Ceremony 8:50-9:00 Photo Taking Chairman: Peter Hanggi 9:00-9:40 Yingcheng Lai PT-11 Two Surprising Phenomena in Relativistic Quantum Chaos 9:40-10:20 Kunihiko Kaneko PT-09 Dimension Reduction and Relevance of Slow Modes in Biological Dynamical Systems 10:20-10:40 Coffee break Chairman: Yingcheng Lai 10:40-11:20 Ravindra Amritkar PT-01 Extreme events on networks and spatially extended regions 11:20-12:00 Sergej Flach PT-03 Dynamical Glass

DDAP10, Huaqiao University, Xiamen China, 2018

 11.2 Session-A 13:30-18:00 WICC G201

Time Speaker No. Title Chairman: Byungnam Kahng 13:30-14:00 Tomoki ST-05 Robust and flexible sequential activities in a multiple Kurikawa timescale neural network 14:00-14:30 Jinzhi Lei ST-06 Dynamics of cancer development 14:30-15:00 Wei Lin ST-07 Dynamical time series analytics: From networks construction to dynamics prediction 15:00-15:30 Hai-Peng Ren ST-09 Nonlinear Dynamics of Soft Shaft Lifting- Rotating system in Silicon Crystal Puller 15:30-16:00 Coffee break Chairman: Tomoki Kurikawa 16:00-16:30 Gang Hu ST-03 Inference of targeted interactions with time delay based on fast-varying noise effects when most of nodes in networks hidden 16:30-17:00 Byungnam ST-14 Metastable group synchronization dynamics in a Kahng competing Kuramoto model and its related phenomena and models 17:00-17:20 Yuki Izumida ST-02 Nonlinear dynamics analysis of a low-temperature-differential kinematic Stirling engine 17:20-17:40 Chunli Huang ST-04 Competition of spiral waves in strong heterogeneous CGLE systems 17:40-18:00 Zhenzhen Lu ST-08 A new computational method for solving fractional diffusion and wave equation

DDAP10, Huaqiao University, Xiamen China, 2018

 11.2 Session-B 13:30-18:00 WICC G204

Time Speaker No. Title Chairman: Jiping Huang 13:30-14:00 Wei Wang ST-37 Fluctuation behaviors of the native states of globular proteins 14:00-14:30 Shiwei Yan ST-38 Statistical Physics and Quantitative Descriptions of Biological Systems 14:30-15:00 Zonghua Liu ST-42 Heat conduction in complex networks 15:00-15:30 Dahai He ST-40 Self-consistent phonon theory for thermal transport through anharmonic systems 15:30-16:00 Coffee break Chairman: Wei Wang 16:00-16:30 Jiping Huang ST-41 Novel controls of heat transfer via metamaterials 16:30-17:00 Jie Ren ST-43 Simulating Topological Phase Transition and Quantum Spin/Valley Hall Effects in Linear Circuit Networks 17:00-17:20 Wang Jian ST-44 Anomalous Energy Diffusion in two-dimensional nonlinear lattices 17:20-17:40 Daxing Xiong ST-45 Using Hilbert transform and classical chains to simulate quantum walks 17:40-18:00 Yong Zhang ST-46 One-dimensional superdiffusive heat propagation induced by optical phonon-phonon interactions

DDAP10, Huaqiao University, Xiamen China, 2018

 11.2 Session-C 13:30-18:00 WICC G101

Time Speaker No. Title Chairman: Ning Xu 13:30-14:00 Kazumasa A. ST-51 Universal scaling laws of growth processes and Takeuchi their implications for chaos instability 14:00-14:30 Zi-Ke Zhang ST-67 Dynamics of Epidemic Spreading and Information Diffusion on Social Networks 14:30-15:00 Wei Li ST-49 Reinforcement Learning for Complementarity Game and Population Dynamics 15:00-15:30 Zhi-Gang Shao ST-50 Contrasting complexity of the recent 2000 and past 122,000 years climates 15:30-16:00 Coffee break Chairman: Wei Li 16:00-16:30 Yilong Han ST-47 Compressing binary crystals into glasses 16:30-17:00 Ning Xu ST-53 Quasicrystal formation in monodisperse soft-core particle systems

17:00-17:20 Hao Liao ST-59 Temporal similarity metrics for network reconstruction based on diffusion dynamics 17:20-17:40 Jilin Jou ST-48 Dynamics of granular rafts on Faraday waves 17:40-18:00 Fenghua Wang ST-52 A robust collective credit allocation in scientific publications

DDAP10, Huaqiao University, Xiamen China, 2018

 11.2 Session-D 13:30-18:00 WICC G103

Time Speaker No. Title Chairman: Tian Qiu 13:30-14:00 Hyeon ST-70 Energetic costs, precision, and transport efficiency of Changbong molecular motors 14:00-14:3 Kyoung J. Lee ST-71 Unusual morphogenesis and cell motility in an in vitro 0 population of cancer cell lines 14:30-15:00 Chunhe Li ST-72 A landscape view on the interplay between EMT and cancer metastasis 15:00-15:30 Weixing Zhou ST-92 Multifractal analysis of financial markets: A review 15:30-16:00 Coffee break Chairman: Weixing Zhou 16:00-16:30 Tian Qiu ST-89 European stock market dynamics based on correlation analysis 16:30-17:00 Huailong Shi ST-91 Momentum and contrarian effects in the Chinese stock market: From perspective of the Adaptive Markets Hypothesis 17:00-17:20 Dehua Shen ST-93 Trading Volume and Return Volatility of Bitcoin Market: Evidence for the Sequential Information Arrival Hypothesis 17:20-17:40 Fei Ren ST-90 Evolution of volatility after extreme events 17:40-18:00 Jie Cao ST-88 The impacts of cross-shareholding network on extreme price movements: Evidence from China

DDAP10, Huaqiao University, Xiamen China, 2018

 11.3 Session-A 8:30-12:00 WICC G201

Time Speaker No. Title Chairman: Chusun Tian 8:30-9:00 Yueheng Lan ST-16 Kuramoto oscillator model in the view of renormalization group analysis 9:00-9:30 Xingang Wang ST-19 Measure synchronization in coupled quantum systems 9:30-10:00 Meng Zhan ST-13 New-generation power system dynamics 10:00-10:20 Coffee break Chairman: Yueheng Lan 10:20-10:45 Chusun Tian ST-21 Chaos-induced spin topological structure in kicked rotor 10:45-11:10 Weiqing Liu ST-17 Amplitude death induced by the periodical coupling 11:10-11:30 Wen Sun ST-10 The event-trigger control method for coupled oscillators with mismatch parameters 11:30-11:45 Wenjiao Sun ST-18 Global asymptotic synchronization and Mittag-Leffler synchronization of fractional-order memristive neural networks with time-varying delays and reaction-diffusion terms 11:45-12:00 Zhe Yin ST-11 An aged-structured SEIR epidemic model with discrete time delay and relapse

 11.3 Session-B 8:30-12:00 WICC G203

Time Speaker No. Title Chairman: Hao Ge 8:30-9:00 Bo Zheng ST-39 Temporal correlation functions of dynamic systems in non-stationary states 9:00-9:30 Penger Tong ST-34 Dynamic heterogeneity in nonequilibrium systems: from turbulence to living cells 9:30-10:00 Jiao Wang ST-35 Equilibrium Equality for Free Energy Differences 10:00-10:20 Coffee break Chairman: Bo Zheng 10:20-10:45 Hao Ge ST-27 Fluctuating-rate model of single-cell dynamics and its applications 10:45-11:10 Jae Dong Noh ST-94 Universal statistics of housekeeping entropy 11:10-11:30 An-Liang Cheng ST-26 Phase Transitions in Optimized Networks: Analytic Results and Efficient algorithm 11:30-11:45 Shanhe Su ST-33 The heat and work of quantum thermodynamic processes with quantum coherence 11:45-12:00 Jun Wang ST-36 Forces on nanoparticles in the free molecule regime

DDAP10, Huaqiao University, Xiamen China, 2018

 11.3 Session-C 8:30-12:00 WICC G101

Time Speaker No. Title Chairman: Takashi Nishikawa 8:30-9:00 Chin-Kun Hu ST-57 Complex networks and universal patterns in historical Chinese dynasties 9:00-9:30 Bing-Hong Wang ST-63 Relation between network structure and controllability transition in complex networks 9:30-10:00 Huijie Yang ST-66 Evolutionary Evidences in Structure of Complex Network 10:00-10:20 Coffee break Chairman: Huijie Yang 10:20-10:45 Takashi Nishikawa ST-60 Vulnerability and Co-susceptibility Determine Large Network Cascades 10:45-11:10 Chulan Kwon ST-29 Various heats in strongly coupled system and bath 11:10-11:30 Po Ho Fai ST-55 The impact of selfish route decisions on optimized transportation networks 11:30-11:45 Tai Tak Shing ST-61 The Benefits of Probabilistic Route Choices in Dynamical Transportation Networks 11:45-12:00 Chenbo Fu ST-56 NES-TL: Network Embedding Similarity-Based Transfer Learning  11.3 Session-D 8:30-12:00 WICC G103

Time Speaker No. Title Chairman: Hong Zhang 8:30-9:00 Ming-Chya Wu ST-73 Time series analysis approach to protein dynamics 9:00-9:30 Pik-Yin Lai ST-81 Cardiac Dynamics in response of periodic external stretching 9:30-10:00 Xiaodong Huang ST-80 Dynamics of the early afterdepolarization in the ventricular myocyte 10:00-10:20 Coffee break Chairman: Ming-Chya Wu 10:20-10:45 Hong Zhang ST-82 Jacobian-Determinant Method of Identifying Phase Singularity During Reentry 10:45-11:10 Zigang Huang ST-84 Network modeling for neuromodulatory intervention 11:10-11:30 Xiang Gao ST-79 A reaction-diffusion model to study the electrophysiological integration of human pluripotent stem cell-derived cardiomyocytes into the myocardial infarction zone 11:30-11:45 Xuhui Huang ST-83 Pattern recognition though spike-timing-dependent plasticity in spiking neural networks 11:45-12:00 Hengtong Wang ST-86 Resonance firing in single autaptic neuron

DDAP10, Huaqiao University, Xiamen China, 2018

 11.3 Plenary-2 13:30-18:00 WICC G204

Time Speaker No. Title Chairman: Zhisong Wang 13:30-14:10 Jurgen Kurths PT-10 Network-induced multistability: Lossy coupling and exotic solitary states 14:10-14:50 Adilson E. Motter PT-13 Symmetric Network States Requiring System Asymmetry 14:50-15:30 Masaki Sano PT-18 Nonequilibrium Phase Transition in Shear Flow 15:30-16:00 Coffee break Chairman: Chin-Kun Hu 16:00-16:40 Zhisong Wang PT-19 Re-think the link: 2nd law of thermodynamics and molecular motors 16:40-17:20 Lin I PT-07 Multi-scale collective micro-excitations and structural rearrangements in cold dusty plasma liquids 17:20-18:00 Eberhard PT-02 Synthetic Biology: A dream come true for the Bodenschatz quantitative study of dynamics and self-organization of life like processes

DDAP10, Huaqiao University, Xiamen China, 2018

 11.4 Plenary-3 8:30-12:00 WICC G204

Time Speaker No. Title Chairman: Changsong Zhou 8:30-9:10 Peter Hanggi PT-05 The ring of Brownian motion: Past ‐ Presence and Future Trends 9:10-9:50 Qi Ouyang PT-14 The free energy cost of Biological Circadian Clock 9:50-10:10 Coffee break Chairman: Qi Ouyang 10:10-10:50 Changsong Zhou PT-20 Cortical dynamics is cost-efficient in neural information representation 10:50-11:30 Peter Robinson PT-16 Nonlinear Dynamics in Systems Neuroscience via Neural Field Theory 11:30-12:10 Hawoong Jeong PT-08 Almighty Google knows almost everything! - Big-data and Network Science

There is no academic activity in the afternoon of 11.4.

DDAP10, Huaqiao University, Xiamen China, 2018

 11.5 Session-A 8:30-12:00 WICC G201

Time Speaker No. Title Chairman: Liang Huang 8:30-9:00 Yoshiyuki Y. ST-22 Collective 1/f fluctuation by pseudo-Casimir-invariants Yamaguchi 9:00-9:30 Tatsuo ST-23 Slow relaxation to equipartition and emergence of Yanagita quasi-equilibrium in slow-fast Hamiltonian systems 9:30-10:00 Ping Ao ST-25 New General Dynamical Structure for Nonlinear Processes: discovery of "fundamental theorem of dynamics" 10:00-10:20 Coffee break Chairman: Ping Ao 10:20-10:45 Liang Huang ST-20 Unconventional Scars in Chaotic Massive Dirac Billiards 10:45-11:10 Peng Ji ST-15 Regularization of synchronization dynamics by inertia, time delay and frustration 11:10-11:30 Li Chen ST-01 Spatial dynamics of interacting contagions 11:30-11:45 Conghui Xu ST-12 Dynamics of Stochastic Predator-Prey System with Age Structure and Crowley-Martin Functional Response

 11.5 Session-B 8:30-12:00 WICC G204

Time Speaker No. Title Chairman: Takahiro Hatano 8:30-9:00 Fabio ST-31 Multipurpose Broadband Cryogenic Mechanical Sensors Marchesoni 9:00-9:30 Sumiyoshi ST-24 Fractional kinetic equations and canonical formalism with Abe temporal nonlocality 9:30-10:00 Alberto ST-32 Self-organization conducted by the dynamics towards the Robledo attractor at the onset of chaos 10:00-10:20 Coffee break Chairman: Sumiyoshi Abe 10:20-10:45 Takahiro ST-28 Fracture of random media: mechanical instability and Hatano enhanced susceptibility 10:45-11:10 Ying Hu ST-95 Nonequilibrium Topological Insulator with Ultracold Atoms 11:10-11:30 Jae Sung Lee ST-30 Stochastic thermodynamics and fluctuation theorems with multiple reservoirs 11:30-11:45 Chen Wang ST-99 Quantum thermal transport in nonequilibrium spin-boson systems

DDAP10, Huaqiao University, Xiamen China, 2018

 11.5 Session-C 8:30-12:00 WICC G101

Time Speaker No. Title Chairman: Jie Sun 8:30-9:00 Deok-Sun Lee ST-58 Understanding the temporal pattern of spreading in heterogeneous networks: Time as a variable 9:00-9:30 Chenping Zhu ST-68 Principle of corresponding states of quasi-particle gases for American passenger flights 9:30-10:00 Tao Zhou ST-54 Predicting discredited behaviors of millions of enterprises 10:00-10:20 Coffee break Chairman: Tao Zhou 10:20-10:45 Gang Yan ST-65 Identifying the scale-free property in real networks 10:45-11:10 Jie Sun ST-96 Sensitive Dependence of Optimal Network Dynamics on Network Structure 11:10-11:30 Sheng-Jun ST-64 Self-organized criticality of earthquake model on Wang networks with randomness 11:30-11:45 Yueying Zhu ST-69 The formation of continuous opinion dynamics based on a gambling mechanism 11:45-12:00 Lijun Pei ST-98 Dynamics and Periodic solutions of two state dependent delay network congestion control models

 11.5 Session-D 8:30-12:00 WICC G103

Time Speaker No. Title Chairman: Xia-qing Shi 8:30-9:00 Leihan Tang ST-77 Front instability in microbial range expansion with chemotaxis 9:00-9:30 Xinliang Xu ST-78 Collective behaviors in active matter systems 9:30-10:00 Yun-Yun Li ST-75 Dynamics behavior of the Microswimmer 10:00-10:20 Coffee break Chairman: Leihan Tang 10:20-10:45 Xia-qing Shi ST-76 Dynamical Sub-classes of Dry Active Nematics 10:45-11:10 Xiongfei Fu ST-74 Spatial structure of bacterial behavior in chemotactic group migration 11:10-11:30 Yuanyuan Mi ST-85 Synaptic Correlates of Working Memory Capacity 11:30-11:45 Dong-Ping Yang ST-87 Unified Analysis of Global and Focal Aspects of Absence Epilepsy via Neural Field Theory of Corticothalamic System 11:45-12:00 Fei Qi ST-97 Bi-motility model to Bacterial Stripe Formation

DDAP10, Huaqiao University, Xiamen China, 2018

 11.5 Plenary-2 13:30-17:40 WICC G204

Time Speaker No. Title Chairman: Miguel A.F. Sanjuan 13:30-14:10 Miguel A.F. PT-17 Basin Entropy and testing for Wada basins to Sanjuan analyze the unpredictability of some physical systems 14:10-14:50 Carlo R. Laing PT-12 Chimeras in two-dimensional domains: heterogeneity and the continuum limit 14:50-15:30 Hiroshi H. PT-06 Nonequilibrium thermodynamics based on Hasegawa information geometry: Optimization problems based on information geometry 15:30-15:50 Coffee break Chairman: Hiroshi H. Hasegawa 15:50-16:30 Zhilin Qu PT-15 Nonlinear dynamics in the heart 16:30-17:10 Celso Grebogi PT-04 CHAOS-BASED QUANTUM CONTROL Chairman: Jibin Li 17:10-17:20 Best Presentation Award Award Presenters: Celso Grebogi, Miguel A.F. Sanjuan, Hiroshi H. Hasegawa 17:20-17:40 Closing Ceremony Zhigang Zheng

DDAP10, Huaqiao University, Xiamen China, 2018

V. Speakers and Abstracts

5.1 List of Speakers

5.1.1 Plenary Speakers

PT-01. Ravindra Amritkar, Institute of Infrastructure, Technology, Research and Management PT-02. Eberhard Bodenschatz, Max Planck Institute for Dynamics and Self-Organization PT-03. Sergej Flach, IBS Center for Theoretical Physics of Complex Systems PT-04. Celso Grebogi , The University of Aberdeen PT-05. Peter Hanggi, University of Augsburg PT-06. Hiroshi H. Hasegawa, Ibaraki University PT-07. Lin I, National Central University PT-08. Hawoong Jeong, Korea Advanced Institute of Science and Technology PT-09. Kunihiko Kaneko, The University of Tokyo PT-10. Jurgen Kurths, Potsdam Institute for Climate Impact Research & Humboldt University PT-11. Yingcheng Lai, Arizona State University PT-12. Carlo R. Laing, Massey University PT-13. Adilson E. Motter, Northwestern University PT-14. Qi Ouyang, Peking University PT-15. Zhilin Qu, University of California at Los Angeles PT-16. Peter Robinson, University of Sydney PT-17. Miguel A.F. Sanjuán, Universidad Rey Juan Carlos PT-18. Masaki Sano, University of Tokyo PT-19. Zhisong Wang, National University of Singapore PT-20. Changsong Zhou, Hong Kong Baptist University

DDAP10, Huaqiao University, Xiamen China, 2018

5.1.2 Session Speakers

Session A: Nonlinear Dynamics ST-01. Li Chen, Shaanxi Normal University ST-02. Yuki Izumida, Nagoya University, Nagoya 464-8601, Japan ST-03. Gang Hu, Beijing Normal University, CHINA ST-04. Chunli Huang, Beijing Normal University ST-05. Tomoki Kurikawa, Kansai Medical University ST-06. Jinzhi Lei, Tsinghua University ST-07. Wei Lin, Fudan University ST-08. Zhenzhen Lu, Beijing Jiaotong University ST-09. Hai-Peng Ren，Xi’an University of Technology ST-10. Wen Sun, Yangtze University ST-11. Zhe Yin, Beijing Jiaotong University ST-12. Conghui Xu, Beijing Jiaotong University ST-13. Meng Zhan, Huazhong University of Science and Technology Special Topic: Synchronization ST-14. Byungnam Kahng, Seoul National University ST-15. Peng Ji, Fudan University ST-16. Yueheng Lan, Beijing University of Posts and Telecommunications ST-17. Weiqing Liu, Jiangxi University of Science and Technology ST-18. Wenjiao Sun, Beijing Jiaotong University ST-19. Xingang Wang, Shaanxi Normal University Special Topic: Hamiltonian Dynamics ST-20. Liang Huang, Lanzhou University ST-21. Chusun Tian, Institute of Theoretical Physics, CAS ST-22. Yoshiyuki Yamaguchi, Kyoto University ST-23. Tatsuo YANAGITA, Osaka Electro-Communication University

Session B: Statistical Mechanics

ST-24. Sumiyoshi Abe, Huaqiao University ST-25. Ping Ao, Shanghai University ST-26. An-Liang Cheng, National Central University ST-27. Hao Ge, Peking University ST-28. Takahiro Hatano, University of Tokyo ST-29. Chulan Kwon, Myongji University ST-30. Jaesung Lee, Korea Institute for Advanced Study ST-31. Fabio Marchesoni, Tongji University ST-32. A. Robledo Universidad Nacional Autonoma de Mexico ST-33. Shanhe Su, Xiamen University ST-34. Penger Tong, Hong Kong University of Science and Technology ST-35. Jiao Wang, Xiamen University 25

DDAP10, Huaqiao University, Xiamen China, 2018

ST-36. Jun Wang, Beijing University of Technology ST-37. Wei Wang, Nanjing University ST-38. Shiwei Yan, Beijing Normal University ST-39. Bo Zheng, Zhejiang University Special Topic: phononics and applications ST-40. Dahai He, Xiamen University ST-41. Jiping Huang, Fudan University ST-42. Zonghua Liu, East China Normal University ST-43. Jie Ren, Tongji University ST-44. Jian Wang, Yangzhou University ST-45. Daxing Xiong, Fuzhou University ST-46. Yong Zhang, Xiamen University

Session C: Complex Systems

ST-47. Yilong Han, The Hong Kong University of Science and Technology ST-48. Jilin Jou, National Central University ST-49. Wei Li, Central China Normal University ST-50. Zhigang Shao, South China Normal University ST-51. Kazumasa A. Takeuchi, The University of Tokyo ST-52. Fenghua Wang, Beijing Normal University ST-53. Ning Xu, University of University of Science and Technology of China ST-54. Tao Zhou, University of Electronic Science and Technology of China Special Topic: Complex networks ST-55. PO Ho Fai, The Education University of Hong Kong ST-56. Chenbo Fu, Zhejiang University of Technology ST-57. Chin-kun Hu, Academia Sinica ST-58. Deok-Sun Lee, Inha University, Incheon 22212, Korea ST-59. Hao Liao, Shenzhen University ST-60. Takashi Nishikawa, Northwestern University ST-61. Tai Tak Shing, The Education University of Hong Kong ST-62. Bing Wang, Shanghai University ST-63. Binghong Wang, University of Science and Technology of China ST-64. Shengjun Wang, Shaanxi Normal University ST-65. Gang Yan, Tongji University ST-66. Huijie Yang, University of Shanghai for Science and Technology ST-67. Zike Zhang, Hangzhou Normal University ST-68. Chenping Zhu, Nanjing University of Aeronautics and Astronautics ST-69. Yueying Zhu, Wuhan Textile University

Session D: Interdisciplinary Dynamics

ST-70. Changbong Hyeon, Korea Institute for Advanced Study ST-71. Kyoung J. Lee, Korea University

DDAP10, Huaqiao University, Xiamen China, 2018

ST-72. Chunhe Li, Fudan University ST-73. Mingchya Wu, National Central University Special Topic: Active matters ST-74. Xiongfei Fu, Shenzhen Institutes of Advanced Technology, CAS ST-75. Yun-Yun Li, Tongji University ST-76. Xiaqing Shi, Soochow University ST-77. Leihan Tang, Beijing Computational Science Research Center ST-78. Xin-Liang Xu, Beijing Computational Science Research Center Special Topic: Cardiac Dynamics ST-79. Xiang Gao, Shaanxi Normal University ST-80. Xiaodong Huang, South China University of Technology ST-81. Pik-Yin Lai, National Central University ST-82. Hong Zhang, Zhejiang University Special Topic: Neural Dynamics ST-83. Xuhui Huang, Institute of Automation, CAS ST-84. Zigang Huang, Xi’an Jiaotong University ST-85. Yuanyuan Mi, Academy of Military Medical Sciences ST-86. Hengtong Wang, Shaanxi Normal University ST-87. Dongping Yang, University of Sydney Special Topic: Econophysics ST-88. Jie Cao, Central South University ST-89. Tian Qiu, Nanchang Hangkong University ST-90. Fei Ren, East China University of Science and Technology ST-91. Huailong Shi, Nanjing University of Information Science & Technology ST-92. Weixing Zhou, East China University of Science and Technology ST-93. Dehua Shen, Tianjin University

ST-94. Jae Dong Noh, University of Seoul ST-95. Ying Hu, Shanxi University ST-96. Jie Sun, Clarkson University ST-97. Fei Qi, Hong Kong Baptist University ST-98. Lijun Pei, Zhengzhou University ST-99. Chen Wang, Zhejiang Normal University

5.1.3 Posters

PS-01. Jianmin Shen, The non-equilibrium universality class of forget-remember mechanism on the (1+1)-d lattice PS-02. Wenjun Zang，The characteristics of Loop occupation base on networks classification by machine learning PS-03. LE Hao, Dynamic Simulation and Structural Optimization of Pulse Separation Device in Dual Pulse Solid Rocket Motor PS-04. Yanmeng Xing, Related research on failure phenomenon in scientific research 27

DDAP10, Huaqiao University, Xiamen China, 2018

PS-05. Jiangxing Chen, Chemically propelled motors navigate chemical patterns PS-06. Wen-Jing Zheng, A New Method of Identifying Vortex Filaments During Ventricular Arrhythmia PS-07. Wei Guo, Ergodicity breaking and ageing of underdamped Brownian dynamics with quenched disorder PS-08. Chen Zhanjun, Nonlinear Dynamic Response of High-aspect-ratio Wings PS-09. Lifang Dong, Square Superlattice Pattern in Dielectric Barrier Discharge PS-10. Nan Yao，Smart chimera masks perturbation through self-adaptive drifting PS-11. Hengzhe Yan, Reciprocity in one-dimensional nonlinear system PS-12. Yang Long, Floquet topological acoustic resonators and acoustic Thouless pumping PS-13. Danmei Zhang, Topology and zero-index in 1D phononic metamaterials with negative density and negative modulus PS-14. Luqin Wang，Graph theory analysis of the quantum thermal switch PS-15. Zi Wang, Non-equilibrium Restrictions on Nano-Scale Heat Pump by Fluctuation Theorem PS-16. Rundong Shi, Detecting Directed Interactions of Networks by Random Variable Resetting PS-17. Ke Cheng, Entropy production and infima statistics in the granular Brownian rotor PS-18. Yisen Wang, Complex Dynamics in the Carbon Nanotubes：A Molecular Dynamics Study PS-19. Xiaohui Dong, Non-Gaussian noise-weakened stability in a foraging colony system with time delay PS-20. Zhigang Zhu，Dynamics of recurrence in Non-integral Nonlinear Schrodinger Equations PS-21. Yafeng Wang, Synchronous patterns in coupled non-identical chaotic oscillators PS-22. Dapeng Zhang An abnormal phenomenon in pinning synchronization PS-23. Liang Wang Pinning cluster synchronization in regular networks PS-24. Lingbo Li, The Ising Model on Signed Networks PS-25. Ge Zhang, A class of chaotic oscillators dependent on the initial selection PS-26. Xinlin Song, Pattern transition of firing in a Morris-Lecar neuron induced autapse PS-27. Yue Liu Anomalous interfacial temperature profile induced by phonon localization PS-28. Ling Zhou, Fractal Aggregates on Geometric Graphs PS-29. Fan Wang, Effects of inhibitory signal on criticality in excitatory-inhibitory networks PS-30. Yanjiang Guo, Heat current through a weak interface PS-31. Changwei Huang, Persistence paves the way for cooperation in evolutionary games PS-32. S.F. Deng and W. Li, The spreading dynamics of forget-remember mechanism PS-33. Huijuan Wang, Cunquan Qu, Chongze Jiao, and Wioletta Rusze, Self-avoiding pruning random walk on signed networks PS-34. Tie Xu, Navigation with and without neural representation PS-35. Taro P. Shimizu and Takeuchi A. Kazumasa, Measuring Lyapunov spectrum of large chaotic systems with global coupling by time series analysis PS-36. Gao Fei, Li Wen-qin, Li Xi-ling, Chaotic dynamics of the Willis aneurysm system with Atangana–Baleanu derivative with fractional order and its control PS-37. Tianai Lu, Chunbiao Li, Coexisting attractors in chaotic system of conditional symmetry

DDAP10, Huaqiao University, Xiamen China, 2018

5.2 Abstracts

5.2.1 Plenary Talks (PT)

PT-01

Extreme events on networks and spatially extended regions Ravindra E. Amritkar Physical Research Laboratory Ahmedabad, India [email protected] Extreme events are very common in nature and many times they are associated with some calamities such as typhoons, earthquakes, floods etc. In many cases, there is some underlying network on which the extreme events take place, e.g. traffic jams on transport network, power black-outs on electrical grids, floods on river network etc. We study extreme events on complex networks using the random walk model where the extreme events are defined as surpassing of the flux above a prescribed threshold. We find that the nodes with smaller number of links are more prone to extreme events than the ones with larger number of links. In the case of biased random walks where the walkers may choose to preferentially hop towards the hubs or small degree nodes, the probability for the occurrence of extreme events on any node in the network depends on its generalized strength, a measure of the ability of a node to attract walkers. The generalized strength is a function of the degree of the node and that of its nearest neighbors. The result reveals that the nodes with a larger value of generalized strength, on average, display lower probability for the occurrence of extreme events compared to the nodes with lower values of generalized strength. It is well known that the frequency of extreme events depends on the location, e.g. some regions are prone to frequent earthquakes. We use the model of a Brownian particle in a potential for analysing the extreme events in spatially extended regions. We study the probability of observing extreme events in a region and this probability shows oscillations and exponential decay. However, if we see the general average behavior then we find that in general, the larger potential has a larger probability of extreme events while smaller potential has a smaller probability of extreme events. In our model, the inherent fluctuations in the model are responsible for giving an extreme event and the extreme event is not obtained by any external driving force. Thus, the extreme event is an integral part of the system and it will always occur.

PT-02

Synthetic Biology: A dream come true for the quantitative study of dynamics and self-organization of life like processes Eberhard Bodenschatz Max Planck Institute for Dynamics and Self-Organization, Germany [email protected]

DDAP10, Huaqiao University, Xiamen China, 2018

Abstract In my talk I will summarize the current initiatives in the study of bottom up synthetic biology. I will show how these approaches promise the development of well controlled experimental systems that allow quantitative tests of theories. This in turn will allow theory to going from descriptive to quantitatively predictive. I then will give examples from my laboratory, where I will address our discovery of cilia based directed fluid transport in the brain and will tie it to our bottom up development of artificial cilia. I will end with a brief discussion of the transformative potential impact of this scientific revolution, which carries with it a huge impact for the future from the mundane to the very specialized.

PT-03

Dynamical Glass Sergej Flach, IBS Center for Theoretical Physics of Complex Systems, KOREA [email protected] [email protected] Abstract: Classical many body interacting systems are typically chaotic and their microcanonical dynamics ensures that time averages and phase space averages are identical. In proximity to an integrable limit the properties of the network of nonintegrable action space perturbations help decide whether i) ergodization time scales stay on the order of the Lyapunov times, or whether ii) the system fragments into regular and chaotic parts and enters a dynamical glass (DG) phase at a finite distance to the integrable limit. This DG phase is induced by coherent localized excitations - generalized discrete breathers - with anomalously large lifetimes. We use a set of observables to quantify the properties of the DG phase. A sectioning of a typical trajectory by equilibrium Poincare manifolds defines events, whose statistics signals the onset of the DG phase. Distributions of finite time averages allow to extract time scales of ergodization, and their dependence on control parameters as compared to the Lyapunov times. I will show short range network results for a classical version of Josephson junction chains, first studied more than 30 years ago, and recently readdressed due to the prediction of a Many Body Localisation phase in the full quantum version of the system. I will also show results for systems with long range networks.

PT-04

CHAOS-BASED QUANTUM CONTROL Celso Grebogi University of Aberdeen, Scotland, UK Xi’an University of Technology, China http://www.abdn.ac.uk/icsmb/people/profiles/grebogi

Quantum chaos is referred to the study of quantum manifestations of systems that are chaotic in the classical limit. Most previous research in the field of quantum chaos focused on the non-relativistic quantum regime. Recently the field of relativistic quantum chaos has emerged, due to the tremendous 30

DDAP10, Huaqiao University, Xiamen China, 2018 development of research on graphene and topological insulators. Phenomena such as relativistic quantum scarring, chaotic scattering, and tunnelling have been explored. The speaker will discuss a number of fundamental issues in relativistic quantum chaos, but from the perspective of quantum control or modulation: how classical chaos can be exploited to harness relativistic quantum behaviours in Dirac fermion and graphene systems? Transport through quantum dot and resonant tunnelling will be used as two prototypical examples to illustrate the principle that chaos-based quantum control can be advantageous and experimentally feasible.

[1] Relativistic Quantum Chaos – An emergent interdisciplinary field, Y.-C. Lai, H.-Y. Xu, L. Huang, and C. Grebogi, Chaos 28, 052101(1-22) (2018) [2] Relativistic Quantum Chaos, L. Huang, H.-Y. Xu, C. Grebogi, and Y.-C. Lai, Physics Reports (2018) https://doi.org/10.1016/j.physrep.2018.06.006 [3] Harnessing Quantum Transport by Transient Chaos, R. Yang, L Huang, Y.-C. Lai, C. Grebogi, and L. M. Pecora, Chaos 23, 013125(1-9) (2013) [4] Effect of Chaos on Relativistic Quantum Tunnelling, X. Ni, L. Huang, Y.-C. Lai, and L.M. Pecora, European Phys. Lett. 98, 50007 (2012) [5] Do Dirac Fermions Scar in Chaotic Billiards? X. Ni, L. Huang, Y.-C. Lai, and C. Grebogi, Phys. Rev. E 86, 016702 (2012)

PT-05

The ring of Brownian motion: Past ‐ Presence and Future Trends Peter Hänggi e‐mail: [email protected]‐augsburg.de

Abstract: Since the turn of the 20‐th century Brownian noise has continuously disclosed a rich variety of phenomena in and around physics. The understanding of this jittering motion of suspended microscopic particles has undoubtedly helped to reinforce and substantiate those pillars on which the basic modern physical theories are resting: Its formal description provided the key to great achievements in statistical mechanics, the foundations of quantum mechanics and also astrophysical phenomena, to name but a few. ‐‐ Brownian motion determines the rate‐limiting step in most transport phenomena via escape events that help to overcome obstructing bottlenecks, does trigger firing events of neurons and in ion channels or, more generally, induces oscillatory dynamics in excitable media. Recent progress of Brownian motion theory involves (i) the description of relativistic Brownian motion and its impact for relativistic thermodynamics, or (ii) its role for fluctuation theorems and symmetry relations that constitute the pivot of those recent developments for nonequilibrium thermodynamics beyond the linear response regime. Although noise commonly is hold as an enemy of order, it in fact also can be of constructive influence. The phenomena of Stochastic Resonance and Brownian motors present two such archetypes wherein random Brownian dynamics together with unbiased nonequilibrium forces beneficially cooperate in enhancing detection and/or in facilitating directed transmission of information. The applications range from innovative information processing devices in physics, chemistry, and in physical biology to new hardware for medical 31

DDAP10, Huaqiao University, Xiamen China, 2018 rehabilitation. Particularly, those additional non‐equilibrium disturbances enable the rectification of haphazard Brownian noise so that quantum and classical objects can be directed along on a priori designed routes (i.e. Brownian motors). Despite its thrilling manifold successes Brownian motion is, nonetheless, not the ''Theory of Everything'', as is revealed by some more doubtful applications. We conclude with an outlook of future prospects, trends and unsolved issues. Home page and access to talks, publications: http://www.physik.uni‐augsburg.de/theo1/hanggi/

PT-06

Nonequilibrium thermodynamics based on information geometry: Optimization problems based on information geometry

Hiroshi H. Hasegawa1,2 and Tomomi Nakamura1

1: Dept. Math. Sci., Ibaraki University, Mito, Ibaraki 310-8512, Japan 2: Center for Complex Quantum Systems, Univ. of Texas Austin, Texas 78712, USA.

Abstract The maximum work formulation of the second law of thermodynamics has been generalized to transitions between nonequilibrium states [1][2]. The generalization involves the relative entropy (Kullback--Leibler divergence) between nonequilibrium states and canonical states. The canonical distribution takes on a novel role in bridging information and energy, even in thermally isolated systems. The relative entropy scaled by the temperature of the canonical state quantifies the work available for extraction from the nonequilibrium state. This scaled relative entropy can be interpreted as an energy- dimensional divergence in information geometry [3]. Amari’s extended Pythagorean theorem on three energy-dimensional divergences gives a geometrical interpretation of the generalized maximum work formulation. In this paper, we give the information-geometric foundation of the generalized second law. Our new concept of “thermodynamic divergence” has an important role in nonequilibrium statistical mechanics. We introduced “thermodynamic divergence” as an energy-dimensional Kullback--Leibler (KL) divergence scaled by temperature. The change of “thermodynamic divergence” gives us the dissipative work. Amari’s extended Pythagorean theorem based on three “thermodynamic divergences” gives us the geometric interpretation of the maximum work formulation of the generalized second law. From the nonequilibrium initial state, the canonical state with the effective temperature is the closest “point” in the “thermodynamic divergence.” The KL divergence scaled by a positive temperature satisfies the mathematical definition of divergence. It plays a role of “thermodynamic divergence” in information geometry. One important recognition in this paper is that the geometry based on the scaled divergences is different from the original. Amari’s extended Pythagorean theorem based on the original KL divergences gives the orthogonality between the one-parameter (the inverse of temperature) canonical-line and the equal-energy surface. On the other hand, Amari’s extended Pythagorean theorem based on the scaled KL divergences by temperature gives the orthogonality between the one-parameter (temperature) canonical- line and the isentropic surface (See Fig.1). We also give some applications of our approach based on information geometry. We introduce another KL divergence scaled by temperature over the fuel price (the number of particles times the

DDAP10, Huaqiao University, Xiamen China, 2018

price par particle in the heat reservoir). This scaled divergence measures the cost performance [4]. Amari’s extended Pythagorean theorem based on this scaled KL divergence gives us new important knowledge in optimizing the fuel consumption.

Fig.1 Amari’s extended Pythagorean theorem based on three scaled KL divergences

[1] H. H. Hasegawa et.al., Phys. Lett. A, 374, 1001 (2010); ibid. 375, 88 (2010); M. Esposite et.al., Europhys. Lett. 95. 40004 (2011); J. Ishikawa et.al. 2014 Entropy 16 (2014) 3471. [2] H. H. Hasegawa, T. Nakamura and D. Driebe, Chaos 27 (2017) 104606; T. Nakamura, H. H. Hasegawa and D. Driebe, (2018) submitted. [3] Amari S and Nagaoka H Methods of Information Geometry (vol 191 of Translations of Mathematical Monographs) (Providence: American Mathematical Society / Oxford: Oxford University Press) (2000); Amari S, Information Geometry and Its Applications, (Japan: Springer) (2016). [4] Kazuma Takara, Ph. D Thesis, Ibaraki University (2012).

PT-07

Multi-scale collective micro-excitations and structural rearrangements in cold dusty plasma liquids Lin I, Chi Yang, Wen Wang, and Hao-Wei Hu Department of Physics, National Central University, Chungli Abstract Microscopically, unlike intuitive expectation, the structure and the motion of the cold liquid around freezing are not completely disordered. Under solid-like dense packing, the interplay of strong particle mutual coupling and moderate thermal agitation leads to structural and dynamical heterogeneities with collective excitations over a wide range of spatiotemporal scales. Spatially, various sized crystalline ordered domains (CODs) coexist with defect clusters around COD interfaces. Temporally, particles alternately exhibit small amplitude rattling in cages formed by surrounding particles, and stick-slip type avalanche cooperative hopping over caging barriers in the form of clusters with various sizes, after accumulating sufficient constructive perturbations. The latter consequently induces structural rearrangement, which can be further complicated by external stress. The uncertain perturbations from the stochastic or slow drives and the complicated many-body interaction through the heterogeneous structure make understanding the excitations and propagations of multi-scale collective motions and their correlation with local structural order challenging issues. 33

DDAP10, Huaqiao University, Xiamen China, 2018

The dusty plasma liquid can be formed by micro-meter sized dust particles negatively charged and suspended in a low pressure gaseous discharge. Its proper spatial and temporal scales allow the direct visualization of micro-structure and individual particle motion over a large area. Its low pressure gaseous background also causes much lower damping rate than that of the colloidal system composed of charged colloids suspended in the dense liquid background, which has been commonly used for studying micro-dynamics of the liquid. The dust plasma liquid is therefore a good platform to mimic and understand the generic microscopic structural and dynamical behaviors for the liquid at the kinetic level. In this talk, our recent studies on the multi-scale collective micro-excitations and structural rearrangements in the quasi-2D cold dusty plasma liquid around freezing are reviewed. Unlike the string-like hopping in the hotter liquid, the multi-scale cooperative hopping through excitations of small co-rotating CODs, which cause cracking/healing/rotation of large CODs, associated with generation/propagation /recombination/annihilation of defects, are found to be the key process for structural rearrangement. It is also the basic mechanism for evolving to the more ordered structure after quenching, and the fractal like COD cracking and healing under weak external stress. The conservation of Burgers vectors is obeyed in the defect dynamical process. Through decomposing particle micro-motion with a continuous multi-scale power spectrum into self-similar modes of acoustic waves with different scales, the thermal induced motions can be viewed as microscopic wave turbulence exhibiting interesting spatiotemporal waveform dynamics. The regions with poor structural orders are the regions have stronger wave amplitudes for slow modes and more vulnerable to structural rearrangement. Key words: dusty plasma liquid, collective excitations, dynamical and structural heterogeneities, waveform decomposition and dynamics.

PT-08

Almighty Google knows almost everything! - Big-data and Network Science Hawoong Jeong Department of Physics, Center for Complex Systems, KAIST, Daejeon, 34141, Korea Abstract: Big-data is defined as "high-volume, high-velocity, and/or high-variety information assets that require new forms of processing to enable enhanced decision making, insight discovery and process optimization." by Gartner, Inc. This field of research has huge potential for practical applications but it also promises new discovery in science. Especially, when big-data from different domains are combined and analyzed together, it usually gives new insight and valuable results, and in this respect, network science can shed a light on making these big-data more useful. In this presentation, I will briefly review what we can do by combining big-data, especially using Google and network science together to study various complex systems such as social network between people, prediction of science and technology trends etc.

DDAP10, Huaqiao University, Xiamen China, 2018

PT-09

Dimension Reduction and Relevance of Slow Modes in Biological Dynamical Systems Kunihiko Kaneko Universal Biology Institute, University of Tokyo, Komaba, Meguro, Tokyo, 153-8902: email [email protected] Abstract A macroscopic theory for cellular states with steady-growth is presented, based on consistency between cellular growth and molecular replication, as well as robustness of phenotypes against perturbations. Adaptive changes in high-dimensional phenotypes are shown to be restricted within a low-dimensional slow manifold, from which a macroscopic law for cellular-states is derived, as is confirmed by adaptation experiments of bacteria under stress. Next, the theory is extended to phenotypic evolution, leading to proportionality between phenotypic responses against genetic evolution and by environmental adaptation. Evolutionary relevance of slow modes in controlling high-dimensional phenotypes is discussed. References 1. Kaneko K., Life: An Introduction to Complex Systems Biology, Springer (2006) 2. K. Kaneko, C.Furusawa, T. Yomo, "Macroscopic phenomenology for cells in steady-growth state", Phys.Rev.X(2015) 011014 3. C. Furusawa, K. Kaneko "Global Relationships in Fluctuation and Response in Adaptive Evolution", J of Royal Society Interface 12(2015), 20150482. 4. C. Furusawa, K. Kaneko " Formation of Dominant Mode by Evolution in Biological Systems” Phys. Rev. E 97(2018)042410 5. K. Kaneko, C. Furusawa “Macroscopic Theory for Evolving Biological Systems Akin to Thermodynamics”, Annual Rev. Biophys. (2018) 47, 273-290

PT-10

Network-induced multistability: Lossy coupling and exotic solitary states Juergen Kurths Potsdam Institute for Climate Impact Research & Humboldt University, Berlin, Germany & HUST, Wuhan, [email protected] Abstract The stability of synchronised networked systems is a multi-faceted challenge for many natural and technological fields, from cardiac and neuronal tissue pacemakers to power grids. In the latter case, the ongoing transition to distributed renewable energy sources is leading to a proliferation of dynamical actors. The desynchronization of a few or even one of those would likely result in a substantial blackout. Thus the dynamical stability of the synchronous state has become a focus of power grid research in recent years. We show that the non-linear stability against large perturbations is dominated and threatened by the presence of solitary state in which individual actors desynchronise. Remarkably, when taking physical losses in the network into account, the back-reaction of the network induces new exotic solitary states in the individual actors, and the stability characteristics of the synchronous state are 35

DDAP10, Huaqiao University, Xiamen China, 2018 dramatically altered. These novel effects will have to be explicitly taken into account in the design of future power grids, and their existence poses a challenge for control. Reference: P. Menck, J. Heitzig, N. Marwan, and J. Kurths, Nature Physics 9, 89 (2013) P. Menck, J. Heitzig, J. Kurths, and H. Schellnhuber, Nature Communication 5, 3969 (2014) P. Schultz, J. Heitzig, and J. Kurths, New J. Physics 16, 125001 (2014) Y. Zou, T. Pereira, M. Small, Z. Liu, and J. Kurths, Phys. Rev. Lett. 112, 114102 (2014) S. Leng, W. Lin, and J. Kurths, Scient. Rep. 6, 21449 (2016) V. Klinshov, V. Nekorkin, and J. Kurths, New J. Physics 18, 013004 (2016) F. Hellmann, P. Schultz, C. Grabow, J. Heitzig, and J. Kurths, Scient. Rep. 6, 29654 (2016) Y. Xu, Y. Li, H. Zhang, X. Li, and J. Kurths, Scient. Rep. 6, 31505 (2016) J. Wu, Y. Xu, H. Wang, and J. Kurths, CHAOS 27, 063105 (2017) J. Nitzbon, P. Schultz, J. Heitzig, J. Kurths, and F. Hellmann, New J. Phys. 19, 033028 (2017)

PT-11

Two Surprising Phenomena in Relativistic Quantum Chaos Ying-Cheng Lai Arizona State University Shaanxi Normal University

Abstract In this talk, two recently discovered phenomena in Relativistic Quantum Chaos will be presented: (1) a relativistic quantum chimera in pseudospin-1/2 Dirac material systems and (2) quantum chaotic scattering that defies Q-spoiling and Klein tunneling in pseudospin-1 Dirac-Weyl systems. The first phenomenon is the simultaneous emergence of classically integrable and chaotic dynamics and their quantum manifestations. The distinct relativistic quantum fingerprints associated with different electron spin states are due to a physical mechanism analogous to the chiroptical effect in the presence of degeneracy breaking. The phenomenon mimics a chimera state in classical complex dynamical systems but here in a relativistic quantum setting - henceforth the term “relativistic or Dirac quantum chimera.” Potential applications include enhancement of spin polarization, unusual coexisting quasi-bound states associated with distinct spin configurations, and spin selective caustics. The second phenomenon is the surprising occurrence of a class of surface plasmon-like modes in pseudospin-1 chaotic cavities, which maintain high-Q values in spite of fully developed chaos and super-Klein tunneling. This type of boundary modes absolutely have no counterpart in nonrelativistic quantum or even in pseudospin-1/2 systems. Experimental observations of the two phenomena are possible through, e.g., optical realizations of ballistic Dirac or Dirac-Weyl material systems.

Key words: quantum chaos, Dirac and Dirac-Weyl materials, relativistic quantum scattering, chimera, spin polarization, Klein tunneling, Q-spoiling, surface plasmon modes.

DDAP10, Huaqiao University, Xiamen China, 2018

PT-12

“Chimeras in two-dimensional domains: heterogeneity and the continuum limit” Carlo Laing Institute of Natural and Mathematical Sciences, Massey University (Albany), Private Bag 102-904, North Shore Mail Centre, Auckland, New Zealand Abstract We consider two different two-dimensional networks of nonlocally coupled heterogeneous phase oscillators. These networks were previously studied with identical oscillators and a number of spatiotemporal patterns found, mostly as a result of direct numerical simulation. Here we take the continuum limit of an infinite number of oscillators and use the Ott/Antonsen ansatz to derive continuum level evolution equations for order parameter-like quantities. Most of the patterns previously found in these networks correspond to relative fixed points of these evolution equations, and we show here the results of extensive numerical investigations of these fixed points: their existence and stability, and the bifurcations involved in their loss of stability as parameters are varied. Our results answer a number of questions posed by previous authors who studied these networks, and provide a better understanding of these networks’ dynamics.

PT-13

Symmetric Network States Requiring System Asymmetry Adilson E. Motter Department of Physics & Astronomy and Northwestern Institute on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, USA Abstract It is generally assumed that individual entities are more likely to exhibit the same or similar behavior if they are equal to each other—imagine animals using the same gait, lasers pulsing together, birds singing the same notes, and agents reaching consensus. In this presentation, I will show that this assumption is in fact false in networks of interacting entities. This surprising observation is rooted in a new network phenomenon we term “asymmetry-induced symmetry” (AIS), in which a stable state of the system can be symmetric only when the system itself is not. Using spontaneous synchronization as a model process, I will discuss scenarios where the state in which all nodes exhibit identical dynamics (i.e., a state of maximum symmetry) can only be realized when the nodes themselves are not identical or are non-identically coupled. AIS can be seen as the converse of the well-studied phenomenon of symmetry breaking and has far-reaching implications for processes that involve converging to uniform states. In particular, it offers a mechanism for yet-to-be-explained convergent forms of pattern formation, in which an asymmetric structure develops into a symmetric one. AIS also has implications for consensus dynamics, where it gives rise to scenarios in which interacting agents only reach consensus when they are sufficiently different from each other. Early references: T. Nishikawa and A.E. Motter, Symmetric states requiring system asymmetry, Phys. Rev. Lett. 117, 114101 (2016); Y. Zhang, T. Nishikawa, and A.E. Motter, Asymmetry-induced synchronization in oscillator networks, Phys. Rev. E 95, 062215 (2017); Y. Zhang and A.E. Motter, Identical 37

DDAP10, Huaqiao University, Xiamen China, 2018 synchronization of nonidentical oscillators: When only birds of different feathers flock together, Nonlinearity 31, R1 (2018).

PT-14

The free energy cost of Biological Circadian Clock Qi Ouyang1,2 1The State Key Laboratory for Artificial Microstructures and Mesoscopic Physics, School of Physics, Peking University, Beijing, 100871, China 2Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, AAIC, Peking University, Beijing, 100871, China Abstract Oscillation is an important cellular process that regulates timing of different vital life cycles. However, under the noisy cellular environment, the oscillation can be highly inaccurate due to phase fluctuation and its sensitivity to external noise can be insufficient. It remains poorly understood how biochemical circuits maintain their phase coherence and phase sensitivity, and what is the free energy cost of increasing the accuracy of timing. Here, we study different types of biochemical circadian clocks. We find that the coherence and sensitivity of the oscillation is enhanced by increasing free energy dissipation. For all the systems we studied, the phase diffusion constant follows the same inverse dependence on the free energy dissipation per period, and the phase sensitivity is proportional to the energy dissipation. We further analyzed the free energy partition of Kai oscillatory system in cyanobacteria. Using two different reaction models, we show that for both models, compared with the free energy that make the system oscillator, a large among of free energy is consumed by the information exchange among individual oscillators so that the molecules can oscillated synchronously. This analysis may explain the qualitative discrepancy between theoretical calculations and experimental measurements. Key Words Circadian clock, Free energy dissipation, Non-equilibrium thermodynamics, Kai system

PT-15

Nonlinear dynamics in the heart Zhilin Qu Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA Abstract The human heart beats once every second or so, totaling about 2–3 billion times in a normal life span. Under diseased conditions, the heart rhythm may suddenly become much faster and irregular, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is governed by nonlinear dynamics, as widely demonstrated by clinical, experimental, and theoretical studies in the last century. In this talk, I will present examples of nonlinear dynamics/ dynamical phenomena occurring in the heart and their links to normal cardiac functions and arrhythmias. I’ll also show how nonlinear dynamics may be useful in guiding the development of 38

DDAP10, Huaqiao University, Xiamen China, 2018 therapeutic strategies. I’ll discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

PT-16

Nonlinear Dynamics in Systems Neuroscience via Neural Field Theory P. A. Robinson School of Physics and Center for Integrative Brain Function, University of Sydney, Australia Abstract Going beyond well-known nonlinear effects at the single-neuron level, nonlinear effects at the systems and whole-brain levels are manifest in epileptic seizures (Hopf bifurcations, limit cycles, saddle cycles, hysteresis, chaos), migraines and visual hallucinations (Turing and Hopf-Turing patterns), strong visual stimulation (harmonic and subharmonic generation, phase locking, entrainment), deep brain stimulation therapy for Parkinson’s disease (harmonic and subharmonic generation, resonance suppression), natural 10 Hz alpha and 14 Hz spindle oscillations (bistability, harmonic generation), normal brain activity (marginal stability, critical slowing), and sleep-wake dynamics (bistability, hysteresis, ghost attractors). Recent analyses of the above nonlinear effects are surveyed here, along with comparisons with experimental results. The analyses presented are based on neural field theory, which averages over microscale dynamics to address mesoscale and macroscale phenomena. The ability of neural field theory to successfully predict a range of experimental outcomes highlights the power of this approach, for which new open-access codes are now available. It also opens up an array of unsolved nonlinear brain-dynamics questions for the attention of the research community.

PT-17

Basin Entropy and testing for Wada basins to analyze the unpredictability of some physical systems Miguel A. F. Sanjuán Nonlinear Dynamics, Chaos and Complex Systems Group, Departamento de Física, Universidad Rey Juan Carlos, Tulipán s/n, 28933 Móstoles, Madrid, Spain Abstract In nonlinear dynamics, basins of attraction are defined as the set of points that, taken as initial conditions, lead the system to a specific attractor. This notion appears in a broad range of applications where multistability is present, which is a common situation in neuroscience, economy, astronomy, ecology, and other disciplines. Nonlinear systems often give rise to fractal boundaries in phase space, hindering predictability. When a single boundary separates three or more different basins of attraction, we call them Wada basins. Usually, Wada basins have been considered even more unpredictable than fractal basins. However, this particular unpredictability has not been fully unveiled until the introduction of the concept of basin entropy. The basin entropy provides a quantitative measure of how unpredictable a basin is. With the help of several paradigmatic dynamical systems, we illustrate how to identify the ingredients that hinder the prediction of the final state. The basin entropy together with two new tests of 39

DDAP10, Huaqiao University, Xiamen China, 2018 the Wada property have been applied to some physical systems such as experiments of chaotic scattering of cold atoms, models of shadows of binary black holes, and classical and relativistic chaotic scattering associated to the Hénon-Heiles Hamiltonian system in astrophysics. [1] A. Daza, A. Wagemakers, M.A.F. Sanjuán, and J. A. Yorke. Testing for Basins of Wada. Scientific Reports 5, 16579 (2015) [2] A. Daza, A. Wagemakers, M.A.F. Sanjuán. Wada Property in Systems with Delay. CNSNS 43, 220–226 (2017) [3] A. Daza, A. Wagemakers, B. Georgeot, D. Guéry-Odelin, and M.A.F. Sanjuán. Basin entropy: a new tool to analyze uncertainty in dynamical systems. Scientific Reports 6, 31416 (2016) [4] A. Daza, B. Georgeot, D. Guéry-Odelin, A. Wagemakers, and M.A.F. Sanjuán. Chaotic dynamics and fractal structures in experiments with cold atoms. Phys. Rev. A 95, 013629 (2017) [5] A. Daza, A. Wagemakers, B. Georgeot, D. Guéry-Odelin, and M.A.F. Sanjuán. Basin entropy, a measure of final state unpredictability and its application to the chaotic scattering of cold atoms. In Chaotic, Fractional, and Complex Dynamics: New Insights and Perspectives. Mark Edelman, Elbert Macau and Miguel A. F. Sanjuán, Editors. Springer, Cham, 2018 [6] Alvar Daza, Jake O. Shipley, Sam R. Dolan and Miguel A. F. Sanjuan. Wada structures in a binary black hole system. Phys. Rev. D (2018) https://arxiv.org/pdf/1807.10741.pdf [7] A. Daza, A. Wagemakers, M.A.F. Sanjuán. Ascertaining when a basin is Wada: the merging method. Scientific Reports 8, 9954 (2018)

PT-18

Nonequilibrium Phase Transition in Shear Flow Masaki Sano 1Department of Physics, University of Tokyo, 7-3-1 Hongo, Tokyo 113-0033, Japan [email protected] Abstract The phase transitions to absorbing states, which the systems may enter but cannot escape from, have attracted considerable interests in recent years. Directed percolation (DP) is the absorbing phase transition with a single absorbing state and constitutes a wide class of universalities in non-equilibrium phase transitions. Recently, the transitions from laminar state to turbulent state in classical fluids, such as pipe flow, channel flow, and Taylor-Couette flow, have been tested experimentally and numerically whether the transition belong to this DP universality class or not. In my presentation, I will explain the recent progress in understanding this problem. I will also demonstrate that (2+1) dimensional DP is robust even in systems having non uniform advection and anisotropy in different experimental systems. References [1] M. Sano and K. Tamai, A Universal Transition to Turbulence in Channel Flow, Nature Physics, 12, 249 (2016). [2] K. A. Takeuchi, M. Kuroda, H. Chaté, and M. Sano, Phys. Rev. Lett. 99, 234503 (2007), Phys. Rev. E 80, 051116 (2009). [3] H. Ohnuma, M. Kuroda, and M. Sano, in preparation (2018).

DDAP10, Huaqiao University, Xiamen China, 2018

PT-19

Re-think the link: 2nd law of thermodynamics and molecular motors Zhisong Wang Physics Department, National University of Singapore, 2 Science Drive 3, Singapore 117542 E-mail: [email protected] Abstract Molecular motors from biology and man-made nanotechnology operate normally under isothermal condition (i.e. at a uniform temperature) due to swift heat transfer across the motor’s small dimension that readily levels any temperature difference. It has been long believed that isothermal molecular motors are not subject to any Carnot-type universal and clear-cut limit that quantifies the best of macroscopic heat engines permitted by the 2nd law. This is not necessarily true as the 2nd law requires a finite energy cost for any sustainable directional motion (i.e., the defining function of a motor) in an isothermal environment. Otherwise, a load would be attached to harvest the environmental heat under a uniform temperature, which violates the 2nd law (not the 1st law due to available environmental heat at any temperature above absolute zero). This well-known finite-cost requirement reminds of an inequality, and is often used as a lower bound to judge whether or not a molecular system may become a motor at all. But this 2nd-law limit implicates an upper bound too. Note that the finite-cost requirement applies even without any load, suggesting that it is a 2nd law-decreed energy price for pure direction rather than any work output. Then the least energy price allowed by the 2nd law for a certain level of ‘pure direction’ (or inversely the highest direction by a certain amount of energy) sets a meaningful upper bound for the best possible molecular motors. Now the question is how to quantify this 2nd-law decreed least energy price versus ‘pure direction’. The answer relies on how ‘pure direction’ is quantified. After a series of studies using the theoretical framework of stochastic thermodynamics, we closed in on a new concept of directional fidelity (instead of widely used forward-to-backward stepping ratio) for quantifying the pure direction, and formulated the 2nd-law-decreed least energy price for a directional fidelity in a surprisingly neat equality involving only the environmental temperature plus Boltzmann constant. This energy-fidelity relation is universal as it involves no model-dependent parameters and holds for arbitrary isothermal objects from self-propelled motors to externally driven particles. The strictest experimental test for the relation is ideally from a ~ 100% efficient molecular motor as it has no room to waste energy more than the 2nd law-decreed least price. Luckily, such a rotational motor exists in biology (F1-ATPase) and a new experiment allowing a straightforward test was just published when our relation was derived. The experiment-theory confrontation, which involves no free parameters, yields almost perfect agreement. The relation is also verified for force-induced directional motion by experimental data from three independent groups using biomotors – for translational as well as rotational motion. In the talk, I shall present a personal but unified view of the theoretical developments, applications to motor protein biology, and parallel nanotechnological efforts to develop artificial molecular motors towards the 2nd law-permitted best. Overall, a much deeper connection is emerging between the 2nd law and molecular motors than used to believed. It essentially elevates the 2nd law from a long perceived bottomline guidance for molecular motors (by a well-known inequality requirement) to a central governing principle even for the best motors (by a Carnot-type universal equality). From this new perspective, isothermal molecular motors become a model system as good as Watt’s heat engines for the purpose of exploring the 2nd law – now at its molecular root important and timely for the thriving nanotechnology as well as wider molecular sciences.

DDAP10, Huaqiao University, Xiamen China, 2018

PT-20

Cortical dynamics is cost-efficient in neural information representation

Changsong Zhou 1 Department of Physics, Hong Kong Baptist University, Kowloon Tong, Hong Kong 2 Centre for Nonlinear Studies and Beijing-Hong Kong-Singapore Joint Centre for Nonlinear and Complex Systems (Hong Kong), Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong 3 Beijing Computational Science Research Center

Abstract The brain is highly energy consuming, therefore is under strong selective pressure to achieve cost-efficiency in both cortical connectivities and activities. However, cost-efficiency as a design principle for cortical activities has been rarely studied. Especially it is not clear how cost-efficiency is related to ubiquitously observed multi-scale properties: irregular firing, oscillations and neuronal avalanches. Here we demonstrate that these prominent properties can be simultaneously observed in a generic, biologically plausible neural circuit model that captures excitation-inhibition balance and realistic dynamics of synaptic conductance. Their co-emergence achieves minimal energy cost as well as maximal energy efficiency on information capacity, when neuronal firing are coordinated and shaped by moderate synchrony to reduce otherwise redundant spikes, and the dynamical clusterings are maintained in the form of neuronal avalanches. Such cost-efficient neural dynamics can be employed as a foundation for further efficient information processing under energy constraint. Acknowledgments: This work was partially supported by Hong Kong Baptist University (HKBU) Strategic Development Fund, NSFC-RGC Joint Research Scheme HKUST/NSFC/12-13/01, NSFC (Grant Nos. 11275027) and NSFC (Grant Nos. 11225526 and 11121403). Reference: Dongping Yang, Haijun Zhou and Changsong Zhou*, Co-emergence of Multi-scale Cortical Activities of Irregular firing, Oscillations and Avalanches Achieves Cost-efficient Information Capacity, PLoS Computational Biology 13, e1005384 (2017).

DDAP10, Huaqiao University, Xiamen China, 2018

5.2.2 Session Talks (ST)

ST-01

Spatial dynamics of interacting contagions Li Chen Shaanxi Normal University Abstract The spread of infectious diseases, rumors, fashions, innovations are complex dynamical processes, embedded both in network and spatial contexts. Here we investigate the spatial dynamics of a sort of realistic contagions, where two agents, say A and B, interact with each other and diffuse simultaneously in the geographic space. The contagion dynamics for each follows the traditional Susceptible-Infected-Susceptible kinetics, and the interaction introduces a potential change in the rate of secondary infection. We found that in the simplest case where the diffusion coefficients for different compartments are uniform, novel wave propagation modes such as receding wave and standing wave are possible. In a more realistic scenario, where diffusion capability depends on the state of the agents, e.g. healthy individuals generally move faster than individual, non-trivial patterns emerges, and pattern hysteresis is revealed. These new dynamics emphasize the important role played by the spatial dimension, which could be underestimated in previous studies."

ST-02

Nonlinear dynamics analysis of a low-temperature-differential kinematic Stirling engine Yuki Izumida Department of Complex Systems Science, Graduate School of Informatics, Nagoya University, Nagoya 464-8601, Japan Abstract Innovative technology that makes us possible to utilize thermal energy in our daily life helps in solving a global energy and environmental problem. The low-temperature-differential (LTD) Stirling engines, which can operate under a small temperature difference between warmth of our hand and an air temperature [1], may serve as one such remarkable technology for our sustainable energy future. It is thus an important task to understand the working mechanism of the engine based on fundamental physics. In particular, simple description without losing essential features of the engine is quite useful. To this end, ideas from nonlinear dynamics could play an important role. In this talk, I present a nonlinear dynamics analysis of an LTD kinematic Stirling engine model proposed in [2]. The kinematic engine is a gamma-type Stirling engine, which is often used as one of the design of LTD Stirling engines. The engine mainly consists of gas as a working substance, a large cylinder with a displacer piston, a small cylinder with a power piston, and a crank. The pistons are

DDAP10, Huaqiao University, Xiamen China, 2018 connected to the crank with a fixed phase difference. The reciprocating motion of the displacer piston controls thermal contact of the gas with heat reservoirs attached to the large cylinder, while that of the power piston is converted into a rotational motion via the crank. The dynamical equations are derived based on mechanical engineering, thermodynamics, and ideas from nonlinear dynamics. The equation of motion of the crank and the time evolution equation of the gas temperature constitute the basic equations. Under the assumption of a time scale separation, the gas temperature can be adiabatically eliminated and the dynamical equations with two degrees of freedom as a dynamical system are obtained. It turns out that the equations describe the engine as a driven nonlinear pendulum motive-powered by the temperature difference. The rotational motion and the stationary state of the engine are described as a stable limit cycle and a stable fixed point of the dynamical equations, respectively. Numerical simulations show that for sufficiently large temperature differences, these two stable states coexist as observed in actual LTD engines. It is also shown that the limit cycle disappears via a homoclinic bifurcation as the temperature difference regarded as a bifurcation parameter decreases, explaining an actual behavior of the engine. In principle, the present bifurcation scenario predicted by the theory could be verified in experiments. The present model is expected to serve as a prototypical model for further researches on LTD Stirling engines. [1] J. R. Senft, An Introduction to Low Temperature Differential Stirling Engines, 4th edition (Moriya Press, Wisconsin, 2000). [2] Y. Izumida, EPL 121, 50004 (2018).

ST-03

Inference of targeted interactions with time delay based on fast-varying noise effects when most of nodes in networks hidden Gang Hu Department of Physics, Beijing Normal University, Beijing 100875, CHINA Abstract Most complex social, biological and technological systems can be described by dynamic networks. Reconstructing network structures from measurable data is a fundamental problem in almost all interdisciplinary fields. Network nodes interact one to other and those interactions often have diversely distributed time delay. Therefore, accurate reconstruction of any targeted interaction to any node requires measurements of all its neighboring nodes together with information of time delay. When networks are large, these data are often not available and time delay facts are deeply hidden. Here, we show fast-varying noises can be of help to solve these difficult problems. By computing suitable correlations we can infer the intensity and time delay of any targeted interaction with data of two nodes (driving and driven nodes) only while all other nodes are hidden. The method is analytically derived and fully justified by numerical simulations.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-04

Competition of spiral waves in strong heterogeneous CGLE systems Chunli Huang, Xiaohua Cui∗, Zengru Di School of Systems Science, Beijing Normal University, Beijing, 100875, China Abstract The effects of spatial heterogeneity in a model of two-dimensional complex Ginzburg-Landau equation (CGLE) are studied. With a large degree of the heterogeneity in two different media, there are three different patterns classified: (a) The spiral waves all be suppressed by the travelling waves that generated from the interface; (b) The spiral wave was swept away in one medium and remnant in another medium; (c) The coexistence of multiple spiral waves in different media. It is found that the interface is a wave source that can generate travelling waves to compete with the original spiral waves in two different media. The competition results depend on the frequencies of the original spiral wave and the generated travelling waves. It is shown that we can suppress and annihilate the spiral waves by selecting appropriate parameters to generate strong heterogeneity.

ST-05

Robust and flexible sequential activities in a multiple timescale neural network Tomoki Kurikawa Kansai medical university, Shinmachi, 2-5-1, Hirakata, Osaka, Japan Abstract In neural system, sequences of neural activities are important for identifying external stimuli, decision making and motor control. Many experimental studies demonstrated that these sequences are robust and reproducible against noise and initial states. Theoretically, there are many proposals to generate such a sequence of activities in neural networks. Usually, these sequential activities shape only simple orbits: once a neural state converges an orbit, sequential patterns are uniquely determined independent of context. However, sequential activities observed in neural system are flexibly changed against context and external input. Thus, how such robust and flexible neural sequential activities are shaped is important question for understanding information processing in neural system. In the present study, we show that simple learning rule requiring only local information, i.e., pre- and post-synaptic neural activities, can generate such a robust and flexible sequence with combination of two recurrent neural networks whose timescales are different. By using the network model, we demonstrated a robust context dependent sequence is formed. To reveal the underlying mechanism, we analyzed phase space of the recurrent network with faster time scale and found that attractors of the faster network are formed and collapsed depending on the neural state of the network with slower timescale. Our study suggests that a neural network with multiple timescales, which is observed in actual neural system, facilitates to generate robust and flexible sequential activities.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-06

Dynamics of cancer development Jinzhi Lei Tsinghua University Abstract Cancer development is a multi-scale long-term dynamics. In this talk, I will introduce a recent consideration of considering the progress of cancer development as a dynamical system of stem cell regeneration. I will show how a simple dynamical model of stem cell regeneration can tell us the insightful pathways for cancer development. Moreover, while we consider the heterogeneous in stem cell regeneration, we established single-cell based computation models that can be used to explore the detail dynamics of cancer development and drug resistance.

ST-07

Dynamical time series analytics: From networks construction to dynamics prediction Wei Lin School of Mathematical Sciences and ISTBI, Fudan University, Shanghai 200433, China Abstract In this talk, I will introduce two model-free frameworks of dynamical time series analytics. One framework is to detect the causation interactions among a large group of dynamical variables, which probably recovers a network hidden in a real-world system we are concerned. The second framework is to make a forecast or future prediction of dynamics based only on short-term and high-dimensional time series, which is usually believed to be a challenging task. Both frameworks use the advantages of Taken's embedding techniques, which reveals that utilization of dynamical system theory is more likely to exploit useful information from time series not only from the models but also from the real-world systems.

ST-08

A new computational method for solving fractional diffusion and wave equation Zhenzhen Lu, Yongguang Yu* Department of Mathematics, Beijing Jiaotong University, Beijing, 100044, P.R.China Abstract In this paper, a new computational method about the Riemann-Liouville fractional derivatives is proposed when 1<α<2. This new computational method is based on the Legendre cardinal functions for two-dimensional nonlinear diffusion and wave equation. In the proposed method, firstly, the solution is expanded by Legendre cardinal functions with unknown coefficients. Then, other functions also are expanded by Legendre cardinal functions with unknown cofficients. Third, the equation are employed to achieve a nonlinear algebraic equation, generalized Sylvester equation. Finally, Matlab is used to 46

DDAP10, Huaqiao University, Xiamen China, 2018 solve this problem. Keywords: Legendre cardinal functions，Generalized Sylvester equation，two dimensional nonlinear diffusion and wave equation.

ST-09

Nonlinear Dynamics of Soft Shaft Lifting- Rotating system in Silicon Crystal Puller Hai-Peng Ren Xian University of Technology Abstract Mathematical model of the soft shaft lifting and rotation system in monosilicon crystal puller is established. The dynamics of this system is investigated. The stability of the system is analyzed, bifurcation and period oscillation are reported in the system. Three routes to chaos are found in bifurcation diagram.

ST-10

The event-trigger control method for coupled oscillators with mismatch parameters Junxia Guan1, Wen Sun1, ZhigangZheng2, Shihua Chen3 1 School of Information and Mathematics, Yangtze University 2 College of Information Science and Engineering, Huaqiao University 3 College of Mathematics and Statistics, Wuhan University Abstract In order to reduce consumption of energy and occupancy of computing power, the event-triggered control methods are widely used in the cooperation and control of coupled identical oscillators. However, oscillators are non-identical in practice. In this paper, we design an event-trigger control protocol for coupled oscillators with mismatch parameters. A triggering function is constructed by using Lyapunov stability. The theoretical shortest trigger interval is established which ensures Zeno behavior does not exist in the system. Simulation results are presented to demonstrate the effectiveness of the proposed event-triggered control methods. Key Words The event-triggered control, Mismatch parameters, Triggering function

DDAP10, Huaqiao University, Xiamen China, 2018

ST-11

An aged-structured SEIR epidemic model with discrete time delay and relapse Zhe Yin and Yongguang Yu Department of Mathematics, Beijing Jiaotong University, Beijing, 100044, P.R.China Abstract This paper is considered with a polycyclic population dynamics of an aged-structured SEIR epidemic model of susceptible, exposed, infected, recovered and relapsed individuals. The model describes about a nonlinear transport system under four initial-boundary value problems for the nonlinear transport equations. Moreover, the equations were given under non-local integral boundary conditions and discrete time delay. And then, it was observed that the explicit recurrent formulae for computing the traveling wave solutions. These model parameters satisfy the restrictions to guarantee the existence and uniqueness of continuous solution. Afterwards, the aged-structured SEIR model with relapse was reduced to the nonlinear autonomous system of delay ODE of by using some insignificant simplifications. Dimensionless indicators and conditions of endemic equilibriums for the infection-induced diseases of the model are studied. In order to guarantee the stability in the neighborhood of the equilibrium, we can obtain the restrictions for the time delay. And we analyze of local asymptotical stability of the SEIR epidemic model with discrete delay and relapse. In the end, the numerical simulation of dynamics of SEIR population have been carried out to illustrate of our series of conclusions. Keywords: SEIR epidemic model; Aged-structured model; Asymptotical stability; relapse

ST-12

Dynamics of Stochastic Predator-Prey System with Age Structure and Crowley-Martin Functional Response Conghui Xu Yongguang Yu* Abstract In this paper, we research a stochastic predator-prey system with Crowley-Martin functional response and age structure, some dynamic properties of the system are researched. First, we obtain the existence and uniqueness of global positive solutions for the system, Then, we study the stochastically ultimately bounded of predator populations and prey populations. Furthermore, we give two sufficient conditions for extinction of predator populations and prey populations, one of them is the sufficient condition for the stochastically extinction of predator populations; the other is the sufficient condition for the extinction of prey populations. in the end, the results of numerical simulation confirm our conclusions. Keywords ： Stochastically ultimately bounded, Stochastically extinction, Brownian motion, Crowley-Martin functional response

DDAP10, Huaqiao University, Xiamen China, 2018

ST-13

New-generation power system dynamics Meng Zhan State Key Laboratory of Advanced Electromagnetic Engineering and Technology, Hubei Electric Power Security and High Efficiency Key Laboratory, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China ([email protected]) Abstract Power system is an important infra-structure, which has been believed as one of the most complicated man-made complex systems. Currently the second great revolution of power system is emerging in history, characterized with traditional power electric devices replaced by power-electronics-based devices gradually and massively. Meanwhile the dominant problems in power system dynamics, which are mainly related to electromechanical behavior of synchronous generator machine, have become more complicated, involved with much fast electromagnetism dynamics induced by power-electronics-based devices. Very recently, new power system oscillation induced by increase of renewable penetration have happened frequently in the world, which have not been fully solved and understood yet. In this talk, facing with these oscillation accidents, we will review the basic existing theories for the convectional power system dynamics and the newly power-electronics-based power system dynamics, and also introduce some of our recent works. We hope that this talk will stimulate broad interest of complex-system researchers on such a realistic important complex system, power system. Key Words: power system dynamics, power-electronics-based power system, power system oscillation

ST-14

Metastable group synchronization dynamics in a competing Kuramoto model and its related phenomena and models Byungnam Kahng Department of Physics and Astronomy, Seoul National University, Seoul 08826, Korea Abstract Here we consider a mixture of positive and negative coupling Kuramoto oscillators. Each oscillator with negative (positive) coupling attracts (repulses) every other oscillator. The synchronized phase of the mixed coupling system is characterized by two groups of oscillators roughly separated by an angle p in the phase space. They can either be static or traveling. When the natural frequency distribution of the oscillators is given uniform, the synchronization phase transition is hybrid. Near the hybrid critical point, we find two-step jump transition dynamics from incoherent to traveling wave state through a long-lasting metastable state with large fluctuations. Our explanation of the metastability is that the dynamic flow remains within a limited region of phase space and circulates through a few active states bounded by saddle and stable fixed points. This complex flow generates a long-lasting critical behavior, a signature of a hybrid phase transition. We show that the long-lasting period can be controlled by varying the density of inhibitory/excitatory interactions. We discuss a potential application of this transition behavior to the recovery process of human consciousness. Reference:

DDAP10, Huaqiao University, Xiamen China, 2018

1. Jinha Park and B. Kahng, Phys. Rev. E. 97, 020203(R) (2018).

ST-15

Regularization of synchronization dynamics by inertia, time delay and frustration Peng Ji Fudan University Abstract The interplay between time delay and phase shifts plays an important role for the regulation of synchronization and gives rise to interesting collective dynamical phenomena. Here, we derive a stability criterion for the synchronous state in a system of identical oscillators in the presence of inertia, time delay and a phase shift. We investigate how the combination of these factors influences the synchronization rates of the system approaching steady states after small perturbations. In particular, we prove that the global synchronization rate is twice the smallest local synchronization rate. We also find that the fastest synchronization rate increases with inertia, but occurs at smaller values of critical delay. Ensembles of oscillators coupled in complex topologies, with time delays, are also investigated.

ST-16

Kuramoto oscillator model in the view of renormalization group analysis Yueheng Lan Beijing University of Post and Telecommunication Kuramoto nonlinear oscillator model is famous for its simplicity and richness in dynamical behavior. In the limit of infinitely many oscillators, the dynamics could be well described by one or more mean-field equations which would be rederived with a renormalization group (RG) method elegantly in comparison with previous derivations. Furthermore, when the number of oscillators is finite but large, the RG method is still applicable and gives a hint to the analytic explanation of the finite-size scaling behavior.

ST-17

Amplitude death induced by the periodical coupling Weiqing Liu Jiangxi University of Science and Technology Abstract Amplitude death behavior of networked oscillators with periodic coupling is investigated. It is observed in simulation that the area of amplitude death could be significantly influenced by tuning the amplitude of the periodical coupling. When the variation of the coupling is positive, the increment of amplitude tends to enlarge the amplitude death domain. However, when the variation of the coupling is 50

DDAP10, Huaqiao University, Xiamen China, 2018 between negative and positive, the increment of amplitude tends to shrink the oscillation death domain. Theoretical results are presented and verify the numerical results well.

ST-18

Global asymptotic synchronization and Mittag-Leffler synchronization of fractional-order memristive neural networks with time-varying delays and reaction-diffusion terms Wenjiao Sun and Yongguang Yu Department of mathematics Beijing Jiaotong University, Beijing, 100044, P.R.China Abstract In this paper, a new memristive neural network system with time-varying delays is proposed by introducing the reaction-diffusion terms and fractional derivative. The global asymptotic synchronization and Mittag-Leffler synchronization of the system is studied. Firstly, in order to incorporate the long term memory into a reaction-diffusion neural network model, a generalized concept of differentiating including fractional values is introduced into the paper. Next, with the help of the Filippov solution defining on the set-valued mapping with differential inclusion, the solution of the discontinuous system make sense. Then, based on the method of Lyapunov and the comparison theorem, two controllers are designed. In the meantime, the sufficient conditions are given to guarantee the global asymptotic synchronization and Mittag-Leffler synchronization of the new system, respectively. Finally, the numerical examples are provided to show the correctness of reduction and validity of the controllers. The results of this paper extend the existing integer-order delayed reaction-diffusion neural networks. Key Words : Fractional order; Neural network; Reaction-diffusion terms; Time-varying delay; Memristive; Mittag-Leffler synchronization; Global asymptotic synchronization

ST-19

Measure synchronization in coupled quantum systems Xingang Wang School of Physics and Information Technology, Shaanxi Normal University, Shaanxi 710119 Abstract By the approach of the measure synchronization, we study the collective behavior of two coupled bosonic Josephson junctions. It is found that as the coupling strength increases, the two systems will be changed from the non-synchronous to synchronous states gradually. Remarkably, it is found that the introduction of a moderate mismatch to the system parameters could enhance synchronization, instead of deteriorate it. The mechanism underlying the observed phenomena is analyzed by tracing the measurement of the system trajectories in the phase space, which shows that under moderate parameter mismatch the trajectories undergo a significant change. Our study shed lights on the collective behavior of complex quantum systems, and might have implications to applications such as quantum information and computing. 51

DDAP10, Huaqiao University, Xiamen China, 2018

ST-20

Unconventional Scars in Chaotic Massive Dirac Billiards Min-Yue Song1, Zi-Yuan Li1, Hong-Ya Xu2, Liang Huang1, and Ying-Cheng Lai2 1Lanzhou University, Lanzhou, Gansu 730000, China 2Arizona State University, Tempe, AZ 85287, USA Abstract Quantum scar is an intriguing phenomenon in quantum (wave) systems that unveils classically unstable periodic orbits that otherwise will not be revealed directly in the classical world due to their vanishing measure [1,2]. This peculiar phenomenon has been observed in massless Dirac fermion systems [3-5]. In particular, due to the infinite mass confinement and the spinor type of wavefunction, there is an extra phase at each reflection at the boundary, which, unlike the Schrodinger case, is different from its time-reversed counterpart. This is originated from the chiral-symmetry-breaking of the system, and thus the name “chiral scars”. The massive Dirac fermion bridges the massless Dirac system and the conventional quantum system, as when the mass goes to infinity, the equation and the boundary condition of the first component returns to that of the Schrodinger billiard. We have investigated extensively the massive Dirac billiards and identified a new phase factor due to the finite mass, which changes the quantization condition completely, leading to unconventional scars with complex behaviors. In this talk, we will show the details of the origin, mechanism, and impact of this phase to the scarring states. Key Words: Relativistic Quantum Chaos, Quantum Scar, Massive Dirac Billiards, Reference [1] S.W. Mcdonald, A.N. Kaufman, ``Spectrum and eigenfunctions for a Hamiltonian with stochastic trajectories,'' Phys. Rev. Lett. 42, 1189–1191 (1979). [2] E.J. Heller, ``Bound-state eigenfunctions of classically chaotic Hamiltonian systems - scars of periodic-orbits,'' Phys. Rev. Lett. 53, 1515–1518 (1984). [3] L. Huang, Y.-C. Lai, D. K. Ferry, S. M. Goodnick, and R. Akis, ``Relativistic quantum scars,'' Phys. Rev. Lett. 103, 054101 (2009). [4] H.-Y. Xu, L. Huang, Y.-C. Lai, and C. Grebogi, ``Chiral scars in chaotic Dirac fermion systems,'' Phys. Rev. Lett. 110, 064102 (2013). [5] L. Huang, H.-Y. Xu, C. Grebogi, and Y.-C. Lai, ``Relativistic quantum chaos,'' Phys. Rep. 753, 1-128 (2018).

ST-21

Chaos-induced spin topological structure in kicked rotor Chusun Tian Institute of Theoretical Physics, CAS Abstract The kicked rotor is a “standard model” in studies of nonlinear dynamics. The kicked rotor without spin degree of freedom nowadays has been well studied. In this talk, I will present our recent result for a kicked rotor with spin degree of freedom. We find a dynamical phenomenon mathematically equivalent to the integer quantum Hall effect occurs, where Planck’s quantum mimics the magnetic 52

DDAP10, Huaqiao University, Xiamen China, 2018 field. I will show that this phenomenon is of chaos origin.

ST-22

Collective 1/f fluctuation by pseudo-Casimir-invariants Yoshiyuki Y. Yamaguchi Department of Applied Mathematics and Physics, Kyoto University, Kyoto 606-8501, Japan Kunihiko Kaneko Department of Basic Science, University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan Abstract The 1/f fluctuation is ubiquitously observed in nature and has attracted general interest over the past few decades. Multiple efforts have been made to understand its origins, but a coherent explanation remains one of the fundamentally unsolved problems. For instance, the superposition of Lorentzian spectra is a classical explanation, but it requires a distribution of multiple timescales and we need to explain an origin of the multiple timescales. In this presentation we restrict ourselves to the finite-size fluctuation in Hamiltonian systems having long-range interaction and provide a universal mechanism for the 1/f fluctuation emerging in collective dynamics of many degrees of freedom [1]. Target systems are self-gravitating systems, plasmas, Euler fluids, and other many-body systems. Dynamics of such systems are described by the Vlasov equation in the large population limit and the Vlasov equation has an infinite number of invariants. When the population is large but finite, the invariants are no longer exact and they become pseudo-invariants. The pseudo-invariants are effective as pseudo-constraints of dynamics up to a certain timescale which is proportional to the population and they provide slow dynamics which gives the 1/f fluctuation. Note that the Vlasov description and this scenario based on the invariants are valid in long-range systems. The proposed scenario is examined by systematic numerical experiments in thermal equilibrium states of the so-called α-Hamiltonian mean-field (α-HMF) model, where the parameter α represents the range of interaction. The α-HMF model contains the mean-field interaction and the nearest neighbor interaction: We can vary the range of interaction from long range to short range. There are two phases of the ordered and disordered phases and we examine the proposed scenario separately in each phase. The disordered phase is accessible for any range of interaction, and therefore, we compute directly the α (range of interaction) dependence of fluctuation of the order parameter. As expected, the 1/f fluctuation tends to disappear around the boundary of long range and short range as the population gets large. In contrast, the ordered phase is not realized with short-range interaction, but it has been reported recently that the fluctuation of the order parameter shows the two-step relaxation due to existence of the pseudo-invariants even if the initial state is put in thermal equilibrium [2]. This two-step relaxation is used to identify the timescale where the pseudo-invariants are effective. The obtained timescale is in agreement with the timescale where the 1/f fluctuation is observed. Finally, by introducing another model, we demonstrate that the existence of phase transitions in the α-HMF model is not essential to have the 1/f fluctuation and that the scenario extracts the hidden 1/f fluctuation. It is worth commenting that the invariants come from Hamiltonian structure of systems, and accordingly, a similar non-Hamiltonian system does not show the 1/f fluctuation. Nevertheless, the essence of the proposed scenario is existence of pseudo-invariants and this idea may be extended in wide fields. References: [1] Y. Y. Yamaguchi and K. Kaneko, Collective 1/f fluctuation by pseudo-Casimir-invariants, Phys. Rev. 53

DDAP10, Huaqiao University, Xiamen China, 2018

E 98, 020201(R) (2018). [2] Y. Y. Yamaguchi, Strange scaling and relaxation of finite-size fluctuation in thermal equilibrium, Phys. Rev. E 94, 012133 (2016).

ST-23

Slow relaxation to equipartition and emergence of quasi-equilibrium in slow-fast Hamiltonian systems Tatsuo Yanagita Osaka Electro-Communication University, Neyagawa, Osaka, 572-8530, JAPAN Abstract We study relaxation to equipartition in a slow-fast Hamiltonian system. By using the so-called bead-spring model which is one of the most popular polymer models, we report on relaxation to equilibrium. When the system has two distinct time scales, i.e, the fast vibration of bonds between beads and a typical time scale of the collision between the beads and solvent, the relaxation time obeys Boltzmann-Jeans law. In the course of the relaxation to equipartition, we found the following phenomena: i) an emergence of quasi-equilibrium state due to the effect of "freezing" for the high-frequency degrees of freedom, ii) induction from the quasi-equilibrium to equilibrium.

ST-24

Fractional kinetic equations and canonical formalism with temporal nonlocality Sumiyoshi Abe1,2,3 1Physics Division, College of Information Science and Engineering, Huaqiao University, Xiamen 361021, China 2Department of Physical Engineering, Mie University, Mie 514-8507, Japan 3Institute of Physics, Kazan Federal University, Kazan 420008, Russia Abstract Fractional kinetic theory has been attracting a lot of attention in the context of continuous time random walks and related semi-Markovian processes, which are of much relevance to the physics of anomalous diffusion phenomenon. Master equations appearing there require the use of fractional calculus, Riemann-Liouville calculus in particular, and have therefore temporal nonlocality. It is generically a hard task to obtain solutions of the equations in their analytic forms. Although numerical analysis is useful in such a situation, approximate analytic solutions are always valuable for getting good physical insights into the properties of the systems. Motivated by such an observation, the variational principle has recently been developed for the fractional Fokker-Planck equation. Since the equation does not possess time-reversal invariance, the state space is extended through introduction of an auxiliary field that is not a probability distribution. Combination of fractional Fokker-Planck equation and the fractional equation for the auxiliary field has been found to exhibit time-reversal structure in a peculiar way. Then, the so-called “Lévy Ansatz” has been proposed, and it has been shown how to explicitly obtain an analytic solution of fractional Fokker-Planck equation.

DDAP10, Huaqiao University, Xiamen China, 2018

Based on the action integral introduced in the above-mentioned variational principle and the associated Lagrangian density, here the field-theoretical Hamiltonian formulation is further developed for the fractional Fokker-Planck equation. For this purpose, Dirac’s generalized canonical formalism is employed, and the equation is recast into the Liouvillian form. However, here appears a new nontrivial issue, which may originate from the temporal nonlocality of the fractional Fokker-Planck equation. More explicitly, there turn out to exist two different Hamiltonians. This nonuniqueness of the Hamiltonian is different from the ordinary one related to total-time-derivative terms and can be interpreted in view of the fact that the canonical formalism is local in time. Thus, an important question arises if these two different Hamiltonians can generate identical time evolution of the system. To answer this question, an approach is developed for the setup of the non-equal-time Dirac-bracket relations with the help of the equations of motion, that is, the fractional Fokker-Planck equation and the fractional equation for the auxiliary field. Then, it is proven that both of the Hamiltonians generate identical time evolution, as desired. References [1] S. Abe, “Variational principles for fractional kinetics and the Lévy Ansatz”, Phys. Rev. E 88, 022142 (2013). [2] S. Abe, “Hamiltonian formulation of fractional kinetics”, to appear in Eur. Phys. J. Special Topics (e-print arXiv:1711.09697).

ST-25

New General Dynamical Structure for Nonlinear Processes: discovery of "fundamental theorem of dynamics" Ping Ao Shanghai Center for Quantitative Life Sciences, Shanghai University Abstract An innovative theoretical framework for stochastic dynamics based on the decomposition of a stochastic differential equation (SDE) into a dissipative component, a detailed-balance-breaking component, and a dual-role potential landscape has been developed, which has fruitful applications in physics, engineering, chemistry, and biology. Such development process is called the discovery of “fundamental theorem of dynamics”. It introduces the A-type stochastic interpretation of the SDE beyond the traditional Ito or Stratonovich interpretation or even the \alpha-type interpretation for multidimensional systems. The potential landscape serves as a Hamiltonian-like function in nonequilibrium processes without detailed balance, which extends this important concept from equilibrium statistical physics to the nonequilibrium region. A question on the uniqueness of the SDE decomposition was recently raised. Our review of both the mathematical and physical aspects shows that uniqueness is guaranteed. The demonstration leads to a better understanding of the robustness of the novel framework. In addition, we discuss related issues including the limitations of an approach to obtaining the potential function from a steady-state distribution. References: 1. SDE Decomposition and A-Type Stochastic Interpretation in Nonequilibrium Processes. RS Yuan, Y Tang, P Ao. Frontiers of Physics 12: 120201 (2017) 2. Potential Function in a Continuous Dissipative Chaotic System: Decomposition Scheme and Role of Strange Attractor. YA Ma, QJ Tan, RS Yuan, B Yuan, P Ao. Intl J Bifurcation and Chaos 24: 55

DDAP10, Huaqiao University, Xiamen China, 2018

1450015 (2014) 3. Beyond Ito vs. Stratonovich. RS Yuan, P Ao. Journal of Statistical Mechanics 2012: P07010 (2012) 4. On the Existence of Potential Landscape in the Evolution of Complex Systems. P Ao, C Kwon, H Qian, Complexity 12: 19-27 (2007) 5. Structure of Stochastic Dynamics near Fixed Points. C Kwon, P Ao, DJ Thouless. Proc. Nat’l Acad. Sci. (USA) 102: 13029-13033 (2005) 6. Potential in Stochastic Differential Equations: Novel Construction. P Ao. J. Phys. A37: L25-L30 (2004)

ST-26

Phase Transitions in Optimized Networks: Analytic Results and Efficient algorithm An-Liang Cheng and Pik-Yin Lai Department of Physics and Center for complex system, National Central University, Taoyuan City 320 Abstract In this research, we focus on the principle of the formation or growth of the networks. We aim at constructing some basic theory for the optimal growing networks and investigate the drastic structural changes of the networks such as the phase transition. We consider the network formation governed by two aspects. The first is the material part, which represents the energy or cost require to construct the network. The second part is the function relevant for a certain properties that is to be maximized in designing the network. The goal is to minimize the cost of material but at the same time trying to get the maximal important network properties such as the network connections. We use the classical network models and statistical mechanics methods to study the growth of some networks under the principle of optimizing some objective function. Our motivation is to understand the relation between the energy cost and the connectivity in a network. By mapping the system to an Ising spin model, we obtain analytic results for two such models. The first model is maximizing total connections and the second model is collective edge formation, both of them show interesting, but different phase transition behaviors for general wiring cost distributions. The phase diagrams for these transitions are also obtained. These results are also extended for networks optimized with the weighted nodes degree for connections. Furthermore, mean-field theory leads to an effective algorithm for finding the fully optimized network in these models. We compare the minimal cost found by the new algorithm with the results found by solving the coupled mean-field equations. The results clearly verify that the new algorithm finds the minimal cost in all cases, but with a much faster computational time. All these results are also verified by Monte Carlo simulations. We investigate the degree distribution and the edge weight distribution of the network that are optimized by our model. Some interesting properties of the fully optimized network are also discussed in the light of small-world-ness, rich-club effects and the 80-20 rules such as the Pareto distribution. Key Words: Optimized networks, Phase transition, Ising model, Mean-field theory, Small-world-ness, rich-club effect, 80/20 rule

DDAP10, Huaqiao University, Xiamen China, 2018

ST-27

Fluctuating-rate model of single-cell dynamics and its applications Hao Ge Peking University Abstract Stochastic processes become more and more popular to model the mesoscopic biophysical dynamics, especially in single-cell biology. We proposed a fluctuating-rate model for the stochastic biochemical dynamics in a single cell, which is indeed piecewise deterministic Markov process. We also found that the fluctuating-rate model yields a nonequilibrium landscape function, which, similar to the energy function for equilibrium fluctuation, provides the leading orders of fluctuations around each phenotypic state, as well as the transition rates between the two phenotypic states. The rigorous proof needs to integrate the well-known Donsker-Varadhan theory and Feidlin-Wentzell theory in such an averaging case. We then apply this model to Lac operon system, and show how and why the stochastic gene-state switching can significantly broaden the environmental parameter ranges for the existence of bistability induced by positive feedback, which can be beneficial dealing with unpredictable environmental changes.

ST-28

Fracture of random media: mechanical instability and enhanced susceptibility Takahiro Hatano1 and Subhadeep Roy2 1. Earthquake Research Institute, the University of Tokyo, Japan 2. Norwegian University of Science and Technology, Norway Abstract Fracture of solids is undoubtedly an important issue in natural sciences, engineering, and industries. It thus has a long history of study, but most theoretical studies in the literature are based on the elasticity theory and involve homogeneous solids. Natural solids such as rocks or sands may compose of various ingredients, the mechanical properties of which have wide distributions. We refer to such solids as random media. In an atomistic view, amorphous materials such as glass are also such random media. In a much larger scale, the upper crust of the Earth can be regarded as a random medium, as it consists of various kinds of rocks. Thus, random media are ubiquitous in many branches of science. Here we discuss the combination of the above two problems: fracture of random media. Due to the randomness, conventional theories based on the elasticity theory have only a limited capability. For instance, application of the elasticity theory to granular solids has been confronted with fundamental difficulties. Instead, discrete modeling has been an effective tool for most aspects of fracture problems in random media: dynamics and statistical properties. In this talk, we introduce a statistical-mechanical model for fracture of random media: the fiber bundle model. This is an assembly of “fibers” that support the external load in parallel. Each fiber has its own failure threshold, which is randomly set according to a specific probability distribution function. The heterogeneity in the elasticity can be also incorporated. Note that the “fibers” here are just a 57

DDAP10, Huaqiao University, Xiamen China, 2018 metaphor, and they can be any ingredients of a system, as long as they support the external load in parallel. Using this model, we explain some power-law behaviors commonly observed in fracture. In particular, we show that mechanical instability of the model is due to the saddle-node bifurcation, and the power-law behaviors result from this bifurcation. We also discuss the relevance of this model to earthquake statistics.

ST-29

Various heats in strongly coupled system and bath Chulan Kwon1, Jaegon Um2, Hyunggyu Park3, Joonhyun Yeo4 1Department of Physics, Myongji University, Yongin, Gyeonggi-Do,17058, Korea 2 Department of Physics, POSTECH, Pohang, Gyeongbuk 37673, Korea 3School of Physics, Korea Institute for Advanced Study, Seoul 02455, Korea 4Department of Physics and Astronomy, Konkuk University, Seoul 08826, Korea Abstract We investigate three kinds of heats produced in system and bath strongly coupled via interaction

Hamiltonian. We define QS and QB as the energy loss of system and the energy gain of bath, respectively, mediated by the change in interaction Hamiltonian through dynamics. We consider bath with a thermostat which prevents it from heating up indefinitely so as to be able to reach either equilibrium or nonequilibrium steady state. Then, we define another kind of heat QSB as the energy dissipation of bath towards super bath surrounding it. The three quantities are equal only in average sense in steady state, but different in transient or time-dependent state. It is contradictory to our common sense about the equivalence of QS and QB. We find that there are many possible fluctuation theorems for entropy production to relate with the Kullback-Leibler divergence between guesssed and true states. We show that the distribution functions of thermal fluctuations for the three heats are different even in steady state.

ST-30

Stochastic thermodynamics and fluctuation theorems with multiple reservoirs Jae Sung Lee and Hyunggyu Park Korea Institute for Advanced Study [email protected] Abstract We reformulate stochastic thermodynamics in terms of noise realizations for Langevin systems in contact with multiple reservoirs and investigated the structure of the second laws of thermodynamics. We derive a hierarchy of ﬂuctuation theorems when one degree of freedom of the system is aﬀected by multiple reservoirs simultaneously, that is, when noise mixing occurs. These theorems and the associated second laws of thermodynamics put stricter bounds on the thermodynamics of Langevin systems. We apply our results to a stochastic machine in noise-mixing environments and demonstrate that our new bounds play a crucial role in determining the potential function and performance of the machine.†. † J. S. Lee and H. Park, New J. Phys. 20 083010 (2018). 58

DDAP10, Huaqiao University, Xiamen China, 2018

ST-31

Multipurpose Broadband Cryogenic Mechanical Sensors Fabio Marchesoni1,2 1 Department of Physics, University of Camerino, Camerino, 62032, Italy 2Center for Phononics and Thermal Energy Science, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China Abstract We study the response of a monolithic broadband folded pendulum operated inside a cryostat as horizontal displacement sensor (without force feedback). Measurements were taken in the temperature range 3-300K and compared with the internal friction in the device’s mechanical structure; we conclude that its sensitivity is mostly limited by the material losses. This suggests the possibility of adopting this class of devices as flexible tools for (i) investigating the mechanical properties of metals at low temperature and low frequencies, (ii) real=time monitoring of human artefacts (buildings and infrastructures) against seismic and other environmental risks; (iii) accurate steering of spacecraft.

ST-32

Self-organization conducted by the dynamics towards the attractor at the onset of chaos Alberto Robledo, Instituto de Física y Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico, [email protected] Abstract We construct an all-inclusive statistical-mechanical model [1] for self-organization based on the hierarchical properties of the nonlinear dynamics towards the attractors that define the period-doubling route to chaos [2-4]. The aforementioned dynamics imprints a sequential assemblage in the model that privileges progressively lower entropies, while a new set of configurations emerges due to the partitioning of collections of degrees of freedom of the original system into secluded portions. The initial canonical balance between numbers of configurations and Boltzmann-Gibbs (BG) statistical weights is drastically altered and ultimately eliminated by the sequential actions of the attractor. However the emerging set of configurations implies a different and novel entropy growth process that eventually upsets the original loss and has the capability of locking the system into a self-organized state with characteristics of criticality, therefore reminiscent in spirit to the so-called self-organized criticality. Some specifics of the approach we develop are: We systematically eliminate access to configurations of an otherwise conventional thermal system model by progressively partitioning it into isolated portions until only remains a subset of configurations of vanishing measure. The thermal system consists of a large number of effective degrees of freedom, each occupying entropy levels with the form of inverse powers of two. The sequential process replaces the original configurations by an emerging discrete scale invariant set of configurations of ensembles. Running this process has the following consequences: 1) The constrained thermal system becomes a close analogue of the dynamics 59

DDAP10, Huaqiao University, Xiamen China, 2018 towards the multifractal attractor at the period-doubling onset of chaos. 2) The new emerging numbers of configurations grow more slowly than exponentially, such that the developing statistical-mechanical properties depart from those of the ordinary Boltzmann-Gibbs form and acquire features from q-statistics. [1] Robledo A., “Self-organization and a constrained thermal system analogue of the onset of chaos”, Europhysics Letters, 123, 40003 (2018). [2] Robledo A., Moyano, L.G., “q-deformed statistical-mechanical property in the dynamics of trajectories an route to the Feigenbaum attractor”, Physical Review E 77, 032613 (2008). [3] Robledo, A., “A dynamical model for hierarchy and modular organization: The trajectories en route to the attractor at the transition to chaos”, Journal of Physics: Conf. Ser. 394, 012007 (2012). [4] Diaz-Ruelas, A., Robledo, A., “Emergent statistical-mechanical structure in the dynamics along the period-doubling route to chaos”, Europhysics Letters 105, 40004 (2014).

ST-33

The heat and work of quantum thermodynamic processes with quantum coherence Shanhe Su Xiamen University Abstract We study quantum thermodynamics for a time-dependent system in a non-adiabatic evolution. The definitions of heat and work are clarified for the general time evolution with quantum coherence characterized by the off-diagonal elements of the density matrix. By integrating the coherence effects in the exactly-solvable dynamics of quantum-spin precession, the internal energy is transferred as work rigorously in the thermodynamic adiabatic process. The present study demonstrates that quantum adiabatic process certainly results in a thermodynamic adiabatic process, but not all thermodynamic adiabatic processes are due to the quantum-mechanical adiabatic processes.

ST-34

Dynamic heterogeneity in nonequilibrium systems: from turbulence to living cells Penger Tong Department of Physics, Hong Kong University of Science and Technology, Hong Kong Correspondence e-mail address: [email protected]

The motion of particles at thermal equilibrium usually has a finite correlation time and will eventually be randomized after a long delay time, so that their displacement follows the Gaussian statistics, as dictated by the central limit theorem. Non-Gaussian fluctuations with an exponential tail in their probability density function (PDF), on the other hand, are often observed in nonequilibrium steady states (NESSs) and one does not understand why they appear so often. In this talk, I will present two examples of NESSs, in which we find the exponential PDF is directly linked to dynamic heterogeneity in the system. One example is turbulent Rayleigh-Bénard convection (RBC), in which warm fluid rises, cold fluid falls, and their mixing produces convective turbulence. The measured PDF P(δT) of 60

DDAP10, Huaqiao University, Xiamen China, 2018 temperature fluctuations δT in the central region of the flow was found to have a long exponential tail [1]. The other example is the lateral motion of acetylcholine receptors (AChRs) on live muscle cell membrane. The AChRs are linked to fluorescent quantum dots and their trajectories are tracked with nanometer accuracy over a broad range of times. The measured PDF of the protein displacements showed a consistent exponential tail, which is robust and universal for cells under different conditions [2]. In both cases, we found that because of the dynamic heterogeneity, the unconditional PDF P(x) can be written as a weighted convolution of a set of fluctuation modes G(x|) conditioned on a constant dynamic variable . For RBC, the correct dynamic variable  is the thermal dissipation rate. For protein diffusion, it is the molecular diffusivity. The conditional PDF G(x|) under a constant  is found to be of Gaussian form and its variance σ2() for different values of  follows an exponential distribution. The convolution of the two distribution functions gives rise to the exponential PDF P(x). Our work thus provides a physical mechanism of the observed exponential PDF in RBC and protein diffusion. It also sheds light on the origin of non-Gaussian fluctuations in other NESSs. The work on RBC was done in collaboration with Xiaozhou He and Yin Wang. The work on protein diffusion was done in collaboration with Wei He, Hao Song, Yun Su, Lin Geng, B. J. Ackerson, and Ben Peng. Both research projects were supported by the Research Grants Council of Hong Kong SAR. [1] X. He, Y. Wang, and P. Tong, Phys. Rev. Fluids 3, 052401(R) (2018). [2]W. He, H. Song, Y. Su, L. Geng, B. J. Ackerson, H. B. Peng, and P. Tong, Nature Communications, 7:11701 (2016).

ST-35

Equilibrium Equality for Free Energy Differences Jiao Wang Department of Physics, Xiamen University, Xiamen 361005, China [email protected] Abstract In 1997, Jarzynski derived an expression for the free energy difference (FED) between two equilibrium states of a system in terms of the work performed during a nonequilibrium process by which the system is driven from one equilibrium state to the other [1]. This expression is known as the Jarzynski equality (JE), which provides a novel nonequilibrium method for measuring and calculating the FED [2]. We find that [3], by adopting a particular nonequilibrium process involving a stage of integrable dynamics, partial of the work can be solved analytically and as a result, the nonequilibrium JE reduces to an equilibrium equality instead. As an application, the FED can be numerically computed based on this equality by sampling the canonical ensemble of one equilibrium state and the computational cost can be significantly saved compared with the method based on the JE. Other related issues will be discussed as well. References: [1] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997) [2] C. Dellago and G. Hummer, Entropy 16, 41-61 (2014) [3] L. Zhu and J. Wang, Phys. Rev. E 98, 022117 (2018)

DDAP10, Huaqiao University, Xiamen China, 2018

ST-36

Forces on nanoparticles in the free molecule regime Jun Wang Beijing University of Technology Abstract The transport properties and the force calculation for nanoparticles in dilute gases are of great interest to many areas, including the aerosol science, chemical engineering, and air pollution. The characteristic size of the nanoparticle is much less than the mean free path of the dilute gas. Thus, the transport properties of nanoparticles are much different from those of larger particles. Based on the gas kinetic theory and by considering the gas-particle interactions, we obtained the formulae for the drag force, thermophoretic force and lift force of nanoparticles in dilute gases. It is found that the gas-particle intermolecular interactions play important roles in evaluating these forces, and the thermophoretic force and lift force might reverse directions in a certain temperature range.

ST-37

Fluctuation behaviors of the native states of globular proteins Wei Wang National Lab of Solid State Microstructure and Department of Physics, Nanjing University, Nanjing 210093, China Abstract The three-dimensional folded structures of proteins, known as native states, make proteins capable of performing related biological functions. Based on the protein structural ensembles of the native states determined by NMR, we study the position fluctuations of residues and find high susceptibility and long-range correlations up to the protein sizes, showing a critical behavior with universal scaling relations. Our analyses are also extended to the studies on the protein conformations during folding.

ST-38

Statistical Physics and Quantitative Descriptions of Biological Systems Shiwei Yan Department of Physics, Beijing Normal University Abstract Systems biology aims to move beyond the study of single biomolecules and the interaction between specific pairs of molecules; its goal is to describe, in quantitative terms, the dynamic systems behavior of complex biological systems that involve the interaction of many components. Two related developments are currently changing traditional approaches to computational systems biology modelling. First, stochastic models are being used increasingly in preference to deterministic models to describe biochemical network dynamics at the single-cell level. Second, sophisticated

DDAP10, Huaqiao University, Xiamen China, 2018 statistical methods and algorithms are being used to fit both deterministic and stochastic models to time course and other experimental data. Both frameworks are needed to adequately describe observed noise, variability and heterogeneity of biological systems over a range of scales of biological organization. I will review how we have applied stochastic thermodynamics and statistical theory on the study of biochemical reactions. As an example, I will present a set of stochastic delay-differential equations which is used to study the temporal behaviors of tumors suppressor p53 protein.

ST-39

Temporal correlation functions of dynamic systems in non-stationary states Bo Zheng Department of Physics, Zhejiang University, Hangzhou 310027, P. R. China Email: [email protected] Abstract Besides the fluctuation-dissipation theorem around the stationary state, how the fluctuations in non-stationary states correlate with the future motion of the dynamic variables remains a challenging problem. Further, most temporal correlations of the dynamic variables are zero or very weak in a large class of dynamic systems if the dynamic effects of non-stationary states are not considered. We propose novel methods to compute the temporal correlation functions taking into account the dynamic effects of the non-stationary states. In various dynamic systems, we reveal that the past dynamic fluctuations do drive the future motion of the dynamic variables. This dynamic effect of the non-stationary states is a robust, intrinsic and important property of the complex dynamic systems. Keywords: non-equilibrium statistical mechanics, complex dynamic systems, econophysics Reference [1] T.T. Chen, B. Zheng, Y. Li, and X.F. Jiang, Temporal correlation functions of dynamic systems in non-stationary states, New J. Phys., 2018, 20, 073005

ST-40

Self-consistent phonon theory for thermal transport through anharmonic systems Dahai He Department of Physics, Xiamen University, China

Abstract: Developing a first-principles based approach for classical/quantum thermal transport across low-dimensional systems not only provides insight to potential nano-device applications, but is also crucial to better understand nonequilibrium statistical physics. To date, it is still a challenging problem to deal with thermal transport in a system of strong anharmonicity, which is nevertheless one of the key ingredients for interesting practical applications. In this talk, I will introduce a feasible and effective approach to study quantum thermal transport through anharmonic systems. As its applications, a high-ratio quantum thermal rectification and bi-resonant thermal transport will be shown. References: [1] D. He, J. Thingna, J.-S. Wang and B. Li, Phys. Rev. B 94,155411(2016).

DDAP10, Huaqiao University, Xiamen China, 2018

[2] D. He, J. Thingna, and J. Cao, Phys. Rev. B 97,195437(2018).

ST-41

Novel controls of heat transfer via metamaterials Jiping Huang Department of Physics, Fudan University, Shanghai 200433, China E-mail address: [email protected] Abstract It is known that humans are faced with a global energy crisis, namely, an increasing shortage of nonrenewable energy resources, such as coal, petroleum, and natural gas. However, much of the energy generated from nonrenewable energy resources is changed into thermal energy, which is hard for humans to re-use freely. Therefore, it is meaningful and challenging to manipulate the flow of heat (thermal energy). Here, by establishing temperature-dependent transformation thermotics for treating materials whose conductivity depends on temperature, we show evidence for switchable thermal cloaking and a macroscopic thermal diode based on the cloaking. Meanwhile, we also establish a theory of temperature trapping, and then propose and fabricate a new thermostat for maintaining constant temperatures (within a temperature gradient in the environment) without the need of consuming additional energy. These results suggest that our theory could be adopted for achieving novel macroscopic heat management by using thermal metamaterials, and it could provide new guidance for energy saving.

ST-42

Heat conduction in complex networks Zonghua Liu East China Normal University Email: [email protected] Abstract I will briefly report our recent results of heat conduction in complex networks.

ST-43

Simulating Topological Phase Transition and Quantum Spin/Valley Hall Effects in Linear Circuit Networks Jie Ren Center for Phononics and Thermal Energy Science, China-EU Joint Center for Nanophononics, Shanghai Key Laboratory of Special Artificial Microstructure Materials and Technology, School of Physics Sciences and Engineering, Tongji University, Shanghai, China Abstract Topological states and quantum spin/valley Hall effect (QSHE/VHE) of electrons has improved the development of condensed matter research nowadays, which describe spin-dependent quantum 64

DDAP10, Huaqiao University, Xiamen China, 2018 transport behavior in solid states. Based on metamaterials, the essential physical mechanism behind topological states and QSHE/VHE of classical waves can be understood deeply and verified easily. I will present, for example, 1) the topological Jackiw-Rebby interface mode between permittivity negative and permeability negative metamaterials, which are topologically distinct. The results generally unify the physics of other phononic and elastic metamaterials. 2) QSHE/VHE in topological circuits. The latter includes the topological Lieb lattice and Kane-Mele lattice in a radio-frequency circuit. We design specific capacitor-inductor connections to resemble the intrinsic spin-orbit coupling and construct the analog spin by mixing degrees of freedom of voltages. As such, one is able to simulate the QSHE/VHE in the topological Lieb lattice and Kane-Mele lattice of linear circuits. We then investigate the spin/valley-resolved topological edge mode and the topological phase transition of the circuit varied with capacitances. Our results may find implications in engineering microwave topological metamaterials for signal transmission and energy harvesting. REFERENCES [1]W. Zhu, S. Hou, Y. Long, H. Chen, and J. Ren, “Simulating quantum spin Hall effect in the topological Lieb lattice of a linear circuit network”, Phys. Rev. B 97, 075310 (2018). [2] W. Zhu, Y. Long, H. Chen, and J. Ren, “Realizing quantum spin/valley Hall effects and spin-valley locking in Kane-Mele circuit network”, under review, 2018 [3] Z. Guo, H. Jiang, Y. Long, K. Yu, J. Ren, C. Xue, and H. Chen, “Photonic Spin Hall Effect in Waveguides Composed of Two Types of Single-Negative Metamaterials,” Sci. Rep. 7, 7742 (2017) [4] L. Yang, J. Ren, H. Jiang, Y. Sun, and H. Chen, Acta Physica Sinica 66, 227803 (2017) [5] J. Jiang, W. Zhu, Y. Long, J. Ren, and H. Chen, “Seeing Topological Order and Band Inversion in Optical Diatomic Chain”, under review 2018.

ST-44

Anomalous Energy Diffusion in two-dimensional nonlinear lattices Liu Tian-xin a), Wang Jian a), Li nianbei b) a) College of Physical Science and Technology, Yangzhou University, Yangzhou 225002, China b) Department of Physics, College of Information Science and Engineering, Huaqiao University Abstract Heat conduction in three-dimensional system satisfies the empirical Fourier law, which microscopically implies a normal diffusion of energy. However, in low-dimensional systems with the conservation of momentum, the Fourier law typically does not hold any more such that thermal conductivity 휅 tends to diverge with the system size 퐿 as 휅~퐿훼 for one-dimensional system or 휅~푙푛퐿 for two-dimensional system. It is demonstrated there is a profound connection between anomalous heat conduction and anomalous energy diffusion. Many investigations have been made to explore the relation between thermal transport and energy diffusion in one-dimensional systems, including the mode-coupling theory, the hydrodynamic method and the energy fluctuation correlation function. However, for two-dimensional systems, both an analytical or numerical approach will be hard due to the largely increased complexity. The energy diffusion in two-dimensional system remains elusive. In the work, we investigate the behavior of energy diffusion in two-dimensional nonlinear lattice through GPU-Accelerated computation. We find that the energy diffusion still shows anomalous process when the dimension of lattice varies from one to two dimension. The exponent 훽 for the mean 65

DDAP10, Huaqiao University, Xiamen China, 2018 square displacement (MSD) 〈∆푟2(푡)〉 ∝ 푡훽 changes from about 1.4 for one-dimension to 1.27 for two-dimension. In contrast with the anomalous energy diffusion, the momentum demonstrates a normal diffusion process. Our findings contribute to understand thermal transport in two-dimensional system such as graphene and silicene, which have been received attention experiments.

Fig. 1 The mean square displacement (MSD) 〈∆푟2(푡)〉 ∝ 푡훽 changes with the increasing width of the nonlinear lattice.

Fig. 2 The heights of momentum peak changes with the increasing width of the nonlinear lattice.

ST-45

Using Hilbert transform and classical chains to simulate quantum walks Daxing XiongA, Felix ThielB and Eli BarkaiB ADepartment of Physics, Fuzhou University, Fuzhou, Fujian, 350108, China BDepartment of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan, 52900, Israel email: phyxiongdxfzu.edu.cn Abstract In this work we report a classical simulation strategy that we have proposed recently used to simulate quantum dynamics, in particular quantum walks. This strategy is based on a classical device of linearly coupled chain of springs together with the technique of Hilbert transform. Through this strategy, we obtain the quantum wave function from the classical evolution. Specially,

DDAP10, Huaqiao University, Xiamen China, 2018 this goal is achieved with the classical momenta of the particles on the chain and their Hilbert transform, from which we construct the many-body momentum and Hilbert transformed momentum pair correlation functions yielding the real and imaginary parts of the wave function, respectively. With such a wave function, we show that the classical chain’s energy and heat spreading densities can be related to the wave function’s modulus square. This relation provides a particular perspective from quantum aspects to understand ballistic energy and heat transport. References D. Xiong, F. Thiel, and E. Barkai, Phys. Rev. E 93, 022114 (2017)

ST-46

One-dimensional superdiffusive heat propagation induced by optical phonon-phonon interactions Yong Zhang, Xiamen University Abstract One-dimensional anomalous heat propagation is usually characterized by a Levy walk superdiffusive spreading function with two side peaks located on the fronts due to the finite velocity of acoustic phonons. In the case when the acoustic phonons vanish, e.g., due to the phonon-lattice interactions such that the system’s momentum is not conserved, the side peaks will disappear and a normal Gaussian diffusive heat-propagating behavior will be observed. Here we show that there exists another new type of superdiffusive, non-Gaussian heat propagation but without side peaks in a typical nonacoustic, momentum-nonconserving system. It implies that thermal transport in this system disobeys the Fourier law, in clear contrast with the existing theoretical predictions. The underlying mechanism is related to an effect of optical phonon-phonon interactions. These findings may open a new avenue for further exploring thermal transport in low dimensions.

ST-47

Compressing binary crystals into glasses Yilong Han, Hong Kong University of Science and Technology [email protected] Abstract TBA.

ST-48

Dynamics of granular rafts on Faraday waves Jilin Jou and Lin I Department of Physics, National Central University, Jhongli Abstract Particles floating on the liquid interface can aggregate and form a granular raft, which affects the interface property and induces several interesting phenomena, like wrinkle pattern or encapsulation 67

DDAP10, Huaqiao University, Xiamen China, 2018 process. However, for particles on the wavy water surface, such as marine debris on a sea surface, the response of wave-particle mutual interaction is still elusive. In this work, we experimentally investigate the spatiotemporal evolution of granular rafts on the surfaces of Faraday waves and the interactions between rafts and waves. Faraday wave is a surface wave parametrically generated by vertically oscillating a tank of fluid. It has been used as a model nonlinear system for investigating many interesting nonlinear wave phenomena such as pattern formation, defect generation, rogue wave, and wave turbulence. Oscillons, localized clumps oscillating with two times of the external driving period, can be formed. Under weak driving, oscillons can be self-organized into ordered patterns such as triangular and square lattice like patterns. With increasing driving and through modulation instability, the lattice like ordered oscillon pattern is deteriorated, to the intermediate state with phase and amplitude modulations of the oscillon pattern, before the transition to the disordered turbulence state. In our work, the raft and raft-wave interaction dynamics in the ordered state and the intermediate state are investigated. Styrofoam particles 1 mm in diameter are sprayed on the surfaces of a Faraday wave of a water tank 18 cm in diameter and 2 cm in depth, vertically driven by a speaker oscillating at 40 Hz. The spatial distribution (in the horizontal xy-plane) of the intensity of light, generated from the bottom uniform LED array, passing through a thin bottom glass diffuser plate and then attenuated by the milk particles in the fluid, is measured by a CCD/telescope system. From the spatial distribution of local light attenuation normalized by the local light intensity without oscillating driving, the spatiotemporal variations of water surface height can be determined. The motion of each particle on water surfaces are digitally tracked from the sequential CCD images. The raft size is changed from 1 to 8 mean inter-oscillon distance. The gravitational force of particles on the wave surface with different tilting angle, the friction between particles and the flowing fluid elements which constitute oscillons are the key factors affecting raft-wave dynamics. It is found that, increasing raft size suppresses underlying wave height and affects the topologies of the surrounding oscillon patterns, which in turn affects the raft dynamics, e.g. inducing splitting or recombination of rafts, causing anomalous diffusions of the translational and rotation motions of rafts. On the ordered oscillon surfaces, raft tends to split along (1, 1) and (1, -1) directions. Increasing wave disordering enhances oscillon height fluctuation and oscillon anomalous diffusion, which enhance the translational and rotational diffusions of rafts. Keywords: granular raft, Faraday wave, wave-particle interaction, cage-hopping motion

ST-49

Reinforcement Learning for Complementarity Game and Population Dynamics Wei Li Central China Normal University [email protected] Abstract We systematically test and compare different reinforcement learning schemes in a complementarity game [J. Jost and W. Li, Physica A 345, 245 (2005)] played between members of two populations. More precisely, we study the Roth-Erev, Bush-Mosteller, and SoftMax reinforcement learning schemes. A modified version of Roth-Erev with a power exponent of 1.5, as opposed to 1 in 68

DDAP10, Huaqiao University, Xiamen China, 2018 the standard version, performs best. We also compare these reinforcement learning strategies with evolutionary schemes. This gives insight into aspects like the issue of quick adaptation as opposed to systematic exploration or the role of learning rates.

ST-50

Contrasting complexity of the recent 2000 and past 122,000 years climates Zhi-Gang Shao Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, SPTE science South China Normal University, Guangzhou, China

Abstract Understanding natural climate variability is essential for assessments of climate change. Paleoclimatic changes through the history of Earth can be evidenced in the high temporal resolution Oxygen isotopic (δ18O) records from ice cores, where the records present a proxy for the temperature. Based on the δ18O records of ice cores, the complexity of the climate of the past 122,000 years and recent 2000 years was investigated by the sample entropy (SampEn) method and Lempel-Ziv (LZ) complexity. In using SampEn method, the climate complexity is measured by the sample entropy, which is a modified approximate entropy defined in terms of the occurring probability of new modes in a record. A larger sample entropy reflects a higher probability to spot a new mode in the data, and in this sense signals a larger complexity of the sample. The δ18O record of the past 122,000-year is found to have smaller SampEn than the recent 2000-year. This result suggests that the climate of the past 122,000-year has less complexity than that of the recent 2000 years, even though the record for the former exhibits stronger fluctuations and multifractality than the latter. This diagnosis is additionally supported by calculations of LZ complexity, which has smaller value for the record of the past 122,000 years than the recent 2000 years. Our theoretical findings may further contribute to ongoing explorations into the nonlinear statistical character of the climate change. Keywords: Climate, Ice core, Isotopic, Complexity Acknowledgments: The National Natural Science Foundation of China (Grant No. 11105054) and the Natural Science Foundation of Guangdong Province (Grant No. 2018A030313322) References: [1] Z.-G. Shao, Contrasting the complexity of the climate of the past 122,000 years and recent 2000 years, Sci Rep 2017, 7, 4143. [2] Z.-G. Shao, & P. D. Ditlevsen. Contrasting scaling properties of interglacial and glacial climates. Nature communications 7:10951 (2016). [3] Z.-G. Shao, & H.-H. Wang, MultifractalDetrended Fluctuation Analysis of the δ18O record of NGRIP ice core, Climate Dynamics 43(7-8), 2105-2109 (2014)

DDAP10, Huaqiao University, Xiamen China, 2018

ST-51

Universal scaling laws of growth processes and their implications for chaos instability Yohsuke T. Fukai1 and Kazumasa A. Takeuchi1,2 1 Department of Physics, The University of Tokyo. 2 Department of Physics, Tokyo Institute of Technology. Abstract Random growth of interfaces is a phenomenon often encountered in systems driven out of equilibrium, which may show universal fluctuation properties classified into universality classes. Among them, the Kardar-Parisi-Zhang (KPZ) universality class is known to arise generically and plays a distinguished role, all the more since many of its properties in one dimension turned out to be exactly solvable [1]. An important outcome of this development is that KPZ interface fluctuations are universal yet dependent on the interface geometry, or equivalently the initial condition, being classified into a few universality “subclasses”. For example, the distribution function was shown to be different between the circular and flat interfaces [1,2] (see figure below). In the talk, the speaker will first overview the geometry-dependent KPZ scaling laws, on the basis of liquid-crystal experiments he and coworkers carried out [1,2]. He will then discuss their implications for dynamical systems theory. Specifically, it has been shown that growing perturbations in spatiotemporal chaos are governed by the KPZ-class exponents [3]. Here we show that growth of the perturbations is also controlled by the KPZ geometry dependence, being determined by the initial form of perturbations [4]. This can be used for example for precise evaluation of the Lyapunov exponent, which is usually hindered by slow convergence specific to spatiotemporal chaos.

Figure: Circular and flat interfaces of growing liquid-crystal turbulence [2]. The interface height h shows the different distributions characteristic of the two geometries, namely the largest-eigenvalue distribution of GUE and GOE random matrices for the circular and flat case, respectively. References [1] For reviews, see, e.g., K. A. Takeuchi, Physica A 504, 77 (2018); I. Corwin, Random Matrices Theory Appl. 1, 1130001 (2012); T. Kriecherbauer and J. Krug, J. Phys. A 43, 403001 (2010). [2] K. A. Takeuchi et al., Sci. Rep. 1, 34 (2011); K. A. Takeuchi and M. Sano, J. Stat. Phys. 147, 853 (2012). [3] A. S. Pikovsky and J. Kurths, Phys. Rev. E, 49, 898 (1994). [4] Y. T. Fukai and K. A. Takeuchi, to be published.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-52

A robust collective credit allocation in scientific publications Fenghua Wang，An Zeng，Xiaohua Cui，Ying Fan，Zengru Di School of Systems Science, Beijing Normal University, 100875 Beijing, China Abstract Collaboration among researchers plays an important role in scientific discoveries, especially in multidisciplinary research. How to allocate credit reasonably to coauthors of a paper is a long-standing problem. The collective credit allocation method which evaluates coauthors' contribution to papers based on network diffusion process provides a new way of solving the problem. Nevertheless, the existing collective allocation methods are sensitive to manipulations. In this article, we propose a robust credit allocation method (RCCA method) considering citing papers' scientific impact based on the state-of-the-art collective credit allocation method (CCA method). We have three steps to validate RCCA method. Firstly, we apply it to Nobel papers to identify Nobel laureates. Then we calculate and compare the number of artificially added papers to make credit share of the original lowest credit share author exceed that of the original largest credit share author in Nobel papers with two methods. Finally, we choose three kinds of papers whose citation range are 0-10, 10-100,100-1000 to testify whether RCCA method is effective to ordinary papers by calculating the number of added papers which make credit share of the original lowest credit share author exceed that of the original largest credit share author. Compared to the existing method, we find that the identification accuracy of Nobel laureates of our method outperforms CCA method. What's more, we find our method can restrain well the influence of credit share by low scientific impact citing papers through applying our method to the Nobel papers dataset and to ordinary papers. It indicates that RCCA method introduced in this paper is more robust comparing to CCA method, because more added papers are required for the low credit share authors to raise their credit share in short term in RCCA method. The parameter of RCCA method can adjust the effect of credit share caused by citing papers in the algorithm. The larger parameter is, the larger high scientific impact citing papers having an effect on credit share is. If the parameter takes an infinite value, which means the credit allocation of a given paper is only determined by the highest impact citing paper. Keywords: Citation network; Credit allocation

ST-53

Qusicrystal formation in monodisperse soft-core particle systems Mengjie Zu1, Peng Tan2, Ning Xu1* 1Hefei National Laboratory for Physical Sciences at the Microscale,CAS Key Laboratory of Soft Matter Chemistry, and Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China 2 State Key Laboratory of Surface Physics and Department of Physics, Fudan University, Shanghai 200433, P. R. China Abstract For soft-core particle systems interacting via a spring-like repulsion, we surprisingly observe the existence of two-dimensional quasicrystals at high densities. Up to now, quasicrystals have been found 71

DDAP10, Huaqiao University, Xiamen China, 2018 in multi-component systems such as alloys and in mono-component systems with multiple length scales in particle interactions or with anisotropic particles such as tetrahedral and patchy particles. In all these systems, multiple length scales are present in particle size, interaction, or anisotropy, which seem to be the consensual condition for quasicrystals to occur. The quasicrystals that we found are unexpectedly formed by monodisperse, isotropic particles interacting via the simple soft-core potential without multiple length scales. We find not only dodecagonal but also octagonal quasicrystals, which have not been found yet in soft quasicrystals. In our quasicrystals, particles tend to form pentagons. A pentagon surrounded by an n-fold polygon constructs the complex structural unit to develop the n-fold quasicrystalline order. In liquid states prior to the liquid-solid transition, we already observe some signs of the development of the quasicrystalline order, including the accumulation of pentagons and the emergence of two competing length scales. We have also verified that our quasicrystals are stable by showing that their formation is independent of history and they have the lowest inherent structure potential energy, compared with other crystalline states. Key Words: Quasicrystals, soft-core, structure, length scale

ST-54

Predicting discredited behaviors of millions of enterprises Tao Zhou Big Data Research Center, University of Electronic Science and Technology of China Abstract Characterizing enterprise credit level and predicting various discredited behaviors play a key role in business cooperation, equity investment, guaranteed loans, investment attraction and other related applications, and are also an indispensable part in building a credit society and improving the business environment. Previous analysis mainly focused on the characteristics of enterprise size, place of operation, industry category, registration and paid-in capital, and lacked in-depth analysis based on massive data. This work builds up a directed investment network including more than 5 million enterprises, among which more than 6% of enterprises have various discredited behaviors. The results show that there is an obvious "network effect" in the discredited behaviors of enterprises. If the target enterprise's shareholders or its invested enterprises have discredited behaviors, the risk of having discredited behaviors of the target enterprise is far greater than the average. Based on this, a simple generalized linear regression algorithm is proposed to predict the discredited behaviors of enterprises, which is far more accurate than the regression method without considering the network effect. Keywords: Discredited Behaviors; Network Effect; Generalized Linear Regression.

ST-55

The impact of selfish route decisions on optimized transportation networks Po Ho Fai1, Chi Ho Yeung1 and David Saad2 1Department of Science and Environmental Studies, The Education University of Hong Kong 2The Non-linearity and Complexity Research Group, Aston University, Birmingham, United Kingdom Abstract Optimizing traffic flow in transportation networks is one way to tackle traffic congestion. 72

DDAP10, Huaqiao University, Xiamen China, 2018

However, even when optimally coordinated routing suggestions have been identified and recommended for each individual road user, some may not follow these recommendations as the suggested routes may incur a larger traffic cost to individuals in order to achieve a lower global cost. Global optimum cannot be achieved when some users selfishly choose their own routes to minimize individual costs. In this work, we analyze a transportation network model where each road user is given a recommended path, but some take a selfish routing strategy to minimize their own costs. Specifically, we employ cavity approach developed in the area of statistical physics to analyze the user behavior in two scenarios of recommended traffic flows, namely a scenario with all users following the shortest path, and another scenario in which their recommended routes are optimally coordinated. In both scenarios, some users choose to optimize their own cost assuming all other users follow the recommended paths. Our analytical and simulation results show that when users are originally given the optimally coordinated routes, the global cost increases only slightly when the number of selfish road users increases. On the other hand, if users originally follow their own shortest path to the destinations, the fractional change of the global cost after taking selfish route decisions is a convex function which first decreases to negative values and then increases back to positive (higher cost), implying that the system benefits when there is only a small number of selfish users but suffers when the number of selfish road users further increases. A phase diagram is drawn which indicates the gain and the loss in this case. Our analytical and simulation results consider users who take their selfish decisions only once, when a recommended route configuration is given; in reality, users continue to take multiple rounds of selfish decisions after the first one. To understand the consequence of multiple rounds of selfish route decisions, we conduct simulations to reveal an interesting dynamical behavior and the potential Nash equilibrium, which may result after multiple rounds of selfish routing decisions. Key Words: Optimization, Selfish Routing, Nash Equilibrium, Transportation network

ST-56

NES-TL: Network Embedding Similarity-Based Transfer Learning Chenbo Fu, Zhejiang University of Technology [email protected] Abstract The transfer learning methodology leverages knowledge from the source domain with fruitful training data to the barren target domain. Recently, new approaches continue to be developed and used to solve different classification tasks, ranging from public news to videos and to many others. Most transfer learning methods are based on the assumption that both training and test data are in the same feature space with the same data distribution, which however is not always true in real applications where it would lead to a negative transfer. In order to overcome this hurdle, the multiple-source transfer learning framework is useful. Since many real systems can be represented by networks, how to utilize the structural similarity between different networks so as to increase the transfer effectiveness becomes important. In this work, the NES specification index is used to quantitatively measure the structural similarity between two networks, based on which a new transfer learning method (named NES-TL) is developed. Experiments on tag popularity prediction in Stack Exchange QA communities verify the effectiveness of the proposed approach, showing that it behaves better than existing baseline methods.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-57

Complex networks and universal patterns in historical Chinese dynasties Chin-Kun Hu Department of Physics, National Dong Hwa University, Hualien 97401 Institute of Physics of Academia Sinica, Taipei 11529 http://proj1.sinica.edu.tw/~statphys/ Abstract Critical physical systems with large numbers of molecules or atoms can show universal and scaling behaviors [1]. It is of interest to know whether human societies with large numbers of people can show universal and scaling behaviors. Here we use network theory [2-4] to analyze the Chinese history in periods 209BC-23 AD and 532AD-619 AD related to Western Han-Xin Mang and Northern Zhou-Sui Dynasty, respectively. The data for the former period were taken from Shi Ji [5] and Han Shu. The data for the later period were taken from Bei Shi and Sui Shu. We find that the networks from two periods separated about 500 years have interesting universal and scaling behaviors [6]. We have also found some interesting universal patterns in historical Chinese dynasties. [1] Chin-Kun Hu, Chinese Journal of Physics 52, 1-76 (2014). [2] R. Albert and A.L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). [3] E. Ravasz and A.L. Barabasi, Phys. Rev. E 67, 026112 (2003) [4] Liming Pan, Tao Zhou, Linyuan Lu and Chin-Kun Hu, Scientific Reports, 6, 22955 (2016). [5] Ting-Ting Chi, Shih-Chieh Wang and Chin-Kun Hu, Physica A, 437, 408 (2015). [6] Yi-Ru Pan, Wen-Jong Ma, Chen-Yu Chang, Hui-Jie Yang ,Ting-Ting Chi, and Chin-Kun Hu, preprint.

ST-58

Understanding the temporal pattern of spreading in heterogeneous networks: Time as a variable Mi Jin Lee1 and Deok-Sun Lee1, ∗ 1Department of Physics, Inha University, Incheon 22212, Korea Abstract Towards a reliable prediction of epidemic or information spreading pattern in complex systems, well-defined measures are essential. In the susceptible-infected (SI) model on heterogeneous networks, the cluster of infected nodes in the intermediate-time regime exhibits too large fluctuation in size to use its mean size as a representative value. We show that the infection of hub nodes is so crucial for global spreading that the number of infected nodes at a given time strongly fluctuates depending on when a hub node was first infected, resulting in a very broad distribution of the cluster size. On the contrary, the distribution of the time taken to infect a given number of nodes is well concentrated at its mean, suggesting the mean infection time to be a better measure. We present a theory for evaluating the mean infection time by using the boundary area of the infected cluster of a given size, and use it to find a non-exponential but algebraic spreading of infection with respect to time in the intermediate stage on strongly heterogeneous networks. Such slow spreading originates in only small-degree nodes left

DDAP10, Huaqiao University, Xiamen China, 2018 susceptible while most hub nodes already infected in the early exponential-spreading stage. Our results offer a way to detour around large statistical fluctuations and quantify reliably the temporal pattern of spread under structural heterogeneity.

ST-59

Temporal similarity metrics for network reconstruction based on diffusion dynamics Ming-Kai Liua, Hao Liaoa*, Manuel Sebastian Marianib,c, Mingyang Zhoua, Xingtong Wua aNational Engineering Laboratory on Big data Application on Improving Government Governance Capabilities, Guangdong Province Key Laboratory of Popular High Performance Computers, College of Computer Science and Software Engineering, Shenzhen University, Shenzhen 518060, PR China. bInstitute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, 610051 Chengdu, People's Republic of China. cURPP Social Networks, Universitat Zurich, CH-8050 Switzerland. Abstract: When investigating the spreading of a piece of information or the diffusion of an innovation, we often lack information on the underlying propagation network. Reconstructing the hidden propagation paths based on the observed diffusion process is a challenging problem which has recently attracted attention from diverse research fields. To address this reconstruction problem, based on static similarity metrics commonly used in the link prediction literature, we introduce new pairwise temporal metrics of similarity between the nodes that take as input the time-series of multiple independent spreading processes. The best-performing similarity metrics take into account both structural and temporal information, based on the hypothesis that two nodes are more likely to be connected if they were often infected at similar points in time. Besides, we find that the best reconstruction accuracy is achieved by a temporal generalization of the Cosine similarity metric. Our findings provide an extensive benchmark of temporal similarity metrics for network reconstruction, and shed new light on the notion of similarity between pairs of nodes in complex networks. Key word: Social networks, Network reconstruction, Temporal similarity, Innovation diffusion

ST-60

Vulnerability and Co-susceptibility Determine Large Network Cascades Takashi Nishikawa Northwestern University Abstract In a network, a local disturbance can propagate and eventually cause a substantial part of the system to fail, in cascade events that are easy to conceptualize but extraordinarily difficult to predict. Here, we first develop a statistical framework that can predict cascade size distributions by incorporating two ingredients only: the vulnerability of individual components and the co-susceptibility of groups of components (i.e., their tendency to fail together). Using the North American power grid as a complex system per excellence, we show that correlations between component failures define 75

DDAP10, Huaqiao University, Xiamen China, 2018 structured and often surprisingly large groups of co-susceptible components. We then identify, quantify, and analyze the set of network components that are vulnerable to cascading failures across multiple conditions in which the system operates. We show that the vulnerable set consists of a small but topologically central portion of the network and that large cascades are disproportionately more likely to be triggered by initial failures close to this set. Aside from their implications for blackout studies, these results provide insights and a new modeling framework for understanding cascades in financial systems, food webs, and complex networks in general.

ST-61

The Benefits of Probabilistic Route Choices in Dynamical Transportation Networks Tai Tak Shing and Yeung Chi Ho Department of Science and Environmental Studies, The Education University of Hong Kong Tai Po, Hong Kong, China Abstract With the increasing number of road vehicles, the transportation networks become more crowded and traffic congestion occur. Most vehicles travel on the shortest paths between the origin and destination of their trips, which are supposed to be the fastest way to go to the destination, but traffic congestion occur when too many shortest paths lie on the major common routes. To further understand traffic congestion, we studied a model of cellular automata in which vehicles probabilistically choose between a move on the shortest path or a move on a random route to their destinations. Specifically, we implement the model of cellular automata on two-dimensional square lattices. In the model, each vehicle starts at a randomly draw origin and travels to a randomly drawn destination. When the vehicle arrives at its destination, a new destination is randomly drawn. On the way to their destinations, vehicles move in four directions (i.e. upward, downward, left and right) and are only probabilistically biased towards the shortest paths to the destinations. We call this probability for vehicles to be on the shortest paths the “path-greediness”, as it represents the tendency for road users to choose the shortest paths on the destination. Our results show that the changes in the path-greediness affect the flow of the traffic. The cases with less path-greediness mean that the road users tend to choose longer paths, and the system transit to congested state at a higher density of cars than the case with high path-greediness. We further compare the cases with a given global path-greediness and an individual adaptive path-greediness to study the impacts by path-greediness, respectively representing scenarios with a centralized control and self-organizations. Our results show that the cases with individual adaptive path-greediness has less congestion than the cases with global path-greediness, but the cases with global path-greediness has less congestions than the cases with individual adaptive path-greediness in when vehicle density is high. It implies that the optimal traffic distribution is dependent on vehicle density. The study also shows if some cars travel longer of choose routes probabilistically, traffic congestion are reduced though the total travelling distance my increase for individual cars. Keywords: Routing, Transportation Networks, Cellular Automata Email address: [email protected], [email protected]

DDAP10, Huaqiao University, Xiamen China, 2018

ST-62

Impacts of mutual selection in random walk process in time-varying networks Mei Yang, Bing Wang*, YueXing Han School of Computer Engineering and Science, Shanghai University, No. 99, Shangda Rd., Baoshan Dis., Shanghai, P.R. China Abstract Random walks are one of the most fundamental types of stochastic processes and have been applied in various domains, such as ranking systems and searching. In this work, we investigate random walk process on a class of temporal networks by considering the mutual selection principle. The mutual selection model has been fully investigated in a static form, however, the study of the mutual selection network model in temporal form, especially for random walks has not been fully investigated. The fundamental factor in the network model is characterized by a linking function that describes the probability of the existence of an edge, which depends mutually on the fitness of the vertices on both ends of that edge. Here, we investigate two typical forms of linking functions. The first case we investigate is based on the assumption that the fitness of the ends on a link is independent, while the other case is assumed that the creation of the link is determined by the sum of the two ends’ fitness. If the sum is larger than a threshold, then the two nodes will connect, otherwise not. We find that the mutual agreement in this temporal network model has different effects on the properties of random walk processes, which is totally different from that obtained in the activity driven model. In the first case, we find that the capability of vertices to gather walkers is not only related to the vertices' activity, but also determined by their propensity of receiving connections. For the second case, we find that for vertices with activity larger than a given threshold. The stronger the activity is, the more walkers they will collect in the stationary state. We confirm our analytical prediction via large-scale numerical simulations. The results presented here contribute to the understanding of the evolution mechanism of this type of temporal network model and give us insight for their effects on random walks processes. Key Words: random walks, mutual selection, temporal networks

ST-63

Relation between network structure and controllability transition in complex networks Bing-Hong Wang (汪秉宏) Department of Modern Physics, University of Science and Technology of China, Hefei , 230026, P.R.China Abstract The network control problem has been attracted increasing attention from avoiding cascading failures of power grids to managing ecological networks. It is proved that the numerical control can be achieved if the number of control inputs exceeds the transition point. Here we investigated the effect of degree correlation on numerical controllability of networks in both real networks and modeling networks. Finding that the transition point of the number of control inputs greatly depends on degree correlation in sparse undirected networks and there exists a minimum value of the transition point with 77

DDAP10, Huaqiao University, Xiamen China, 2018 increasing of degree correlation. More interestingly, the effect of degree correlation on the transition point cannot be observed in dense networks for numerical controllability, which is contrary to the result of structural controllability. In the condition of directed random networks and scale-free networks, the influence of degree correlations is determined by the types of correlations. Our approach provides the understanding of control problems in complex sparse networks.

ST-64

Self-organized criticality of earthquake model on networks with randomness Sheng-Jun Wang Department of Physics, Shaanxi Normal University, Xi’An Abstract We study the nonconservative earthquake model on a random spatial network. The spatial networks are composed by sites on a 2D plane which are connected locally. Different from the regular lattice, the randomness of sites are modelled in the way that sites are randomly placed on the plane. Use the same connectivity degree as the 2D lattice, however, the spatial network cannot exhibit critical earthquake behavior. Mimicking the long range energy transfer, the connection radius is increased and the connectivity degree of the spatial network is increased. Then we show that the model exhibits self-organized criticality. The mechanism of the structural effect is presented. The spatial network includes many modules when connectivity degree is very small. The effect of modular structure on the avalanche dynamics is to limit the spreading of avalanches in the whole network. When the connectivity degree is larger, the long range energy transfer can overcome the effect of local modularity and the criticality can be reached. Key Words: self-organized criticality, earthquake, network

ST-65

Identifying the scale-free property in real networks Gang Yan School of Physics Science and Engineering, Tongji University, 1239 Siping Road, Yangpu, Shanghai, 200092 Abstract The scale-free property, i.e. node degree distribution follows a power law, has been one of the most important discoveries in the field of network science. Yet, there was a debate very recently on whether scale-free networks widely exist, prompting us to revisit the identification of the scale-free property in real networks. In this talk, we will show that the existing methods is not able to distinguish lower law distributions and lognormal distributions, and present a new method of identifying the scale-free property. We validate our method in both model networks and a variety of real networks, finding that scale-free networks do exist widely.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-66

Evolutionary Evidences in Structure of Complex Network Shiguo Deng1, Jin Zhang1,2, and Huijie Yang1 1Department of Systems Science, University of Shanghai for Science and Technology, Shanghai 200093, China 2 School of Science, Jinan University, Jinan 250022, Shandong Province, China Abstract Most proteins of an organism function through a complex web of interactions with other proteins, which is known as a protein interaction network (PIN). Research on the evolution of PINs not only sheds light on the principles driving the evolution of living organisms but also helps us understand the biological organization and function of cellular proteins, just as Dobzhansky stated that ‘nothing in biology makes sense except in light of evolution. In one of our works, we proposed a method to improve the accuracy of age predicting by embedding multisource biological information in each iteration of an archaeology algorithm for yeast PIN. We reduce the probability of reversing errors by decreasing the non-duplication protein pairs, which are obtained from 460 gene trees constructed by means of a multiple sequence alignment and the neighbor joining algorithm. The reliable crossover standard from different biological information sources is also used to decrease local random errors of alternative treatment. Our research strongly suggests that integrating non-biological restrictions with biological context will bear more favorable results. In another of our works, we introduce a new concept called reducibility of complex network to find historical evidences in network structures. In many evolutionary networks, redundancy is produced by duplication and mutation rules. For instances, in a protein-protein interaction network, a duplication of gene may produce an exact copy of an existing protein. These two identical proteins share the same set of nodes as their neighbors. At the same time, mutation of a gene will change its corresponding protein, which may convert two identical nodes to two distinguishable nodes. Some existing linkages may be broken, while some new linkages may be constructed. In an academic cooperation network, at the beginning a new scientist may inherit relationships from its supervisor. Due to some reasons such as change of interesting, he/she may construct new linkages and lose some existing linkages with other people. These two competitive factors lead to complicated redundant structures in the resulting networks. Hence, redundancy of a network is closely related with its evolutionary behaviors. Rich information on evolutionary behaviors are embedded in redundancy of networks, but how to measure quantitatively redundancy of a network and to discover evolutionary clues from redundancy are still open problems. In the present paper, a new concept, called reducibility of network, is introduced to evaluate in a global way redundant degree of a network. The basic idea is to regard all the nodes as a set of variables describing a complex system and each column of adjacency matrix as a set of measurement records for the corresponding variable in different experimental cases. A copy of a node will have the same records with the template node. Mutation will leads to difference between the records of the two nodes. Hence, the covariance matrix for all the nodes stores the information on competition of duplication and mutation. Technically, ranking the eigenvalues of the covariance matrix in a descending order, the ratio of summation of a specified number of eigenvalues over that of the total spectrum can be employed to measure quantitatively the reducibility.

DDAP10, Huaqiao University, Xiamen China, 2018

This concept detects successfully evolutionary information stored in the empirical human metabolic network and the modeling protein-protein interaction networks, which may be extend further to other problems such as finding leadership groups in a social network and searching critical communities in a cellular network. Keywords: Reducibility of complex network; Evolutionary evidence Reference [1] J Zhang, H J Yang, et al., IEEE Acess 5,15893-15900(2017). [2] S G Deng, H J Yang, et al., Reducibility of Complex networks, submitted. [3] S. Navlakha, and C.Kingsford, PLoS Comput. Biol. 7(4):e1001119(2011).

ST-67

Dynamics of Epidemic Spreading and Information Diffusion on Social Networks Xiu-Xiu Zhan, Chuang Liu, Zi-Ke Zhang*,Gui-Quan Sun, Jonathan J.H. Zhu, Zhen Jin Alibaba Research Center for Complexity Sciences, Hangzhou Normal University, Hagnzhou 311121 Email: [email protected] Abstract Recently, the coupling effect of information diffusion (or awareness) and epidemic spreading has facilitated an interdisciplinary research area. When a disease begins to spread in the physical society, the corresponding information would also be transmitted from among individuals, which in turn influence the spreading pattern of the disease. In this paper, we have studied the coupling dynamics between epidemic spreading and relevant information diffusion. Empirical analyses from representative diseases (H7N9 and Dengue fever) show that the two kinds of dynamics could significantly influence each other. In addition, we propose a nonlinear model to describe such coupling dynamics based on the SIS (Susceptible-Infected-Susceptible) process. Both simulation results and theoretical analyses show the underlying coupling phenomenon. That is to say, a high prevalence of epidemic will lead to a slow information decay, consequently resulting in a high infected level, which shall in turn prevent the epidemic spreading. Further theoretical analysis demonstrates that a multi-outbreak phenomenon emerges via the effect of coupling dynamics, which finds good agreement with empirical results. The findings of this work may have various applications of network dynamics. For example, as it has proved that preventive behaviors introduced by disease information can significantly inhibit the epidemic spreads, and information diffusion can be utilized as a complementary measure to efficiently control epidemics. Therefore, the government should make an effort to maintain the public awareness, especially during the harmonious periods when the epidemic seems to be under control. In addition, in this work, we only consider the general preventive behavioral response of crowd. However, the dynamics of an epidemic may be very different due to the behavioral responses of people, such as adaptive process, migration, vaccination, and immunity. This work just provides a start point to understand the coupling effect between the two spreading processes, a more comprehensive and in-depth study of personalized preventive behavioral responses shall need further efforts to discover.

DDAP10, Huaqiao University, Xiamen China, 2018

Figure 1. Empirical analysis of the epidemic spreading (blue circles) and the information diffusion (pink diamonds) for: (a) H7N9; (b) Dengue fever. (b1) and (b2) are details of partial enlargement of dengue fever).

Figure 2. Evolution pattern of the density of informed and infected with different values of the infected probability β for unaware individuals. (a) β=0.05; (b) β=0.2; (c) β=0.8. The bottom is an interval indication Of β, which corresponds to three phases of informed level: (i) healthy state for 0 ≤β<0:092; (ii) oscillatory state for 0:092 <β ≤ 0:239; (iii) unimodal for 0:239 <β≤ 1. The result is obtained by 10,000 independent realizations.

ST-68

Principle of corresponding states of quasi-particle gases for American passenger flights Chenping Zhu Nanjing University of Aeronautics and Astronautics Abstract The present paper describes a non-equilibrium process with a state equation in equilibrium state, which is realized by mining big data of departure/arrival records of domestic passenger flights of whole the United States for 20 years. We set up a quasi-particle gas model for passenger flights, put temporal effect into phenomenological parameters, and erect a common gaseous constant R for the species of all the years. Such state equations are reduced into a universal form by reducing quasi-thermodynamic quantities (\mu, V, T) with the critical values respectively. In this way, we extend the principle of corresponding states in thermodynamics outside the realm of physics. A practical strategy to enhance the operational efficiency of passenger flights emerges quantitatively: to traverse the valley of universal equation t'(\mu') by increasing the number of Arrival-Departure type of passenger flights, where t' and \mu' are reduced temperature and reduced mole number, respectively. 81

DDAP10, Huaqiao University, Xiamen China, 2018

This is the first work to describe the states of practical domestic passenger flights with large scale long term big data. It presents the first state equation for passenger flights based on big data. The first time we extend the principle of corresponding states outside traditional physics.

ST-69

The formation of continuous opinion dynamics based on a gambling mechanism Yueying Zhu, Qiuping A. Wang, Wei Li, and Xu Cai No. 1, Yangguang Avenue, Jiangxia district, Wuhan, China, 430020. Abstract The formation of continuous opinion dynamics is investigated based on a virtual gambling mechanism where agents fight for a limited resource. We propose a model with agents holding opinions between −1 and 1. Agents are segregated into two cliques according to the sign of their opinions. Local communication happens only when the opinion distance between corresponding agents is no larger than a pre-defined confidence threshold. Theoretical analysis based on the mean field approximation, regarding special cases, provides a deep understanding of the roles of both the resource allocation parameter and confidence threshold in the formation of opinion dynamics. Results provide that, when the system reaches a steady state, the distribution of positive opinions and that of negative ones will gradually tend to be symmetric as resource allocation parameter increases. Furthermore, the increase of confidence threshold will lead the underlying system more close to a consensus state (all agents hold the same opinion). For a sparse network (the communications between agents are quite limited), the evolution of opinion dynamics is negligible in the region of low confidence threshold when the mindless agents are absent. Numerical results also imply that, in the presence of economic agents, a very high confidence threshold is required for apparent clustering of agents in opinion. Moreover, a consensus state is generated only when the following three conditions are satisfied simultaneously: mindless agents are absent, the resource is concentrated in one clique, and confidence threshold tends to a critical value(= 1.25 + 2/ka; ka > 8/3, the average number of friends of individual agents). For fixed confidence threshold and resource allocation parameter, the most chaotic steady state of opinion dynamics happens when the fraction of mindless agents is about 0.7. We also study the probability of both economic and mindless agents to win in the gambling mechanism. Numerical results demonstrate that economic agents are more likely to win at gambling, compared to mindless ones. This is highly consistent with the actual phenomenon in social system. Finally, the global uncertainty and sensitivity analysis is introduced, helping understand the uncertainty of characteristics of opinion dynamics controlled by the underlying system, and also the importance of roles acted by three involved parameters. In terms of Latin hypercube sampling-based method, we quantify the uncertainty of three different observations, including the order parameter of system, dispersion index of agents clustering in opinion, and the fraction of winners in gambling mechanism. The importance of three involved parameters is also estimated, in establishing the uncertainty of three observations. Key Words: Agent-based models, Algorithmic game theory, Evolutionary game theory, Nonlinear dynamics, Uncertainty analysis, Sensitivity index.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-70

Energetic costs, precision, and transport efficiency of molecular motors Hyeon Changbong Korea Institute for Advanced Study, 02455 Seoul, Korea Abstract An efficient molecular motor would deliver cargo to the target site at a high speed and in a punctual manner while consuming a minimal amount of energy. However, according to a recently formulated thermodynamic principle, known as the thermodynamic uncertainty relation, the travel distance of a motor and its variance are constrained by the free energy being consumed. Here we use the principle underlying the uncertainty relation to quantify the transport efficiency of molecular motors for varying ATP concentration ([ATP]) and applied load (f). Our analyses of experimental data find that transport efficiencies of the motors studied here are semi-optimized under the cellular condition. The efficiency is significantly deteriorated for a kinesin-1 mutant that has a longer neck-linker, which underscores the importance of molecular structure. It is remarkable to recognize that, among many possible directions for optimization, biological motors have evolved to optimize the transport efficiency in particular.

ST-71

Unusual morphogenesis and cell motility in an in vitro population of cancer cell lines Kyoung J. Lee Dept. of Physics, Korea University, Seoul, Korea Abstract Monoclonal cell line MDA-MB-231, which is highly malignant breast cancer cells, has been widely used and studied in laboratory in various cancer researches over many years. Yet, recently we observed two very intriguing new phenomena associated with them that have never been addressed previously. The first phenomenon is about the unexpected, beneficial, role of cellular senescence (a permanent cell-cycle arrest) transforming tissue structures around them serving as aggregating centers for non-senescent cells in their vicinity: The senescent MDA-MB-231 cells actively organize localized 3D cell-clusters in a confluent 2D tumor layer. We find that the biophysical mechanism underpinning the surprising phenomenon primarily involves mitotic cell-rounding , dynamic and differential cell attachments, and cellular chemotaxis. By incorporating these few biophysical factors, we were able to recapitulate the experimental observation via a cellular Potts Model. The second phenomenon is about the motilities of freely crawling MDA-MB-231 cells as well as those of the same cells within a confluent population on 2D system. In both cases, cells exhibit directional persistence in their crawling trajectories. Counterintuitively, we find that MDA-MB-231 cells in a dense population have a better motility compared with that of freely crawling cells. In other words, in a dense population the cells are not jammed to each other but their cell-to-cell contact interactions rather enhance the persistence in their movements. Again, via a cellular Potts model, we show that this property is closely related to the rotational motion of two adhering cells on 2D substrate. We explore the phase space of different motile behaviors spanned by two key parameters, namely, the strength of directional persistence of individual cells and its memory time, and discover an interesting crossover. 83

DDAP10, Huaqiao University, Xiamen China, 2018

ST-72

A landscape view on the interplay between EMT and cancer metastasis Chunhe Li Fudan University, [email protected] Abstract The epithelial-mesenchymal transition (EMT) is a basic developmental process that converts epithelial cells to mesenchymal cells. Although EMT might promote cancer metastasis, the molecular mechanisms for it remain to be fully clarified. To address this issue, we constructed an EMT-metastasis gene regulatory network model and quantified the potential landscape of cancer metastasis-promoting system computationally. We identified four steady state attractors on the landscape, which separately characterize anti-metastatic (A), metastatic (M), and two other intermediate (I1 and I2) cell states. The tetrastable landscape and the existence of intermediate states are consistent with recent single cell measurements. We identified one of the two intermediate states I1 as the EMT state. From a minimum action path (MAP) approach, we found that for metastatic progression, cells need to first undergo EMT (enter the I1 state), and then become metastatic (switch from the I1 state to the M state). Specifically, for metastatic progression, EMT genes (such as ZEB) should be activated before metastasis genes (such as BACH1). This suggests that temporal order is important for the activation of cellular programs in biological systems, and provides a possible mechanism of EMT promoting cancer metastasis. To identify possible therapeutic targets from this landscape view, we performed sensitivity analysis for individual molecular factors, and identified optimal interventions for landscape control. We found that minimizing transition actions more effectively identifies optimal combinations of targets that induce transitions between attractors than single factor sensitivity analysis. Overall, the landscape view not only suggests that intermediate states increase plasticity during cell fate decisions, providing a possible source for tumor heterogeneity that is critically important in metastatic progress, but also provides a way to identify therapeutic targets for preventing cancer progression.

ST-73

Time series analysis approach to protein dynamics Ming-Chya Wu1,2 1Research Center for Adaptive Data Analysis, National Central University, Zhongli, Taoyuan 32001 2 Institute of Physics, Academia Sinica, Nankang, Taipei 11529 Abstract Nanomechanics of molecular mechanoactivation is an interesting feature of muscle proteins. This report introduces implementations of time series analysis approaches to study dynamics of proteins involved in muscle tension generation. The results indicate the proposed analysis is promising in identifications of protein functional domains and in studies of functions of intrinsically disordered proteins.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-74

Spatial structure of bacterial behavior in chemotactic group migration Xiongfei Fu, [email protected] Shenzhen Institutes of Advanced Technology, CAS Abstract Chemotactic bacteria can migrate in a group in capillary tube by consuming chemoattractant collectively. Previous study showed this collective migration can coordinate individuals of different abilities by spatially sorting them within traveling group. Although this coordination mechanism was shown in theory, how individual cells resolve the conflict between heterogeneity and unified performance is still lack of direct evidence. Moreover, the model predicts the behavior feedback escalates the difference the of operational tumble bias between the front and the back. By analyzing single bacterial trajectories in the traveling wave, we show that not only self-organized, bacteria behavior is spatially structured: cells at the back of the wave run less tumble more, but they are better polarized to the direction of collective migration. The polarity advantage of cells in the back helps the higher tumble cells to exhibit a more efficient directional migration so as to keep up with the front cells.

ST-75

Dynamics behavior of the Microswimmer Yun-Yun Li (李云云) [email protected] Tongji University, Abstract How far away the fluctuations can drive the micro-engines and nano-engines beyond the limits of the macroscopic laws. This is the key question in modern thermodynamics. To answer this question, we discussed the design of high-efficiency autonomous machines for active (or self-propelling) Micro-swimmer. Meanwhile, based on Oseen approximation, the long-ranged hydrodynamic interaction is investigated when the active Micro-swimmers are trapped in a two-dimensional film. References: 1. T. Debnath, Y. Li*, P. K. Ghosh and F. Marchesoni, Hydrodynamic interaction of trapped active Janus particle in two dimensions, Phys. Rev. E 97. 042602 (2018) 2. Y. Li, F. Marchesoni, T. Debnath and P. K. Ghosh, Two-dimensional dynamics of a trapped active Brownian particle in a shear flow, Phys. Rev. E 96, 062138 (2017) 3. D. Debnath, P. Ghosh, Y. Li*, F. Marchesoni and B. Li, Diffusion of eccentric microswimmers, Soft Matter 12(7), 2017-2024(2016) 4. Y. Li, Pulak K Ghosh, F Marchesoni, Baowen Li, Manipulating chiral microswimmers in a channel, Phys. Rev. E 90, 062301, (2014)

DDAP10, Huaqiao University, Xiamen China, 2018

ST-76

Dynamical Sub-classes of Dry Active Nematics Li-bing Cai,1, 2 Hugues Chaté,3, 4, 2 Yu-qiang Ma,1, 2 and Xia-qing Shi2, 3 1.National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing 210093, China 2.Center for Soft Condensed Matter Physics and Interdisciplinary Research, Soochow University, Suzhou 215006, China 3.Service de Physique de l'Etat Condens_e, CEA, CNRS Universit_e Paris-Saclay, CEA-Saclay, 91191 Gif-sur-Yvette, France 4.Computational Science Research Center, Beijing 100094, China Abstract We show that the dominant mode of alignment plays an important role in dry active nematics, leading to two dynamical sub-classes defined by the nature of the instability of the nematic bands that characterize, in these systems, the coexistence phase separating the isotropic and fluctuating nematic states. In addition to the well-known instability inducing long undulations along the band, another stronger instability leading to the break-up of the band in many transversal segments may arise. We elucidate the origin of this strong instability for a realistic model of self-propelled rod and determine the high-order nonlinear terms responsible for it at the hydrodynamic level.[1] [1] Li-bing Cai, Hugues Chaté, Yu-qiang Ma and Xia-qing Shi, Dynamical sub-classes of dry active nematics, arXiv: 1807.05725

ST-77

Front instability in microbial range expansion with chemotaxis Leihan Tang, Beijing Computational Science Research Center Abstract Range expansion is a well-known phenomenon in the development of microbial communities. The speed of the expansion is the result of the often complex dynamics within the community as well as the presence of attractants in the surrounding medium. By coupling bacterial motility to the local cell density using a genetic circuit, Liu et al. observed ring patterns of varying wavelength in grown E. coli colonies [1]. The patterning process is caused by an instability when the motility drop in the frontal region becomes sufficiently sharp [2]. Here, we extend previous theoretical work by including the effect of a chemo-attractant in the medium. Consequently, a chemotactic group of cells appear in the frontal region which not only sets the front speed, but also affects the cell density profile. Depending on the ratio of cell densities in the chemotactic ring and in the motility transition region, different types of dynamic behavior can be observed. Through extensive numerical studies, we map out the phase diagram of the system and identify stability conditions of a steadily propagating front. Various characteristics in the pattern-forming region, such the front speed and oscillating frequency, will be discussed. Work supported in part by the RGC of the HKSAR under grant HKBU 12324716. [1] Chenli Liu et al., Science 334, 238 (2011). [2] Xiongfei Fu et al., Phys. Rev. Lett. 108, 198102 (2012). 86

DDAP10, Huaqiao University, Xiamen China, 2018

ST-78

Collective behaviors in active matter systems Xinliang Xu Beijing Computational Science Research Center Abstract Active matter systems are driven out of equilibrium as each individual propels itself, e.g. each fish swims in a school. These systems are of significant interest for their demonstration of fascinating collective behaviors that emerge. Here two studies about these emergent behaviors are reported. One is about the relation between two of these emergent behaviors (clustering and flocking, where our results show that flocking can suppress the formation of clusters in an active granular system. The other is about the response of the flocking behavior of bacterial suspensions, to the external drive of a shear flow, where the physics underlying the observed shear banding structure is unveiled.

ST-79

A reaction-diffusion model to study the electrophysiological integration of human pluripotent stem cell-derived cardiomyocytes into the myocardial infarction zone Xiang Gao Shaanxi Normal University, Xi’an, China Abstract Implanting human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) is a promising approach in cardiac regeneration therapy. However, in the early days of posttransplantation, the graft causes serious side effects, i.e. arrhythmia. Because arrhythmia can be life threatening, it needs to be understood and managed before safe clinical translation. We use a reaction-diffusion model to investigate the electrophysiological integration of hPSC-CMs into the host myocardial infarction zone. Especially, two typical cases in many experimental trials are idealized and considered in our simulations. The one is the injection of hPSC-CMs into transmural myocardial infarction which involves a three-dimensional zone through three-layers (epicardial, endocardial, and midmyocardial cells) heart wall. The other is the epicardial implantation of a cardiac patch which is a two-dimensional hPSC-CM sheet. We consider the electrophysiological properties and their differences between the engraft and host cardiomyocytes, such as cell morphology, gap junction, and action potential. Then, we study the mechanism of arrhythmia initiation caused by these electrophysiological inhomogeneities, and identify the basins of arrhythmia through a wide range of physiologically relevant conditions. The size of these basins can be used to determine the main causes of arrhythmia initiation from numerous possibilities. The model we use can quantitatively fit to the experimental results of action potential morphology and rate dependence of action potential duration on conduction velocity for both hPSC- and host cardiomyocytes. Yet former experiments were carried in different designs and preparations. The experimental data consistent for both grafted hPSC- and host cardiomyocytes is missing. Thus, our numerical results are semiquantitative, but still instructive enough to find the main causes of 87

DDAP10, Huaqiao University, Xiamen China, 2018 arrhythmia initiation among the numerous possibilities caused by the engraftment of hPSC-CMs. It will be helpful to design subsequent experiments which are efficient to verify our predictions. Key Words: Reaction-Diffusion System, Cardiac Arrhythmia, Cardiac Regeneration Therapy

ST-80

Dynamics of the early afterdepolarization in the ventricular myocyte Zhikun Chu and Xiaodong Huang Department of physics, South China University of Technology, Guangzhou, China

Abstract Early afterdepolarization (EAD) is a pathologic behavior in the ventricular myocyte, which is closely associated to the long QT sydrome and the relevant arrhythmogenesis [1, 2]. In physical studies it is revealed that such a kind of dynamical behavior is due to the Hopf bifurcation [3]. In order to study the details of the EAD dynamics, we simplify the phase I Luo-Rudy (LR1) ventricular myocyte model (a complex Hodgekin-Huxley type physiological model obtained from the guinea pig ventricular cells [4]) by polynomial approximation. By this way we can make the system analytically solvable. Based on the model we can discuss the following problems. i) The parameter conditions for the occurence of EAD. It was found that low potassium or high calcium current through the membrane may induce EAD, while sometimes over reduced potassium or increased calcium could eleminate it [5]. Also the varying time constant of the activation gate of the potassium current is deterministic for EAD [3]. By the solvable model we can analytically obtain the parameter conditions and explain the phenomena. ii) The amplitude and period of the EAD oscillations. The numerical results of the dependence of the EAD properties on the parameters will be shown. The results need theoretical explanation. We will propose our idea in the coming work on getting the amplitude equation of the EAD dynamics based on the the polynomial approximated LR1 model. As we know in the vicinity of the Hopf bifurcation point there exists a general equation (the Ginzberg-Landau equation) describing the evolvements of the limit cycle amplitude and the oscillatory period, which is dependent of the parameter set of the system. The amplitude equation may reveal and predict the characteristics of the EAD properties. Studying the complex biological dynamics via the way of physics has provided deeper insights into the subject. We hope this work may be helpful for the treatments of the relevant arrhythmogenesis. Key words: early afterdepolarization, nonlinear dynamics, amplitude equation References: [1] Z. Qu, G. Hu, A. Garfinkel and J. Weiss, Physics Reports, 543, 61 (2014). [2] X. Huang, T. Y. Kim, G. Koren, B.-R. Choi and Z. Qu, Am J Physiol Heart Circ Physiol, 311: H1470 (2016) [3] D. Tran, D. Sato, A. Yochelis, J. Weiss, A. Garfinkel, and Z. Qu, Phys. Rev. Lett., 102, 258103 (2009). [4] C. H. Luo and R. Rudy, Circ. Res., 68, 1501 (1991). [5] Z. Qu, L.-H. Xie, R. Olcese, H. Karagueuzian, P.-S. Chen, A. Garfinkel and J.Weiss, Cardiovascular Research, 99, 6 (2013).

DDAP10, Huaqiao University, Xiamen China, 2018

ST-81

Cardiac Dynamics in response of periodic external stretching Po-Yu Chou1, Wei-Yin Chiang1, C.K. Chan2 and Pik-Yin Lai1 1Dept. of Physics and Center for complex system, National Central University, Taoyuan City 320 2Institute of Physics, Academia Sinica, Nankang, Taipei 115 Abstract Effects of mechanical coupling to cardiac dynamics are studied by monitoring the beating dynamics of a cardiac tissue which is being pulled periodically. The tissue is taken from the heart of a bullfrog. Since the tissue includes pacemakers, it beats spontaneously; giving an almost constant interbeat interval (IBI). However, the IBI can be observed to vary under a periodic external drive. Interestingly, when the period of the external drive is about two times the intrinsic IBI of the tissue without pulling, the IBIs as a function of time appear in a form similar to a wave packet. Our experimental results can be understood theoretically by a phase-coupled model under external driving. In particular, the theoretical prediction of the wave-packet period as a function of the normalized driving period agrees excellently with the observation. Furthermore, the mechanical coupling constant can be extracted from the experimental data using our model and found to be insensitive to the external driving period. Key Words: Cardiac dynamics, phase model

ST-82

Jacobian-Determinant Method of Identifying Phase Singularity During Reentry Hong Zhang Zhejiang University, [email protected] Abstract Reentrant spiral waves (also called rotors) have been observed in the heart muscle during cardiac arrhythmias and are nowadays targeted by ablation therapy in order to cure certain heart rhythm disorders. Phase singularity (PS) is considered to represent the organizing center of the spiral wave and the spatiotemporal behavior of the spiral wave can be extensively quantified by tracking PS. Recent clinical studies suggested that ablating the tissue at PS locations may cure atrial fibrillation. However, for experimental data, typically only one state variable, i.e. the voltage, is recorded. Accordingly, the calculation of phase and the identification of PSs must been carried out using one variable. Here, a Jacobian-determinant method using one variable in reconstructed state space is proposed and the advantage of the method is demonstrated.

ST-83

Pattern recognition though spike-timing-dependent plasticity in spiking neural networks Xuhui Huang Institute of Automation, CAS Abstract 89

DDAP10, Huaqiao University, Xiamen China, 2018

In this talk, we will report that how spiking neural networks (SNNs) with spike-timing-dependent plasticity (STDP) could realize both visual and auditory classification tasks in a biologically plausible (bio-plausible) way. For the visual task, we are interested in how supervised learning (SL) could be achieved based on both biologically realistic learning rules and the spike-based computation, since there is a big challenge for training SNNs in a supervised way due to the non-differentiable spike activities. We proposed a novel bio-plausible SNN model for SL, based on a symmetric STDP rule, synaptic scaling as well as dynamical threshold. Our model can achieve high recognition performance in the benchmark MNIST task. For the auditory task, we had explored if and how an STDP-based convolutional SNN could solve the automatic speech recognition (ASR) efficiently. Currently, SNN models for ASR are either lack of high-efficiency learning algorithms for training, or trained in biologically implausible ways. We proposed a biologically inspired convolutional SNN model for ASR with the following four key characteristics: STDP for unsupervised learning, the time-to-first-spike coding scheme for fast and efficient information processing, the strategy of local weight sharing for better feature extraction. With the linear support vector machine (SVM) as a classifier, our model could achieve state-of-the-art classification performance on the TIDIGITS dataset comparing to the best results of artificial neural networks. Furthermore, the validity of our model can also be well maintain ed by using the spike-based classifier---tempotron as well as on a more difficult speech classification task (17 categories of words) based on the TIMIT dataset. Our studies indicate that SNNs equipped with STDP are capable of efficiently solving both visual and auditory recognition tasks in a bio-plausible way. Reference: 1. Yunzhe Hao*, Xuhui Huang*, Meng Dong, Bo Xu. A Biologically Plausible Supervised Learning Method by Dopamine-modulated STDP in Spiking Neural Networks. Submitted to AAAI 2019 (* Co-first authors) 2. Meng Dong, Xuhui Huang*, Bo Xu*. Unsupervised speech recognition through spike timing-dependent plasticity in a convolutional spiking neural network. PLOS ONE (Accepted, * Co-corresponding authors)

ST-84

Network modeling for neuromodulatory intervention Zigang Huang Xi’an Jiaotong University Abstract Neuromodulation can be understood from two aspects. It refers to the natural physiological process in the nervous system, and also the therapeutic electromagnetic or chemical stimulation of nerve cells. The two different subjects of neuromodulation have intimate correlation since the later one works based on the former one. In this talk, we give a brief introduction to the two concepts of neuromodulation, and then show how the techniques from complex system theory are adopted in the related studies.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-85

Synaptic Correlates of Working Memory Capacity Yuanyuan Mi Academy of Military Medical Sciences [email protected] Abstract Psychological studies indicate that human ability to keep information in readily accessible working memory is limited to 4 items for most of the people. This extremely low capacity severely limits execution of many cognitive tasks, but its neuronal underpinnings remain unclear. Here we show that in the framework of synaptic theory of working memory, capacity can be analytically estimated to scale with characteristic time of short-term synaptic depression relative to synaptic current time constant. The number of items in working memory can be regulated by external excitation, enabling the system to be tuned to the desired load and to clear the working memory of currently held items to make room for new ones.

ST-86

Resonance firing in single autaptic neuron Hengtong Wang Shaanxi Normal University Abstract Autapses, synapses between a neuron and itself, are believed to serve as a regulator of information processing in the nervous system. Noise, as random or unpredictable fluctuations, affects nearly all aspects of nervous function. Here, I report the resonance firing induced by autaptic connection with or without external inputs.

ST-87

Unified Analysis of Global and Focal Aspects of Absence Epilepsy via Neural Field Theory of Corticothalamic System [1]

Dong-Ping Yang1,2*, P. A. Robinson1,2

1 School of Physics, University of Sydney, NSW 2006, Australia 2 Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia * Corresponding author Email: [email protected] Abstract Absence epilepsy is characterized by a sudden paroxysmal loss of consciousness accompanied by oscillatory activity propagating over many brain areas. Although primary generalized absence seizures are supported by the global corticothalamic system [2-6], converging experimental evidence supports a focal theory of absence epilepsy [7-9]. Here we propose a dynamical mechanism to unify the global and 91

DDAP10, Huaqiao University, Xiamen China, 2018 focal aspects of absence epilepsy, with focal absence seizures associated with seizure localization, and the global ones associated with seizure generalization. A physiology-based corticothalamic model [4-6] is investigated with focal spatial heterogeneity, to unify global and focal aspects of absence epilepsy. Numerical and analytical calculations are employed to investigate the emergent spatiotemporal dynamics induced by focal activity as well as their underlying dynamical mechanisms. The spatiotemporal dynamics can be categorized into three scenarios: suppression, localization, and generalization of the focal activity, as summarized from a phase diagram vs. focal width and characteristic axon range. The corresponding temporal frequencies and spatial extents of wave activity during seizure generalization and localization agree well with experimental observations of global and focal aspects of absence epilepsy, respectively. Our corticothalamic model is used to investigate how seizure rhythms and spatial extents are related in these two different aspects of absence epilepsy. The emergent seizure localization provides a biophysical explanation of the temporally higher frequency, but spatially more localized cortical waves observed in genetic rat models that display characteristics of human absence epilepsy. The results account for the difference of the experimentally observed seizure rhythms and spatial extents between humans and genetic rat models, which has previously been used to argue against the validity of such rats as animal models of absence epilepsy in humans [3, 8, 9]. Predictions are also presented for further experimental test. Reference: [1] On-line preprint: doi: https://doi.org/10.1101/339366 [2] Destexhe A. Spike-and-Wave Oscillations Based on the Properties of GABAB Receptors. J Neurosci. 1998; 18(21): 9099–9111. [3] Destexhe A. Can GABAA conductances explain the fast oscillation frequency of absence seizures in rodents? Eur J Neurosci. 1999; 11(6): 2175–2181. [4] Robinson PA, Rennie CJ, Rowe DL. Dynamics of large-scale brain activity in normal arousal states and epileptic seizures. Phys Rev E. 2002; 65: 041924. [5] Breakspear M, Roberts JA, Terry JR, Rodrigues S, Mahant N, Robinson PA. A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis. Cereb Cortex. 2006; 16(9): 1296–1313. [6] Yang DP, Robinson PA. Critical dynamics of Hopf bifurcations in the corticothalamic system: Transitions from normal arousal states to epileptic seizures. Phys Rev E. 2017; 95: 042410. [7] Meeren HKM, Pijn JPM, van Luijtelaar ELJM, Coenen AML, Lopes da Silva FH. Cortical Focus Drives Widespread Corticothalamic Networks during Spontaneous Absence Seizures in Rats. J Neurosci. 2002; 22(4): 1480–1495. [8] Meeren HKM, van Luijtelaar G, Lopes da Silva FH, Coenen AML. Evolving concepts on the pathophysiology of absence seizures: The cortical focus theory. JAMA Neurol. 2005; 62(3): 371–376. [9] Depaulis A, Charpier S. Pathophysiology of absence epilepsy: Insights from genetic models. Neurosci Lett. 2018; 667: 53 – 65.

ST-88

The impacts of cross-shareholding network on extreme price movements: Evidence from China Jie Cao1, FenghuaWen1,2,3 1 School of Business, Central South University, Changsha, Hunan 410083, China, 92

DDAP10, Huaqiao University, Xiamen China, 2018

2 Supply Chain and Logistics Optimization Research Centre, Faculty of Engineering, University of Windsor, Windsor, ON, Canada 3 Centre for Computational Finance and Economic Agents, University of Essex, Colchester CO4 3SQ, UK Abstract Extreme price movement of listed company is highly concerned by the regulatory authorities and investors, which is influenced by its information transparency. Crossshareholding relationship links the companies into a complex network structure, which not only provides channel for information flow, collection and diffusion, but also has important impacts on information transparency and extreme price movement. Due to the linkages and position, more central nodes in the network are thought to be able to collection and dissemination information more efficiently. On the basis, this paper explores the impacts of cross-shareholding network centrality on extreme price movements by using ownership information of Chinese listed companies from 2004 to 2016. The cross-shareholding network centralities associated with each firm are measured by several widely used centralities of nodes, i.e., the eigenvector centrality, the degree centrality and the betweenness centrality, because they can not only measure the connections that nodes linked but also assess the ability to collect information. By means of panel regression technology, we analyze the impacts of network structure on extreme price movements. Furthermore, considering the different attitudes to positive news and negative news of managers, the transparency of these two kinds of news may be different. In order to test the different effects of network centrality on extreme price upward movement and downward movement, we also test asymmetric impacts of network centrality on both price limit up and price limit down, respectively. Our empirical findings can be summarized as follows. To begin with, the fat-tailed phenomenon of cross-shareholding networks centrality of Chinese stock market is observed, i.e., a few listed companies are at a high-level in centrality, they have more advantages in information publishing and collecting, while other companies with much less centrality are just the opposite. Secondly, we estimate the relationships between the network structures and extreme price movements by regressing the network centralities and the total numbers of price limit of each company. The results show that network centralities have negative relationship with extreme price movements, suggesting that the centrality of company improve the information transparency, firms positioned more centrally therefore have less extreme price movements. Furthermore, we test the asymmetry effects of network centrality on both extreme price upward movements and extreme price downward movement. Regression results show that network centralities have a stronger negative relationship with extreme price upward movement than downward movements because good news about stock price is more likely to be published by managers. Keywords: Extreme price movements; cross-shareholding networks; centrality; information transparency.

ST-89

European stock market dynamics based on correlation analysis Tian Qiu, Nanchang Hangkong University Abstract The European stock market dynamics has been investigated, focusing on understanding the European union effect on European stock index dynamics. More than ten European countries are 93

DDAP10, Huaqiao University, Xiamen China, 2018 studied, which are distinguished as the formerly joining the European Union (EU) countries (F-EU), later joining the European Union countries (L-EU), and not joining the European Union countries (N-EU). By investigating the correlations and partial correlations between different countries, the union is found to have a great impact on the stock markets of F-EU countries, but show little influence on those of the L-EU countries. Moreover, the American stock markets and the stock markets of the F-EU countries present a high influence on the European stock index in the global stock markets. The correlations of the stock markets of L-EU and N-EU countries are found to increase in the financial crisis period.

ST-90

Evolution of volatility after extreme events Fei Ren, East China University of Science and Technology Abstract This paper studies the evolution of volatility after extreme events, including endogenous and exogenous events. We decompose the volatility of index into three components, i.e., the market mode, sector mode and random mode. Different from the most disordered systems which follow either stretched exponential or power law decays, we find that the volatility after extreme events for the index and its market mode follows the stretched exponential decay at the first stage and then the power law decay at the second stage. Furthermore, we propose a new dynamic model based on the Wiener process, which can successfully explain the empirical results of volatility evolution after endogenous and exogenous events for both index and its components. We use the hypotheses of sequential arrival of information hypothesis (SAIH) and mixture of distribution hypothesis (MDH) to further examine our dynamic model, which comprise of two stages that investors transform from uninformed to informed at the first stage and then informed investors dominate the market at the second stage.

ST-91

Momentum and contrarian effects in the Chinese stock market: From perspective of the Adaptive Markets Hypothesis Huailong Shi, Nanjing University of Information Science & Technology Abstract This work firstly investigates two versions of momentum and contrarian effects. Using the data of A-share stocks listed on the Shanghai Stock Exchange and Shenzhen Stock Exchange, we examine the profitability of cross-sectional contrarian portfolios and find evidence of short-term and long-term contrarian profitability in the whole sample period when the look-back and holding horizons are 1 month or longer than 12 months and the annualized returns of contrarian portfolios increases with the look-back and holding horizons. Generally, loser portfolios and winner portfolios have positive returns and loser portfolios perform much better than winner portfolios. These results still hold after robustness checks including skipping one-month between the look-back and holding periods, and performing subperiod analysis. We also study the time series momentum or contrarian effects in the Chinese stock market. We evaluate the performance of the time series momentum strategy applied to major stock 94

DDAP10, Huaqiao University, Xiamen China, 2018 indices in mainland China and explore the relation between the performance of time series momentum strategies and some firm-specific characteristics. Our findings indicate that there is a time series momentum effect in the short run and a contrarian effect in the long run in the Chinese stock market. The performances of the time series momentum and contrarian strategies are highly dependent on the look-back and holding periods and firm-specific characteristics. At last, in light of Adaptive Markets Hypothesis (AMH), the study is conducted to investigate the evolving arbitrage opportunities within the framework of momentum and contrarian effects. The results reveal that the arbitrage opportunities based on the contrarian portfolios wax and wane over time. The performance of contrarian portfolios is highly dependent on several market conditions. The periods with upward trend of market state, higher market volatility and liquidity, lower macroeconomics uncertainty are related to higher contrarian profitability. These findings are consistent with the Adaptive Markets Hypothesis and have practical implications for market participants.

ST-92

Multifractal analysis of financial markets: A review Weixing Zhou, East China University of Science and Technology Abstract Multifractality is ubiquitously observed in complex natural and socioeconomic systems. Multifractal analysis provides powerful tools to understand the complex nonlinear nature of time series in diverse fields. Inspired by its striking analogy with hydrodynamic turbulence, from which the idea of multifractality originated, multifractal analysis of financial markets has bloomed, forming one of the main directions of econophysics. We review the multifractal analysis methods and multifractal models adopted in or invented for financial time series and their subtle properties, which are applicable to time series in other disciplines. We survey the cumulating evidence for the presence of multifractality in financial time series in different markets and at different time periods and discuss the sources of multifractality. The usefulness of multifractal analysis in quantifying market inefficiency, in supporting risk management and in developing other applications is presented. We finally discuss open problems and further directions of multifractal analysis.

ST-93

Trading Volume and Return Volatility of Bitcoin Market: Evidence for the Sequential Information Arrival Hypothesis Dehua Shen, Tianjin University Abstract This paper gives the first empirical evidence on the relationships between trading volume and return volatility of the Bitcoin market for 15 foreign currencies by investigating two competing hypotheses, i.e., Mixture of Distribution Hypothesis (MDH) and Sequential Information Arrival Hypothesis (SIAH). Allowing for both linear and nonlinear correlation and causality tests, the empirical results mainly show that: first, trading volume and return volatility are negatively correlated, implying a lack of support for the MDH; second, we document significant lead–lag relationships between trading volume and return volatility, which support the SIAH; third, the results are robust to alternative measurements 95

DDAP10, Huaqiao University, Xiamen China, 2018 of trading volume and sub-period analysis. Generally speaking, these findings have practical implications for investors, who are interested in investing in Bitcoin market.

ST-94

Universal statistics of housekeeping entropy Jae Dong Noh, Department of Physics, University of Seoul, Seoul 02504, Korea Abstract The entropy production of a nonequilibrium system with broken detailed balance is a random variable whose mean value is nonnegative. Among the total entropy production, the housekeeping entropy production is associated with the heat dissipation in maintaining a nonequilibrium steady state. We derive a Langevin-type stochastic differential equation for the housekeeping entropy production. The equation allows us to define a housekeeping entropic time $\tau$. Remarkably, it turns out that the probability distribution of the housekeeping entropy production at a fixed value of $\tau$ is given by the Gaussian distribution regardless of system details. The Gaussian distribution is universal for any systems, whether in the steady state or in the transient state, whether they are driven by time-independent or time-dependent driving forces. We demonstrate the universal distribution numerically for model systems.

ST-95

Nonequilibrium Topological Insulator with Ultracold Atoms Ying Hu State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Laser Spectroscopy, Shanxi University, Taiyuan, Shanxi 030006, China Abstract I will present some recent work on topological insulators far from equilibrium. The first part is relevant to ongoing experiments aimed at realizing synthetic quantum Hall materials with ultracold atoms with dynamical means. At zero temperature, the direct correspondence between the Chern number of the ground state and Hall conductance is well established. Yet, non-equilibrium scenarios generically occur in atomic setups, where starting from a topologically trivial initial state, the Hamiltonian is driven into a topological regime. However, the system can never enter a topological insulator state under coherent dynamics, raising the challenge as to which manifestations of topology can be actually observed. Addressing this issue, we show that a quantized nonequilibrium Hall response can build up without a topological state [1]. I will proceed to discuss how periodic driving can be harnessed to transcend fundamental constraints in equilibrium systems. Specifically, I show that the phenomenon of perfect spin momentum locking, which in equilibrium can only occur in 2D static topological insulators, emerge in the stroboscopic dynamics of 1D lattice systems [2]. [1]. Ying Hu, Peter Zoller, Jan Carl Budich, Dynamical Buildup of a Quantized Hall Response from Non-Topological States, Phys. Rev. Lett. 117, 126803 (2016) [2]. Jan Carl Budich, Ying Hu, Peter Zoller, Helical Floquet Channels in 1D Lattices, Phys. Rev. Lett. 118, 105302 (2017) 96

DDAP10, Huaqiao University, Xiamen China, 2018

Key words: Non-Equilibrium, Coherent dynamics, Topology

ST-96

Sensitive Dependence of Optimal Network Dynamics on Network Structure Jie Sun Department of Mathematics, 8 Clarkson Ave, Potsdam, NY 13699-5815, USA Abstract The relation between network structure and dynamics is determinant for the behavior of complex systems in numerous domains. An important long-standing problem concerns the properties of the networks that optimize the dynamics with respect to a given performance measure. Here, we show that such optimization can lead to sensitive dependence of the dynamics on the structure of the network. Specifically, using diffusively coupled systems as examples, we demonstrate that the stability of a dynamical state can exhibit sensitivity to unweighted structural perturbations (i.e., link removals and node additions) for undirected optimal networks and to weighted perturbations (i.e., small changes in link weights) for directed optimal networks. As mechanisms underlying this sensitivity, we identify discontinuous transitions occurring in the complement of undirected optimal networks and the prevalence of eigenvector degeneracy in directed optimal networks. These findings establish a unified characterization of networks optimized for dynamical stability, which we illustrate using Turing instability in activator-inhibitor systems, synchronization in power-grid networks, network diffusion, and several other network processes. Our results suggest that the network structure of a complex system operating near an optimum can potentially be fine-tuned for a significantly enhanced stability compared to what one might expect from simple extrapolation. On the other hand, they also suggest constraints on how close to the optimum the system can be in practice. Finally, the results have potential implications for biophysical networks, which have evolved under the competing pressures of optimizing fitness while remaining robust against perturbations.

ST-97

Bi-motility model to Bacterial Stripe Formation Fei Qi, Hong Kong Baptist University Abstract Nature has built a variety of fascinating temporal and spatial patterns that are vital to the movement, growth, and colonies of living organisms. A well-known example is the periodic stripe formation in Escherichia coli colonies growth and self-organization. [Science 334, 238, 2011]. In order to decipher the mechanisms of this pattern formation process, we introduced a Bi-motility model which couples cell density and motility to modulate them propagating colony front. The numerical simulations of the model show a phase transition between steadily-propagating and oscillatory fronts. Moreover, In nutrient surroundings, the consequence of motility genes are activated and led to the wave-front propagation could transition from pull to pushed wave via chemotaxis effect along self-generated chemotactic gradients. The influence of various parameters for those transitions is discussed.

DDAP10, Huaqiao University, Xiamen China, 2018

ST-98

Dynamics and Periodic solutions of two state dependent delay network congestion control models Lijun Pei ,Zhengzhou University Abstract Network congestion is an important problem that can cause data packet loss and communication delay. It restricts the overall performance of the network. In order to avoid the appearance of network congestion, network control becomes the first problem to be solved. In this paper, two wireless network congestion control models are considered. One is the traditional TCP control protocol. It is a model with the state variables of the size of the window, the average queue length and the state-dependent round trip time-delay, in which the round trip time depends on the transmission delay and the average queue length. It mainly avoids congestion by controlling network efficiency and packet loss probability. The other is an improved FAST TCP protocol. It takes congestion windows and queue delays as state variables, and its round-trip time is state-dependent variable, in which round-trip delay depends on queue delay. The queue delay is a congestion control variable, which is more accurate than the loss probability estimation. The problem of the approximate solution of periodic solution for dynamic system can be solved by shooting method, finite difference method, etc. However, these methods are time-consuming and have strong dependence on the initial value, so we use the semi-analytical and semi-numerical method, i.e., harmonic balance frequency time conversion (harmonic balance and altering frequency domain, HB-AFT) method to solve this problem. The AFT technique in this method can avoid the integration or series expansion of complex nonlinear term, and this method can obtain the high order harmonic solution of good precision for the system's steady state response. The main process is as follows. First, harmonic solution and harmonic balance process are first carried out. Then, according to the harmonic balance, we can obtain algebraic equations that confirm the relationship of all ordered harmonic coefficients. The unknown variables are given by the AFT process. Last, we use Newton iteration to solve their fixed point, that is, the coefficients matrix of periodic solution. It is proved by Winpp’s numerical simulation that this method is effective and correct for state-dependent delayed dynamic systems. And their other dynamics are also considered in this paper. Its practical significance is that we can avoid periodic oscillations by adjusting parameters to improve the network performance. Key words: nonlinear system; state-dependent delayed system; congestion model of wireless network; periodic solution.

ST-99

Quantum thermal transport in nonequilibrium spin-boson systems Chen Wang Department of Physics, Zhejiang Normal University, Jinhua 321004, Zhejiang P. R. China Abstract Nonequilibrium thermal transport in small quantum systems has been attracting much attention in nonequilibrium statistics. As a representative system, the spin-boson model has been extensively investigated in quantum transport, quantum thermodynamics and quantum optics. In particular, the 98

DDAP10, Huaqiao University, Xiamen China, 2018 system-bath interaction is unraveled to the crucial to modulate heat current and fluctuations[1,2]. In this report, we mainly discuss how strong spin-bath coupling affect the mechanism of nonequilibrium transfer processes. The main content includes two parts: i) we propose a nonequilibrium polaron-transformed Redfield equation (NE-PTRE) method, and combine it with the full counting statistics to unify heat current and fluctuations in the nonequilibrium spin-boson model from an analytical perspective[3]. Moreover, we extend the NE-PTRE to analyze the geometric-phase-induce heat transport[4]. ii) based on three-terminal setup, we investigate the negative differential thermal conductance and thermal transistor behaviors[5]. Keywords: quantum thermal transport, geometric-phase-induced heat current, thermal transistor, nonequilibrium spin-boson model [1] J. Ren, P. Hanggi and B. W. Li, Berry-Phase-Induced Heat Pumping and its Impact on the Fluctuation Theorem, Phys. Rev. Lett. 2010, 104, 170601. [2] T. Chen, X. B. Wang and J. Ren, Dynamic control of quantum geometric heat flux in a nonequilibrium spin-boson model, Phys. Rev. B 2013, 87, 144303. [3] C. Wang, J. Ren and J. S. Cao, Nonequilibrium energy transfer at nanoscale: A unified theory from weak to strong coupling, Sci. Rep. 2015, 5, 11787. [4] C. Wang, J. Ren and J. S. Cao, Unifying quantum heat transfer in a nonequilibrium spin-boson model with full counting statistics, Phys. Rev. A 2017, 95, 023601. [5] C. Wang, X. M. Chen, K. W. Sun and J. Ren, Heat amplification and negative differential thermal conductance in a strongly coupled nonequilibrium spin-boson system, Phys. Rev. A 2018, 97, 052112.

DDAP10, Huaqiao University, Xiamen China, 2018

5.2.3 Posters

PS-01

The non-equilibrium universality class of forget-remember mechanism on the (1+1)-d lattice Jianmin Shen Huazhong Normal University Email:[email protected] It is very compelling to study the universe behaviours of non-equilibrium phase transitions, which may lead to the deep understandings of their interconnections and controls of these complex systems. This work was based on the Monte Carlo(MC) simulations for the study of the universality class of forget--remember mechanism(FRM) model of epidemic spreading on lattice, aiming to prove that the forget--remember mechanism(FRM) model of epidemic spreading belongs to the directed percolation (DP) universality class.

PS-02

The characteristics of Loop occupation base on networks classification by machine learning Wenjun Zhang Huazhong Normal University Email:[email protected] In study the different kinds of networks process, we notice that some special properties are different between them. But there are no exact expression to describe it. So we define it as Loop Vertices Occupation (LVO), and given it’s algorithm. Then we study it in basic networks’ model and real networks. The results in networks’ model is intuitional. The LVO is keep constant in different size Erdős–Rényi networks with same average degree, and increase with average degree meanwhile appear critical phenomena. We compare the LVO between networks’ model and some real networks and found the value of LVO are smaller than networks’ model in general. The three kinds of transport networks was studied by calculated LVO. There are have different results in different kind of transport networks. This is very useful in machine learning characteristic choose in network recognition.

PS-03

Dynamic Simulation and Structural Optimization of Pulse Separation Device in Dual Pulse Solid Rocket Motor LE Hao1,QIN Hai Yi2,YANG Ming1,LOU Yong Chun2,CUI Zhen Shan1 （1.Shanghai Space Propulsion Technology Research Institute,2. Shanghai Jiao Tong University, No. 1777 Zhong Chun Road, Minhang District, Shanghai ,201109,[email protected]） Abstract Recently, in order to improve the spacecraft’s control performance, viability and propulsive range,

100

DDAP10, Huaqiao University, Xiamen China, 2018 dual pulse solid rocket motor has gradually become the main power device of spacecraft. Usually, the dual pulse motor consists of two burning chambers, separated by a bulkhead, designated as pulse separation device (PSD) which achieves the two start-up operation and provides intermittent thrust. Pulse separation device, the most important component of motor, is composed of spoke structure and metal diaphragm with prefabricated notches. Actually, the PSD protects the propellant grain in the second pulse chamber against high temperature and pressure impact during the first pulse operation. At initiation of the second pulse, the PSD reliably opens due to the gas pressure. It opens along predetermined weakened lines without shedding any parts which could choke or damage the nozzle. Then, the combustion gas from the second pulse chamber can pass the empty first pulse chamber and flow through the nozzle to atmosphere. Compared with traditional one, dual pulse motor has more complex structure such as PSD, which makes it difficult to study its dynamic feature. Therefore, a method was set up here to solve this problem. Firstly, the Hyper-mesh was used to model the PSD. Then, based on LS-DYNA, the instantaneous opening process of PSD, the formation and movement of debris, and the impact damage to motor were simulated completely. Meanwhile, the dynamic characteristics of the PSD were obtained by cold-flow test and firing test. With the comparison, it was know that there was a high correlation between simulation and test, which proved that the model established was correct. Consequently, it was regarded as the standard method for analysis. By using the above method, the influence of different structural attributes such as thickness, groove depth, quantity, length on the opening pressure as well as failure mode of PSD were obtained. Later, based on the results, the optimal design of the PSD was determined. Moreover, the influence of gas temperature and ignition shock on the opening process of the PSD were analyzed. The results showed that both of them had little effect on the opening property of PSD in operative conditions. Finally, based on the simulation analysis and experimental study, the method suitable for analyzing the whole opening process of the PSD was systematically proposed, which has great significance on the structural optimization of pulse separation device as well as dynamic characteristics of solid rocket motor. Keyword: dual pulse solid rocket motor; pulse separation device; dynamic simulation; structural optimization

PS-04

Related research on failure phenomenon in scientific research Yanmeng Xing School of Systems Science, Beijing Normal University, Beijing, 100875, PR China Email:[email protected]

Abstract The development of science depends on the unremitting efforts of scientific researchers. Attracting more new scientific researchers to continue to engage in scientific work is the driving force of scientific development. A newcomer to scientific research - loosely defined as the author who at least publishes one journal paper in the first year of his or her career. We collected large-scale career histories of physicists in the APS data set and analyzed the failure of the newcomers in scientific research .We find that the length of research career shows a power law distribution, which indicates that a large number of newcomers stop doing research after publishing a paper. Individual science career can be divided 101

DDAP10, Huaqiao University, Xiamen China, 2018 into four stages: first career, fluctuating period, stable period and retirement.38% newcomers quit scientific research after their first year. A relative 25% of researchers quit during the second period. The relative failure rate of stable period per year is less than 5%.After 35 years of scientific careers, scientists are retiring from science. Next, we explore the factors related to failure at different stages. We find that the individual performance (the amount of published articles and the quality of articles) and the team level at the first 2 stages respectively play a different role. There is no obvious correlation between the latter two stages and the earlier stage. Failure may be caused by random factors or other complex mechanisms. To sum up, our findings provide clear guidance for newcomers on how to continue their scientific work, under where they should join a high-level team, improve the quality and quantity of articles and maintain the initiative, which can fundamentally promote the development of contemporary science to a certain extent. Keywords: New scientific researchers, Failure, Science career

PS-05

Chemically propelled motors navigate chemical patterns Jiangxing Chen Hangzhou University of Electronic Science and Technology Email:[email protected] Abstract Very small synthetic motors that use chemical reactions to drive their motion are being studied widely because of their potential applications, which often involve active transport and dynamics on nanoscales. Like biological molecular machines, they must be able to perform their tasks in complex, highly fluctuating environments that can form chemical patterns with diverse structures. Motors in such systems can actively assemble into dynamic clusters and other unique nonequilibrium states. It is shown how chemical patterns with small characteristic dimensions may be utilized to suppress rotational Brownian motions of motors and guide them to move along prescribed paths, properties that can be exploited in applications. In systems with larger pattern length scales, domains can serve as catch basins for motors through chemotactic effects. The resulting collective motor dynamics in such confining domains can be used to explore new aspects of active particle collective dynamics or promote specific types of active self-assembly. More generally, when chemically self-propelled motors operate in far-from equilibrium active chemical media the variety of possible phenomena and the scope of their potential applications are substantially increased.

PS-06

A New Method of Identifying Vortex Filaments During Ventricular Arrhythmia Wen-Jing Zheng School of Physics and Information Technology, Shaanxi Normal University, Xi'an 710062, China Abstract ： Pattern is a non-uniform macroscopic structure with some regularity in space or time, which is ubiquitous in nature. Scroll waves are the most common type of non-equilibrium patterns and can be widely observed in oscillating systems, excitation systems, and bistable systems. The self-organized dynamics of the scroll wave is the basis for the formation of complex spatio-temporal patterns in many 102

DDAP10, Huaqiao University, Xiamen China, 2018 excitable states of chemical and biological systems. Scroll waves can cause ventricular tachycardia, and their breakup to spatiotemporal chaos will lead to the transition of the heart from tachycardia to ventricular fibrillation. The filament is the main component of the scroll wave. It has special properties and governs the motion of the scroll wave. Therefore, the dynamic behavior of the scroll wave can be reflected by the filament, and is fundamental to effectively control tachycardia. In the cardiac ventricle, a three-dimensional coiled filament is the self-organizing center of scroll wave, or called vortex, which is corresponding to life-threatening arrhythmia. In the event of arrhythmia, scroll waves are usually observed. In order to eliminate the scroll wave, it is often necessary to determine the location and drift of scroll wave filament. Researches on the kinetics, instability and control of scroll wave are of guiding significance for the prevention and treatment of heart disease such as ventricular Arrhythmia. In this paper, we use the reaction diffusion system. In the numerical simulation, we found that in the three-dimensional gradient system, the gradient size and the direction of that can effectively determine the spiral wave in the adjacent time. So we propose a Jacobian-determinant method to identify the filament location and its spatio-temporal variation just using action potential measured in experiments. Using this method, we can describe the instantaneous drift of the scroll wave filament.At the same time, there is no need to calculate the local phase, which also facilitates a simpler determination of the spiral wave filament. We conclude that the use of Jacobian-determinant method allows for a novel means to identify the localization of filament.

Fig. 1 (color online)In Luo-Rudy 1 model an instant time t=490s and a time delay τ=10s are chosen. A snapshot of the spatial distribution of the fast variable V. The red line is the filament of current wave. Key Words scroll wave, vortex filaments,, arrhythmia

PS-07

Ergodicity breaking and ageing of underdamped Brownian dynamics with quenched disorder Wei Guo1, Lu-Chun Du2 1 Department of Physics, Kunming University, Kunming 650214, China 2 Department of Physics, Yunnan University, Kunming 650091, China E-mail: [email protected] (W. Guo) and [email protected] (L. C. Du) The dynamics of an underdamped Brownian particle moving in one-dimensional quenched disorder under the action of an external force is investigated. Within the tailored parameter regime, the 103

DDAP10, Huaqiao University, Xiamen China, 2018 transiently anomalous diffusion and ergodicity breaking, spanning several orders of magnitude in time, have been obtained. The ageing nature of the system weakens as the dissipation of the system increases for other given parameters. Its origin is ascribed to the highly local heterogeneity of the disorder. Two kinds of approximations (in the stationary state), respectively, for large bias and large damping are derived. These results may be helpful in further understanding the nontrivial response of nonlinear dynamics, and also have potential applications to experiments and activities of biological processes. Keywords: Ergodicity Breaking; Ageing; Transiently anomalous diffusion Acknowledgments: We are grateful to Professor Guo-Jun Jin for stimulating and illuminating discussions. This work was supported by the National Natural Science Foundation of China (Grant No. 11547027 and No. 11505149). Reference: Wei Guo*, Yong Li, Wen-Hua Song, and Lu-Chun Du*, Ergodicity breaking and ageing of underdamped Brownian dynamics with quenched disorder, J. Stat. Mech. 2018, 033303 (2018). Wei Guo*, Lu-Chun Du, Zhen-Zhen Liu, Hai Yang, and Dong-Cheng Mei, Uphill anomalous transport in a deterministic system with speed-dependent friction coefficient, Chin. Phys. B 26, 010502 (2017) .

PS-08

Nonlinear Dynamic Response of High-aspect-ratio Wings Chen Zhanjun, Liu Ziqiang China Academy of Aerospace Aerodynamics, 100074, Beijing, China Email: [email protected] Abstract: As one of main green energy, solar energy is attracted more attentions. According to the improvements of solar panel efficiency, composite material, many vehicles powered by solar energy are raised. In all the vehicles that impulse by solar energy, High Altitude, Long Endurance (HALE) aircraft gets more attentions because its practicability and useable both in military and civilian areas. In order to get more solar energy and better flight performance, HALE aircraft are normally built with relatively long and slender high-aspect ratio wings, with high lift-to-drag ratios to maximize aerodynamic efficiency. Because the limited available power, the dynamic pressure is relatively small, the weight of aircraft is limited so the wings of HALE aircraft are very flexible, with large wing deformations during flight. The natural frequencies are low enough so it tends to exhibit an overlap in the lowest elastic mode frequencies and the rigid-body flight dynamic frequencies and the coupling is serious. Researchers have made plenty of creative investigations on the subject with frequency and time-domain analysis by coupling nonlinear aerodynamic, structure and flight dynamics since the mishap of Helios since 2003. In the research of flexible wings, it is very important to get unsteady airloads accompanied with the deformation of wings with little calculated amount. Strip theory or vortex lattice method are used to get unsteady airloads in above works. The former combined with the dynamic stall model can simulate the aerodynamic force after the wing stall, but cannot give the three-dimensional aerodynamic effect of the wing. The latter is opposite. In this paper, two methods coupled geometric nonlinear beam model are established to establish the time domain analysis method of flexible wing respectively, to study the effect of dynamic stall and three-dimensional effect on the dynamic response of flexible wing. A medium aspect-ratio wings, Goland wing, which the aspect-ratio is 3.3 is studied to verified the 104

DDAP10, Huaqiao University, Xiamen China, 2018 models first. The velocity is increased until the displacement of wing is divergent. It can be found that the flutter speed is 143.1 m/s and 169.4 m/s separately, which coincide with others’ solutions . Because the velocity is small, the wing doesn’t reach the stall, the main reason of different flutter speed is the three-dimensional effect of UVLM. Then nonlinear aeroelasticity of a High-aspect-ratio wing, HALE wing, which the aspect-ratio is up to 16 is studied. For static equilibrium, the displacement of the wing in the tip is up to 44% of the span of the wing, which indicate the notable geometric nonlinear characteristic. In time domain simulation, the displacement of wing is not divergent like Goland wing when the inflow speed is larger than the flutter velocity, but vibrates with equal amplitude, which called limited loop oscillation (LCO). But it is the same with Goland wing that the the flutter of ONERA is lower than UVLM for HALE wing. The aerodynamic from ONERA over-valuate the outplane displacement, which is associated to Flutter velocity. ONERA modes underestimate the flutter velocity. Key words: Flexible aircraft, High-aspect-ratio wing, Aeroelasticity, Geometrical nonlinearity, Dynamic stall, UVLM, Flutter

PS-09

Square Superlattice Pattern in Dielectric Barrier Discharge Lifang Dong2 Yuyang Pan1 1College of Quality and Technical Supervision, Hebei University, Baoding 071002, China 21College of Physics Science and Technology, Hebei University, Baoding 071002, China E-mail: [email protected] and [email protected] Abstract Pattern formation is a self-organized phenomenon in nonlinear systems, which not only can be widely observed in natural systems such as animal coat markings, cloud formations and so on, but also can be studied in laboratory systems including thermal convection, nonlinear optics, chemical reactions, gas discharge, and so on. It is believed that the study of the pattern formation in gas discharge would advance the nonlinear science and plasma physics. Besides, there are plenty of potential applications about pattern formation, such as the plasma photonic crystal, the localized material growth, and so on. There are two water-electrodes in the experimental setup sealed with glass plates. A metallic ring is immersed in each container and connected to a sinusoidal ac power supply. A glass frame is clamped between the two parallel glass plates, serving as a lateral boundary. The whole cell is placed in a big chamber, where the gas parameters can be changed. An intensified charge-coupled device (ICCD) camera (HSFC pro) is applied to record frames from the end view of the electrodes. The camera has three similar, intensified photographing channels. An optical beam splitter is placed between the input lens and the intensified channels. An input beam is split into three similar beams by the beam splitter and then received by each intensified channel. Through the computer-controlled software, the exposure time of each channel can be changed, and the delay times between channels can also be set to snap frames at different times. With this setup, many types of superlattice patterns have been observed for the first time. By using high speed cameras it is found that the pattern in dielectric barrier discharge is a spatial-temporal pattern. Usually it is an interleaving of several sublattices, which discharge at different time. These results should greatly advance the science of pattern formation and plasma physics

105

DDAP10, Huaqiao University, Xiamen China, 2018

Key Words：pattern formation, dielectric barrier discharge, superlattice

PS-10

Smart chimera masks perturbation through self-adaptive drifting Nan Yao Xi'an University of Technology Email:[email protected] Chimera state is a remarkable spatiotemporal dynamical state since the structurally and dynamically identical oscillators in a coupled networked system spontaneously break into two groups, with one exhibiting coherent motion and another incoherent. In this paper, we report on a striking phenomena that chimera can smartly protect itself against the influence of perturbation through the so-called self-adaptive drifting. In the collective behavior of self-adaptive drifting, it always selects the optimal direction and the optimal final location for the chimera. We have systematically analyzed the spatiotemporal behavior of the drifting process of the smart chimera. Moreover, the exponential relaxation of chimera towards the optimal stable state collapses independently independent of the location of perturbation, which indicates a simple law underlying the collective dynamics. We have proposed an empirical model taking into account the damping and restoring force up on the chimera, which reproduces the numerical results and uncovers the dynamics of self-adaptive drifting in spatiotemporal pattern.

PS-11

Reciprocal Conditions in One-Dimensional Nonlinear System Hengzhe Yan, Gaomin Wu and Jie Ren Tongji University, Shanghai Email:[email protected] Abstract: Reciprocity means the invariance under source and detector exchange. Reicprocity has been proved in linear system with time-reversal symmetry while nonlinearity doesn’t always violate it. We investigate reciprocal conditions in a one-dimensional system consisting of multiple nonlienar δ-function potential seperated from each other by linear potential. A general sufficient and necessary reciprocal condition is proposed, in terms of an equation about transfer matrix. Based on this, we derive a geometric resonant condition. It indicates that reciprocity appears periodically by changing width of linear potential. Under resonant condition, the whole scatter behaves like a single nonlinear delta-barrier. We also discuss the existence of other reciprocal condition. For system with two nonlinear delta-barriers, we prove that it’s reciprocal only if it’s spatial symmetric or under resonant condition. And for system with more nonlienar barriers, reciprocity may apppear even if neither spatial symmetry nor resonant condition is fulfilled. We also explore magnitude reciprocal conditions that transitions amplitude’s magnitude is invariant while its phase change for system consisting of two nonlinear barriers. Both analytical and numerical results demonstrate the existence of such magnitude reciprocal conditions. Especially, for wide linear region with potential higher than energy, magnitude reciprocity appears for particular transmitted power, the approxmate value of which can be obtained using series expansion. 106

DDAP10, Huaqiao University, Xiamen China, 2018

Key Words:Reciprocal Conditions, Nonlinear

PS-12

Floquet topological acoustic resonators and acoustic Thouless pumping Yang Long, Jie Rena, Hong Chen Center for Phononics and Thermal Energy Science, China-EU Joint Center for Nanophononics, Shanghai Key Laboratory of Special Artificial Microstructure Materials and Technology, School of Physics Sciences and Engineering, Tongji University, Shanghai 200092, China [email protected] Abstract: Constructing the topological states can serve as novel approaches for robust acoustic wave transportations and manipulations. Here, we develop a scheme to realize acoustic topological states and adiabatic Thouless pumping in acoustic Floquet resonator systems. The acoustic wave can be manipulated adiabatically from one side to opposite side due to non-trivial topological nature. The physical mechanism behind these physical phenomena can be understood by effective one-dimensional Aubry-Andre Harper Hamiltonian under Floquet spatial periodic modulation, which is attributed to the projection from a high dimensional topological Hofstadter model for the Integer quantum Hall effect. Our scheme offers a new method for synthesizing acoustic topological states for efficient acoustic wave manipulations. Key Words: Acoustics, Topological states, Adiabatic pumping

PS-13

Topology and zero-index in 1D phononic metamaterials with negative density and negative modulus Danmei Zhang, Tianxiong Zhou and Jie Ren* Center for Phononics and Thermal Energy Science, China-EU Joint Center for Nanophononics, Shanghai Key Laboratory of Special Artificial Microstructure Mate- rials and Technology, School of Physics Sciences and Engineering, Tongji University, Shanghai 200092 *E-mail: [email protected]. Abstract: Metamaterials usually exhibit extraordinary properties that can not be find in nature like the negative permittivity, negative magnetic permeability of electromagnetical metamaterials which can realize reversed Doppler effect, reversed Cerenkov radiation and negative refraction index. Along with the electromagnetical metamaterials for photonics, the phononic metamaterials have attracted extensive attention since they are flexibly adjustable to control the phonon transmission. The mass-spring structure with different coupling and spatial distribution can also have negative effective parameters because of Bragg scattering and resonant mechanism in different frequency regimes. In 2007, Milton and Willis proposed a mass-in-mass system and showed that the dynamic effective mass of composite materials, defined in the framework of Newton’s law of motion,exhibits the existence of single or double negative properties. And in 2011, Liu et al. proposed an elastic model with double negative parameters by integrating a tri-chiral lattice with softly coated inclusions and an other double negative 107

DDAP10, Huaqiao University, Xiamen China, 2018 systems constructed by chiral mass-spring unit. Since then, a plenty of phononic metamaterials based on chiral and spiral structures are proposed. On the other hand, there are many other properties of phononic metamaterials attracting increasing interest recent years, such as the zero refraction index (zero-index for short) and topological bands. In zero-index metamaterials, waves does not carry any spatial phase changes and the wave length is effectively infinite long which can be used in wave surface modulation and bending waveguides. Many interesting topological phenomena, such as topological interface and edge states, have been observed in phononic metamatericals. In this work, we study a one-dimensional phononic metamaterial, made of mechanical resonant oscillators and chiral couplings. We show that by design the oscillator mass and inter-oscillator coupling, although both are positive naturally, can be either single negative or double negative effectively within a certain frequency range. At the frequency where the effective mass and coupling are both infinite, a flat band emerges that will induce an extremely high density of states which indicates that the band is in accordance with that caused from dark state. At the critical point of band degeneracy, a Dirac-like point emerges where both effective mass and the inverse of effective coupling are simultaneously zero, so that zero index is realized for phonons. Moreover, the phononic topological phase transition is observed that the phononic band gap switches between single mass-negative and single coupling-negative regime. As a consequence, a topological interface state is identified, well explained by the theory. These novel properties will provide good guidance for designers to obtain devices could have great application in wave propagation, for example shock absorption, damping noise and anti-interference and so on. Key Words: double-negative,

PS-14

Graph theory analysis of the quantum thermal switch Luqin Wang Tongji University 上海市四平路 1239 号 Email：[email protected] We demonstrate that a thermal switch can be composed by a quantum system of three interacting two-level systems, coupled to a thermal reservoir each. This thermal switch is analogous to an electronic bipolar one, which has the ability to control thermal currents at the collector and the emitter with the imposed current at the base. By tuning the temperature of the middle bath, we observe the switch effect and the heat amplification effect. Furthermore, we apply the graph theory to mainly analyze cyclic heat flux, edge flux and other features. They clearly exhibit which transition dominates in transport processes, which may guide to design thermal nano-devices.

PS-15

Non-equilibrium Restrictions on Nano-Scale Heat Pump by Fluctuation Theorem Zi Wang(王子), Jie Ren(任捷)* Center for Phononics and Thermal Energy Science, China-EU Joint Center for Nanophononics, School of Physics Sciences and Engineering, Tongji University, Shanghai 200092 *E-mail: [email protected] 108

DDAP10, Huaqiao University, Xiamen China, 2018

Abstract: A great deal of physical applications and insights has been provided by the elegant and successful theory of linear response, which provides us an understanding of the ubiquitous phenomena regarding action and response in the linear regime. The surge of marching into the more generic non-equilibrium region has been on for nearly two decades, since the discovery of Jarzynski’s work relation, providing information beyond the linear response approximation. The series of established constructions is called fluctuation theorem in general. The non-equilibrium thermodynamics concerning heat can be studied from this perspective. We adopted a stochastic formalism, including Langevin equation, Fokker-Planck equation, et al., to study the heat pump system both in contact of heat reservoir and in manipulation of work protocol, giving the relation of heat flow and work injected. Some restrictions will be given in model systems, which is indifferent to the specific details of system studied. These relations pertain robustly in the presence of far non-equilibrium state and strong non-linearity. Key Words: Fluctuation Theorem, Heat Conduction, Nano-pump

PS-16

Detecting Directed Interactions of Networks by Random Variable Resetting Rundong Shi, Changbao Deng and Shihong Wang School of Sciences, Beijing University of Posts and Telecommunications - Beijing 100876,China Email:[email protected] Abstract We propose a novel method of detecting directed interactions of a general dynamic network from measured data. By repeating random state variable resetting of a target node and appropriately averaging over the measurable data, the pairwise coupling function between the target and the response nodes is constructed. This method is applicable to a wide class of networks with nonlinear dynamics, hidden variables and strong noise. The numerical results have fully verified the validity of the theoretical derivation. Key Word Time series analysis, network reconstruction, nonlinear system, complex network

PS-17

Entropy production and infima statistics in the granular Brownian rotor Ke Cheng, Lanzhou University Tianshui South Road No. 222, Lanzhou, Gansu province Email:[email protected] Abstract: We experimentally studied the relation of trajectory entropy and thermal entropy, the fluctuation theorem and the infima law of entropy production in a non-equilibrium system. It is well known that in an isolated system, entropy will not decrease on average, but the entropy's increment is not smooth because of the fluctuations in the environment. Those fluctuations will lead to the negative entropy production within a short time. Therefore, the following question is important: what are the extreme values of negative entropy production? A recent work [1] shows an out-of-ordinary conclusion. In their investigation, entropy production in a non-equilibrium state is regarded as a martingale. The increasing 109

DDAP10, Huaqiao University, Xiamen China, 2018 entropy has an infimum, which is the negative limit in a time interval t. One universal conclusion is that the average value of this negative limit over a large number of time intervals is greater or equal than minus the Boltzmann constant. This conclusion gives a character of entropy: no matter how fast the entropy production rate is, the global infimum maintains a constant. They also provide an over-damped drift-diffusion process and make an analytical solution for the infimum of the entropy production. We set up a horizontal Brownian rotor on a thin vibrated container, and calculate the trajectory entropy, the thermal entropy, and total entropy production. The experiment is Our experiment gives a result that the global infima of entropy production in an under-damped Brownian rotor would exceed minus the Boltzmann constant if the sampling time is much smaller than the system characteristic time, and then find a scaling law of global infima value and sampling time. Key words: entropy production, fluctuation theorem, infima law.

PS-18

Complex Dynamics in the Carbon Nanotubes：A Molecular Dynamics Study Yisen Wang, and Liang Huang Lanzhou University, Lanzhou, Gansu 730000, China Email:[email protected] Abstract In recent years, the nonlinear dynamics of carbon nanotubes has drawn much attention of scientists/researchers due to its unique mechanical properties and significant application potentials ranging from ultrasensitive mass sensor to single-electron spin detection. In this work, the nonlinear and chaotic motions of the single wall carbon nanotubes (SWCNT) subjected to an external driven force are investigated with extensive molecular dynamics (MD) simulations. The SWCNT has a length of L = 30nm, and is assumed to clamped at both ends, with the next four layers being taken as thermostat atoms and the remaining are treated as Newtonian atoms. After the initial relaxation, a uniform driven force is applied to the SWCNT, and the trajectory of the center of the SWCNT is used to investigate the dynamical behavior of the system. The numerical results show that, for the small driven force, the dynamics appear to be linear and the SWCNT oscillates within a fixed plane. However, as the driven force is sufficiently large to yield substantial bending and consequently strong axial stress, the planar motion is no longer stable and a whirling or “jump-rope” like motion occurs due to the nonlinear interaction between orthogonal vibrational modes, which is consistent with the previous theoretical results [1, 2]. Furthermore, rich dynamical behaviors have been identified towards the development of chaos and the transition between different types of chaotic motions. Key Words: Carbon Nanotubes, Nonlinear Dynamics, Nonplanar Motions, Chaotic Motion Reference [1] W.G. Conley, A. Raman, C.M. Krousgrill, and S. Mohammadi, Nano Lett. 8, 1590 (2008) [2] Q. Chen, L. Huang, Y.-C. Lai, C. Grebogi, and D. Dietz, Nano Lett. 10, 406 (2010)

PS-19

Non-Gaussian noise-weakened stability in a foraging colony system with time delay 110

DDAP10, Huaqiao University, Xiamen China, 2018

Xiaohui Dong1,2(董晓慧)，Chunhua Zeng1(曾春华) 1.Faculty of Science, Kunming University of Science and Technology, Kunming 650093, Yunnan, China 2.Center of soft matter physics and its applications, Beihang University, Beijing 100191, China Email：[email protected] Abstract: In this paper, we studied the dynamical properties of a foraging colony system with time delay and non-Gaussian noise. Using the delay Fokker-Planck approach, the stationary probability distribution (SPD), the associated relaxation time (ART) and normalization correlation function (NCF) are obtained, respectively. The results show that: (i) the time delay and non-Gaussian noise can induce a transition from a single peak to double peaks in the SPD. The further numerical simulations are consistened with the approximate theoretical results; (ii) there exists a maximum in the ART as a function of the noise intensity. This maximum of the ART is identified as the characteristic of the non-Gaussian noise-weakened stability of the foraging colonies in the steady state; and (iii) the time delay and departure from Gaussian noise play opposite roles in ART or NCF. Key Words: time delay; non-Gaussian noise; foraging colony system; dynamical properties

PS-20

Dynamics of recurrence in Non-integral Nonlinear Schrodinger Equations Zhigang Zhu Lanzhou University 甘肃省兰州市天水南路 222 号兰州大学 Email：[email protected] The mode exciting dynamics of nonlinear Schrodinger equations in non-integral systems exhibit exponential instability. We give an analytical expression to characteristic this instability. It also exhibit recurrence phenomenon at the beginning of excitation. We also provide an reduced model equation to formulate this phenomenon.

PS-21

Synchronous patterns in coupled non-identical chaotic oscillators Yafeng Wang, and Xingang Wang School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710062, China Email:[email protected] Abstract: A universal feature of the spatiotemporal nonlinear systems in nature is that under few certain conditions they can self-organized into various patterns, which has arouse the extensive studies of pattern formation over the past few decades. In general, the considered systems are homogeneous in space and the generated patterns are the regular structures, e.g., the model of coupled map lattice. Recently, inspired by the rapid progress of network science, more attentions have been also paid to the formation of patterns in inhomogeneous systems represented by complex networks. Synchronizaiton, as a well-studies topic, describes an inherntly collective dynamic behaviors of coupled nonlinear systems and has attracted continuous interest in nonlinear science. Complete synchronization evolves 111

DDAP10, Huaqiao University, Xiamen China, 2018 on the same chaotic trajectory starting with the random initial conditions in the phase space when the coupling strength among the oscillators is larger than the critical value. Simultaneously, partial synchronization has been also broadly interested and extensively studied, for its strong connections to the functioning of many biological and physiological systems. In a partial synchronization state, oscillators inside a synchronous cluster are highly correlated, while they are loosely or not correlated if the oscillators belong to different clusters. The specific form of the synchronous clusters, namely the synchronous pattern, depends not only on the dynamics of the oscillators, but also on the coupling function and coupling strength. Many studies have demonstrated the generation of the synchronous patterns is closely dependent on the network topological symmetries. In particular, using the technique of computational group theory, Pecora et al. are able to identify all the permutation symmetries of a complex network and, based on the method of eigenvalue analysis, predict the stable synchronous patterns that can generated form the random initial conditions. In the synchronization of coupled nonlinear oscillators, a general finding is that as the oscillators become non-identical, the propensity for synchronization will be deteriorated. Yet, parameter mismatch is unavoidable in realistic complex systems and, in some cases, the oscillators are even of different types of dynamics. The contradiction between theory and reality poses the following paradox for oscillator synchronization: how non-identical oscillators are synchronized in a complex network? Here, by the model of coupled chaotic oscillators, we are able to demonstrate numerically and argue analytically that, despite the parameter mismatch (different dynamics) of the oscillators, stable synchronous patterns can still be generated in complex networks. A general framework is proposed for the identification and stability analysis of synchronous patterns in the general complex networks, with the theoretical predications in good agreement with the numerical results. The study sheds lights on the collective behaviors of coupled non-identical oscillators, and might have implications the functionality and operation of many realistic systems, e.g., the power-grid and the human brain network. Key Words: Synchronous patterns, Non-identical chaotic oscillators, complex network

PS-22

An abnormal phenomenon in pinning synchronization Dapeng Zhang1, Xingang Wang1，* *: [email protected] 1 School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710062, China Abstract Synchronization, the most common phenomenon in nature, have been extensively studied in nonlinear science including physics, optics, chemistry and biology and investigated in realistic networks or man-made systems. To answer questions of complete synchronization (CS) from the research of behaviors of coupled identical chaotic oscillators, which can trace back to the work of Pecora and Carroll in 1990. In addition, the critical condition for CS can be analyzed by the method of master stability function (MSF). MSF, as a power tool, can transform the coupled oscillators into isolated form and synchronization of system can be determined by the following two factors: the MSF curve (got by the oscillator dynamics and the coupling function) and the eigenvalues (calculated by the coupling matrix of system). With the discovery and study of small-world networks and scale-free networks on synchronization, which means most of the synchronization observed in regular networks 112

DDAP10, Huaqiao University, Xiamen China, 2018 can occur in complex networks. As times goes on, the study of synchronization is more detailed and divided into different types, such as partial (cluster) synchronization (PS). PS means that the evolving of oscillators within the same cluster converges to synchronous manifold of the cluster, but not for other oscillators from the network. The primary problem in the study of PS in complex network is the identification of clusters, while a large amount of evidences show that PS is directly related to the symmetry of the network. In Ref. [2], Pecora et al. identify all the permutation symmetries numerically in a given complex network by means of computational group theory, based on the generalized method of eigenvalue analysis. To achieve the purpose of synchronization for any cluster, based on the stability analysis of PS, the necessary condition is that the perturbations transverse to the synchronous manifold of cluster should be damping to 0 with time. According to MSF, PS can be treated as a problem of global synchronization and analyzed when one cluster is not related to another. Thus, the key point is that all Lyapunov exponents of transverse perturbations in cluster are negative. In recent years, inducing PS states in complex networks by pinning control (an extra controller) have been studied. In this paper, focusing on inducing PS by pinning targeted clusters, we will present an abnormal phenomenon: synchronization occurs for controlled cluster, at the same time, which may occur for uncontrolled cluster. Here, based on the method of Lyapunov exponents of transverse manifold for synchronous clusters, we are able to figure out above-mentioned cases. In particular, we find that this method is consistent with simulation and inconsistent sometimes. Aiming at this phenomenon, we further give the detailed information for the synchronous mechanism of cluster. Our results demonstrate that the emergence of synchronization for uncontrolled cluster is from tremendous changes of manifold, by the influence of the parameters of system. Moreover, similar phenomenon can be identified in map. Key words: partial synchronization; pinning control; Lyapunov exponents

PS-23

Pinning cluster synchronization in regular networks Liang Wang1, Xingang Wang1,* 1 School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710062, China * [email protected] Abstract: Synchronous behaviors are commonly observed in natural and man-made systems and widely recognized as important to the system functionality and operations. The complete synchronization in the network is that all nodes have reached the same state. In theory, it can be analyzed by the method of MSF, which can give a good prediction. In the recent research about chaos synchronization in networks, partial synchronization has attracted wide attention, which is a special synchronous state before complete synchronization. It is straightforward to find that partial synchronization of network depends on its topological symmetry, that is to say, every symmetry potentially supports a synchronous pattern. Although of practical significance, the state of partial synchronization, potentially supported by the symmetry, is not necessarily able to appear due to its instability, such as the study of Chimera state, it helps us to understand the unihemispheric sleep of dolphins and other birds. If a network is unable to spontaneously reach synchronization, an additional controller can be introduced to help network achieve synchronization. The controller affects the whole network. Therefore, we introduce a controller, based on the symmetry of the original network, to control some nodes by the method of pinning control 113

DDAP10, Huaqiao University, Xiamen China, 2018 and effectively generate the cluster synchronous patterns. We take the global network and ring network as examples. The regular networks generally have high symmetry, which support abundant synchronous patterns theoretically , however, with the increase of coupling strength, many potential synchronous patterns can not appear because of instability. Therefor, a controller is introduced to unidirectionally control all nodes within a cluster supported by a kind of symmetry of the original network, and the controller is considered as a special node in an enlarged network. Due to the influence of the controller, the original symmetry of the network will decline, and the new phenomenon of cluster synchronization supported by the symmetry will appear. This phenomenon can be predicted by the corresponding theoretical method. The feasibility and efficiency of the control method are verified by numerical simulations. In addition to, we also find that it is impossible to make the pinned oscillators synchronize with the controller except the synchronization of the whole network. Of course, our simulation and theory will show large deviation in some conditions, which is caused by the change of nodal manifold through analysis, and the greater the controlling strength, the greater the manifold change. Our research provides a way to generate the cluster synchronous patterns. Although it is proposed on the regular networks, it can also be extended to the complex networks. Key words: synchronization, pinning control, regular network

PS-24

The Ising Model on Signed Networks Lingbo Li, Beijing Normal University Ying Fan, Beijing Normal University Email:[email protected] Abstract Complex networks can accurately outline and describe real systems in nature and human society, and provide more effective methods to express the structure and relationship of natural and social networks than other methods in the past. In order to meet more specific and adaptive needs of description, the complex network theory is also developing and improving. For instance, in signed networks, positive and negative edges have different effects on dynamic behaviors among individuals, but the specific role that the negative edges play in the process is still unclear, which is the focus and difficulty of the research about signed networks. The dynamic methods have been used to study the infection and evolution mechanism of public opinion on the complex network structure, which can effectively describe and explain the interacting process of individual views, so as to observe the macro emergence of public opinion. The Ising model is the beginning of public opinion evolution research by dynamic methods. This paper applies the classical Ising model to the study of public opinion infection and evolution on signed networks. The particles in the Ising model are placed on the nodes of the signed network to represent the opinions of nodes, and the opinions of nodes are often influenced or even changed by the friends around them. The couple of nodes connected by the positive edge have positive relation, who tend to hold the same opinion, while the couple of nodes connected by the negative edge have negative relation, who are more inclined to form the opposite view. According to the rules of the Ising model, the total energy of the system is proposed to measure the state of the system. When the total energy of is lowest, the system is in the most stable state. Therefore, whatever the initial state of the system is, it 114

DDAP10, Huaqiao University, Xiamen China, 2018 will constantly evolve and change the specific distribution of opinions of nodes in the network, so as to constantly reduce the total energy value of the system. The mutual influence between the nodes results in the final macro emergence, thus showing the phenomenon of public opinion polarization in the actual system. This study found that network structure characteristics affected by the proportion and distribution of negative sides, have a fundamental effect on the result of public opinion infection and evolution on signed networks. There exists the critical ratio, which amount to a physical transformation point. When the proportion of negative edges in the network over the critical ratio, the evolution results appear phase transition, and the distribution of negative edges affects the value of the critical ratio. In addition, this report also complements the research on the network balance in the process of the public opinion evolution on the signed networks. Key Words: signed network, Ising Model, negative edges, opinion spreading, structural balance

PS-25

A class of chaotic oscillators dependent on the initial selection Jun Ma Ge Zhang Fuqiang Wu Zhilong Liu Wenkang Xu Department of Physics, Lanzhou University of Technology, Lanzhou 730050, China [email protected] Abstract Nonlinear dynamics composed of the state variables and control parameters have been investigated by the bifurcation theory. The transition between chaotic and periodic atractors in the system depends on the selection of the control parameters and initial values. In this paper, we propose a class of nonlinear system with memory effect by adding a nonlinear term such as z2y. The nonlinear system can trigger a different profile of attractors by setting different initial values. The emergence of chaos and state selection has been confirmed by methods, including the phase portrait, Lyapunov exponent spectrum, Hamilton energy and bifurcation analysis. In the case of network synchronization and pattern selection, a chain network is reconstructed of the oscillator with hidden attractor and memory effect, and its statistical factor of synchronization is calculated to predict the synchronization stability. It is found that the synchronization stability shows some dependence on initial setting for state variable. When some nodes in a chain network have the memory effect, enhancement of memory function can decrease the synchronization, while memory function in a small region can contribute to the synchronization. And the Hamilton energy is also dependent on initial setting and control parameters based on Helmholtz theorem. We also design a corresponding circuit based on the field programmable gate array (FPGA). The nonlinear system with memory dependence on initial setting was verified in digital signal processing (DSP) builder block under Matlab/Simulink. Key words: Synchronization; Bifurcation; Network; Field programmable gate array (FPGA)

PS-26

Pattern transition of firing in a Morris-Lecar neuron induced autapse Xinlin Song (宋欣林)1,2, Hengtong Wang (王恒通)3, Yong Chen (陈勇)1,2 115

DDAP10, Huaqiao University, Xiamen China, 2018

1.Center of soft matter physics and its applications, Beihang University, Beijing 100191, China 2.School of Physics and Nuclear Energy Engineering, Beihang University, Beijing 100191, China 3.School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China Email:[email protected] Abstract: It is well-known that there are three firing patterns of a Morris-Lecar (ML) neuron under the external DC input. In this work, we investigated the transition between the firing patterns induced by autapse in the ML neuron. We presented the frequency band ranges dependent on the autaptic time delay and the phase response curves. Thus, we observed that the firing pattern could transform from the Class-1 to the Class-2 under the excitatory or inhibitory autaptic feedback. Importantly, It was founded that for the Class-3 firing pattern, the excitatory autapse based on the single spike could exist at some autaptic time delay in some range of phase plane. This means that the Clas-3 ML neuron could behave the Class-2 firing pattern, not Class-1 firing. Key Words: Autapse; Morris-Lecar neuron, firing pattern, time delay

PS-27

Anomalous interfacial temperature profile induced by phonon localization Yue Liu and Dahai He Department of Physics and Institute of Theoretical Physics and Astrophysics, Xiamen University, Xiamen 361005, Fujian, China Email:[email protected] Abstract: Through the integration of the power spectral density, we obtain temperature profiles of both multisegment harmonic and anharmonic systems, showing the presence of an anomalous negative temperature gradient inside the interfacial segment. Via investigating patterns of the power spectral density, we found that the counterintuitive phenomenon comes from the presence of interfacial localized phonon modes. Two out-band localized modes of the harmonic model, which make no contributions to local temperature due to the absence of phonon interactions, result in the concave temperature profile and overcooling effect. For the anharmonic model, thanks to the phonon-phonon interactions, the localized modes are excited and make considerable contributions to interfacial temperature, which is clearly shown by examining the temperature accumulation function. When anharmonicity is considerably large, the negative temperature gradient is absent since the localized phonon modes are fully mixed. The presence of localized modes are evidently demonstrated by the inverse participation ratio and normal mode analysis for the isolated harmonic model. The localized modes make contribution to interfacial temperature profiles of the harmonic system when they are excited in initial conditions of simulations. Key words: Nonequilibrium statistical mechanics, phonon localized, spectral decompositions of temperature.

PS-28

Fractal Aggregates on Geometric Graphs 116

DDAP10, Huaqiao University, Xiamen China, 2018

Ling Zhou, Sheng-jun Wang School of Physics and Information Technology Shaanxi Normal University, Xi’an 710119, China Email: [email protected] Abstract: The aggregation process is of interest in many areas of science. The most known growth model is the diffusion-limited aggregation (DLA). The DLA and its variants can recur in many fractal patterns in diffusive systems. In the DLA model, the particle jumps from the current site to one of their nearest ( lattice site at each step. However, some realistic systems have random media, that is, sites do not array regularly on lattice but randomly distribute in space. A typical example is the porous media. In this case, the randomness can be modeled by a geometric graph. a) Geometric graphs offer a way to model a lot of networks and infrastructures found in the real world. Water distribution, telephone services and water flow in a stream are all examples of resource flows that can be modeled and analyzed using a geometric graph.The geometric graph is probably the simplest model of the spatial network, is a graph in which the sites or edges are spatial elements associated with geometric objects. The sites are located in a space equipped with a certain metric and the distribution has a certain randomness, which has important effects on their topological properties. There is a large body of literature studying the effect of randomness of geometric graph. Some dynamical behaviors are sensitive to the randomness. In our work, we will focus on the aggregation process on the geometric graph by the method of computer simulation analysis, and investigate the effect of randomness of the network structure. The numerical simulation results indicate that the structure of the aggregate is similar to that of the DLA model on the regular lattice. Density correlations within the aggregates fall off with distance with a fractional power law. It indicates that the aggregate pattern on the geometric graph is fractal. The fractal is robust against the randomness in the structure. We calculated the fractal dimension of aggregates. The coarse-grained grid affects the fractal dimension of aggregates. Furthermore, the connection degree increases monotonically as the radius of local area increases. We have found that the aggregate and the exponent of the density correlation function depend on the connection degree of network. The fractal dimension can be adjusted by changing the connection degree of the geometric graph. Keywords: Fractal, Aggregate, Geometric Graph, DLA

PS-29

Effects of inhibitory signal on criticality in excitatory-inhibitory networks ( Fan Wang and Sheng-Jun Wang School of Physics and Information Technology Shaanxi Normal University, Xi’an 710119, P. R. China b) Email: [email protected] Abstract: Self-organized critical states are found for many complex systems in nature, from earthquakes to neuronal avalanches. Several lines of evidence point to the existence of such critical states in brain ( activity. The excitable networks can be used as a simplified model to study the criticality in neural networks. In realistic neural systems, the excitatory and inhibitory neurons are coexisting, and inhibitory is always included in neural dynamical models. The behavior of excitatory-inhibitory (EI) a) 117

DDAP10, Huaqiao University, Xiamen China, 2018 network is critical for understanding how neural circuits produce cognitive function. We study the criticality in excitable networks consisting both excitatory and inhibitory elements. In general, the change of E/I signal strength was modeled by increasing or decreasing the branching ratio. In this present work, we concern the effect of the strength of the inhibitory elements on criticality and sensibility of network. The effects of the inhibitory coupling strength on the criticality are investigated numerically and analytically. We consider a model of EI network on the Erdös-Rényi (ER) graphs using the model proposed by Kinouchi and Copelli. The signal transferred from an excitatory node increases the probability that the neighbors of this node are excited, while the signal from an inhibitory node decreases this probability. To study the effects of the strength of inhibitory coupling in the model, we set the maximal inhibitory probability is several times than maximal excitatory probability. We show that the inhibitory coupling strength does not affect the critical point. The dynamic range is decreased as the inhibitory coupling strength increases. Through an analytic treatment with mean-field approximation, we show that in excitable networks the inhibitory signal is a feedback of excitatory signal and only play an important role in high activity state. In the critical state the mean activity is small and the inhibitory feedback is ignorable. Our results show that the behaviors of excitable network model with inhibitory elements are different from the biological experiments [1], and more sophisticated modeling studies are needed to understand the criticality in EI networks. These studies provide new insight into the role of inhibitory elements in realistic neural networks. The analysis is also interesting for other many systems which can be studied using excitable networks. [1] Shew W L, Yang H, Petermann T, Roy R and Plenz D 2009 Journal of Neuroscience 29 15595– 15600 Keywords: Self-organized criticality, excitable network, dynamic range

PS-30

Heat current through a weak interface Yanjiang Guo Department of Physics, Renmin University of China Email:[email protected] We systematically study heat current J that flows through a few one-dimensional nonlinear lattices, each of which consists of two identical segments that are coupled by a weak interface. And we found that J follows a power-law decay with the strength and the power exponent depends on the details of the lattices and the interface interaction.

PS-31

Persistence paves the way for cooperation in evolutionary games Changwei Huang Email：[email protected] Beijing University of Posts and Telecommunications Abstract: Cooperation is an effective way to maximize collective benefits, especially in modern human society. The issues on the emergence and maintenance of cooperation have attracted much attention in recent years. Here, we introduce the persistence parameter \tau to characterize the time duration of 118

DDAP10, Huaqiao University, Xiamen China, 2018 choices held by individuals and consider the effects of \tau on cooperation. We find that persistence could promote cooperation in a population no matter what the network structure is. Furthermore, the results on heterogeneous networks show that individuals with larger \tau are more inclined to cooperate than those with smaller \tau. Moreover, we investigate the effects of correlations between degree and persistence in scale-free networks and find that assortative matching could remarkably enhance cooperation whereas disassortative matching has adverse impacts on the evolution of cooperation. Key words: Strategy persistence, Cooperation, Complex networks.

PS-32

The spreading dynamics of forget-remember mechanism S.F. Deng and W. Li Hua-Zhong Normal University, Wuhan 430079, China Abstract: We study extensively the forget-remember mechanism (FRM) for message spreading, originally introduced in Eur. Phys. J. B 62, 247 (2008). The freedom of specifying forget-remember functions governing the FRM can enrich the spreading dynamics to a very large extent. The master equation is derived for describing the FRM dynamics. By applying the mean field techniques, we have shown how the steady states can be reached under certain conditions, which agrees well with the Monte Carlo simulations. The distributions of forget and remember times can be explicitly given when the forget-remember functions take linear or exponential forms, which might shed some light on understanding the temporal nature of diseases like flu. For time-dependent FRM there is an epidemic threshold related to the FRM parameters. We have proven that the mean field critical transmissibility for the SIS model and the critical transmissibility for the SIR model are the lower and the upper bounds of the critical transmissibility for the FRM model, respectively. Keywords: spreading dynamics; percolation; phase transition Reference: PHYSICAL REVIEW E 95, 042306 (2017)

PS-33

Self-avoiding pruning random walk on signed networks Huijuan Wang2, Cunquan Qu1,2,*, Chongze Jiao2, and Wioletta Rusze2 1) School of Mathematics, Shandong University, 27 Shanda Nanlu, Ji’nan, P.R.China, 250100 2) Faculty of Electrical Engineering, Mathematics, and Computer Science,Delft University of Technology, Mekelweg 4, Delft, The Netherlands, 2628 CD.

Abstract A signed network represents how a set of nodes are connected by two logically contradictory types of links: positive and negative links. In a signed products network, two products can be complementary (purchased together) or substitutable (purchased instead of each other). Such contradictory types of links may play dramatically different roles in the spreading process of information, opinion, behavior etc. In this work, we propose a Self-Avoiding Pruning (SAP) random walk on a signed network to model e.g. a user's purchase activity on a signed products network. A SAP walk starts at a random node.

119

DDAP10, Huaqiao University, Xiamen China, 2018

At each step, the walker moves to a positive neighbour that is randomly selected and its previously visited node together with its negative neighbours are removed. We explored both analytically and numerically how signed network topological features influence the key performance of a SAP walk: the evolution of the pruned network resulted from the node removals, the length of a SAP walk and the visiting probability of each node. These findings in signed network models are further verified in two real-world signed networks. Our findings may inspire the design of recommender systems regarding how recommendations and competitions may influence consumers' purchases and products' popularity. Key words: Signed Network, random walk, self-avoiding pruning walk, network pruning, walk length, visiting probability

PS-34

Navigation with and without neural representation Tie Xu, Israel Institute of Technology Abstract The neural dynamics and navigation，we study the property of coupled dynamical system of recurrent neural network and environment, showing the emergence of spatial representation in random network.

PS-35

Measuring Lyapunov spectrum of large chaotic systems with global coupling by time series analysis Taro P. Shimizu1 and Takeuchi A. Kazumasa2,1 1 Department of Physics, Tokyo Institute of Technology. 2 Department of Physics, The University of Tokyo. Chaos is often characterized by the Lyapunov exponent. The Lyapunov exponents are related to fundamental physical quantities of large chaotic systems, such as the attractor dimension and the metric entropy. Moreover, the Lyapunov exponents and related quantities play an important role in chaos control and data assimilation. Experimentally, one has to estimate the Lyapunov exponents from time series since the equation of motion is usually unknown. The most standard approach to evaluating the Lyapunov exponents from time series is regarding recurrences of time series as perturbations on trajectories and measuring the exponential growth rate of their norms. However, for large systems, it remains practically impossible to evaluate even the maximal Lyapunov exponent by this recurrence method. This is because recurrence becomes extremely rare if the number of degrees of freedom, hence the dimensionality of phase space, is large. Recently, by using a machine learning technique, Pathak et al. succeeded in predicting trajectories and the Lyapunov exponents of a large spatially extended system [1]. However, adjusting various parameters involved in this method without guiding principles is an issue to be solved in practice. In this work, we attempt to extend the recurrence method and to solve the issue of the lack of recurrences in large systems. Here we restrict our target to systems with global coupling. We focus on the fact that the time evolution of these systems is determined by the local variable and the mean field. This enables us to reduce the effective dimensionality for detecting recurrent events. We, therefore, 120

DDAP10, Huaqiao University, Xiamen China, 2018 gather recurrences with this local set of variables and show that it is sufficient to construct the Jacobian matrix, which is necessary to compute the full spectrum of the Lyapunov exponents. We apply our method to two globally-coupled systems, namely, logistic maps and limit-cycle oscillators with global coupling, and demonstrate that this method is able to evaluate the Lyapunov spectrum reasonably well. In this poster presentation, we will show you the details of our method and simulation results [2]. References [1] J. Pathak et al., Chaos 27, 121102 (2017); Phys. Rev. Lett. 120, 024102 (2018). [2] T. P. Shimizu and K. A. Takeuchi, arXiv:1810.01146 (2018).

PS-36

Chaotic dynamics of the Willis aneurysm system with Atangana–Baleanu derivative with fractional order and its control Gao Fei, Li Wen-qin, Li Xi-ling School of Science, Wuhan University of Technology, Wuhan 430070 Abstract: The Willis aneurysm system has some limitations in the description of the complex hemodynamic mechanism of blood with viscoelasticity. The fractional calculus has been used to depict some complex and disordered processes in organisms. Thus, we propose a fractional Willis aneurysm system (ABCWAS) by using the Atangana-Beleanu fractional differential and its theory in the present article. we prove that the ABCWAS has a chaotic characteristic by analyzing the phase portraits and Poincaré section, and it is a rational extension of its integer order form. We investigate the influences of fractional order on the ABCWAS by means of bifurcation diagram and period doubling bifurcation. The results show that small changes of fractional order can lead to a remarkable effect on the motion state of the ABCWAS. As the chaotic ABCWAS indicates that the brain blood flow is unstable, and the cerebral aneurysms are more likely to rupture in a very chaotic velocity field. Therefore we use a methods to control the chaotic ABCWAS. The methods is to design a suitable controller based on the stability theorem of fractional nonlinear non-autonomous system. The numerical simulations show that the proposed the methods can control the blood flow velocity and speed up the periodic fluctuation within a small range, which shows that the cerebral aneurysm is not easy to rupture. The theoretical results in our article can provide some theoretical guidance for controlling and utilizing the actual ABCWAS system. Key words: fractional Willis aneurysm system, Atangana-Beleanu fractional differential, Poincaré section, chaos control

121

DDAP10, Huaqiao University, Xiamen China, 2018

PS-37

Coexisting attractors in chaotic system of conditional symmetry Tianai Lu1, 2 Chunbiao Li1, 2 a) 1 Jiangsu Collaborative Innovation Center of Atmospheric Environment and Equipment Technology (CICAEET), Nanjing University of Information Science & Technology, Nanjing 210044, China 2Jiangsu Key Laboratory of Meteorological Observation and Information Processing, Nanjing University of Information Science and Technology, Nanjing 210044, China Sajad Jafari Biomedical Engineering Faculty, Amirkabir University of Technology, 424 Hafez Ave, 15875-4413, Tehran, Iran Fuhong Min School of Electrical and Automation Engineering, Nanjing Normal University, Nanjing 210042, China ABSTRACT Conditional symmetry is known as a new regime for providing coexisting attractors. The polarity of some of the feedback in a dynamical system can be revised for hosting conditional symmetry. In this paper, new chaotic systems of conditional symmetry are coined from 1-D and 2-D offset boosting based on a previous polarity adjustment. Conditional symmetric or asymmetric coexisting attractors are found accordingly. Moreover, the distance between two coexisting attractors of conditional symmetry can be modified linearly by the offset boosting constant while the size of the coexisting attractors can also be controlled by a newly introduced slope. Coexisting attractors of conditional symmetry are also demonstrated by the implementation based on STM32. Keywords: Multistability, conditional symmetry , offset boosting, STM32

a) Corresponding author, E-mail address: [email protected]; [email protected].

122

DDAP10, Huaqiao University, Xiamen China, 2018

VI. Contacts

The conference email address is [email protected]. You can also contact the following members for any information about the conference.

Name Email Zhigang Zheng 郑志刚 [email protected] Nianbei Li 李念北 [email protected] Hongbin Chen 陈宏斌 [email protected] Yun Zhai 翟云 [email protected]

For Hotels, please contact: Hotel Contact North Bay Wanda Hotel 北海湾惠龙万达酒店叶腾华(Ye Teng-hua) 17805980420 Ling Ling Hotel 灵玲大酒店林惠丽(Lin Hui-li)15980989148 Mingzhu Harbor Hotel 明珠海湾大酒店叶丹萍(Ye Dan-ping)15280119161

Institute of Systems College of Information Science and Engineering Huaqiao University

华侨大学系统科学研究所信息科学与工程学院

123

DDAP10, Huaqiao University, Xiamen China, 2018

VII. Acknowledgement

The DDAP Committee acknowledges the following sponsors for the support：

 Chaos, AIP Publishing

 Nonlinearity, IOP Publishing

 Ye-Bin Jian-Zhan Bowls Culture Co. Ltd, Nanping, Fujian

(福建南平叶镔建盏文化传播有限公司)

124