CN111721693A - 基于分子动力学模拟的沥青抗老化性能评价方法 - Google Patents

基于分子动力学模拟的沥青抗老化性能评价方法 Download PDF

Info

Publication number
CN111721693A
CN111721693A CN202010645196.5A CN202010645196A CN111721693A CN 111721693 A CN111721693 A CN 111721693A CN 202010645196 A CN202010645196 A CN 202010645196A CN 111721693 A CN111721693 A CN 111721693A
Authority
CN
China
Prior art keywords
asphalt
aging
simulation
molecular dynamics
dynamics simulation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010645196.5A
Other languages
English (en)
Inventor
顾兴宇
胡栋梁
崔冰彦
周洲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202010645196.5A priority Critical patent/CN111721693A/zh
Publication of CN111721693A publication Critical patent/CN111721693A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N17/00Investigating resistance of materials to the weather, to corrosion, or to light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N5/00Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N2021/3595Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR

Abstract

基于分子动力学模拟的沥青抗老化性能评价方法,本发明的目的在于克服现有试验方法评价沥青抗老化性能研究效果的不足,通过开展基于反应力场的分子动力学模拟,提出一种理论依据可靠、结果简洁明确的方法评价沥青的抗老化性能。本发明采用基于反应力场的分子动力学模拟方法,通过精确量化沥青老化后含氧官能团的数量,直接反映沥青抗老化性能,过程规范,结果可靠,克服了实验研究方法无法准确量化或者不能直接反映沥青抗老化性能的不足之处;并且本发明采用的分子动力学模拟方法所用力场参数根据第一性原理计算的结果优化得到,反映了沥青在抗老化方面本身的物理化学属性,不受试验条件的干扰。

Description

基于分子动力学模拟的沥青抗老化性能评价方法
技术领域
本发明属于沥青抗老化性能评价技术领域,具体涉及一种基于分子动力学模拟的沥青抗老化性能评价方法。
背景技术
沥青路面在使用阶段长期处于大气环境下,因此沥青容易发生老化,沥青老化是影响沥青路面使用性能的重要因素之一。沥青老化后的宏观表现为针入度下降,软化点上升,延度减小,这些物理指标的劣化导致沥青路面强度和耐久性降低,容易引起沥青路面开裂与剥落。由于路面所处的自然环境无法改变,因此只能通过提高沥青的抗老化性能以防止沥青路面过早地出现老化相关病害。
为了提高沥青的抗老化性能,一种方法是向沥青中掺加炭黑、橡胶粉、纳米材料等抗氧化剂以制备改性沥青,另一种方法是直接就选用具有较好的抗老化性能的沥青,其中后一种方法的成本相对较低,且能确保沥青具有较好的存储稳定性。因此,准确评价沥青的抗老化性能并选择具有较好的抗老化性能的沥青用于沥青路面的工程实践中,这将削弱沥青路面使用阶段的老化效应,预防沥青路面的早期损坏,提升沥青路面整体质量,减少沥青路面全寿命周期建设成本。
评价沥青抗老化性能的常用研究方法是对沥青进行旋转薄膜烘箱老化(RTFOT)和压力老化(PAV)后,采用傅里叶红外光谱试验直接测试沥青老化后含氧官能团的变化,或者进行沥青三大指标试验、黏度试验等物理试验间接反映沥青的抗老化性能。然而,这两种试验测试方法均有不足,前一种方法能直接表征沥青的老化程度,但无法精确量化不同沥青的含氧官能团数量的差异;后一种方法虽然能对沥青老化后性能进行量化,但由于这些试验指标并不直接等价于沥青的老化程度,不同试验指标间往往会出现相互矛盾的现象。
与试验方法不同的是,采用基于反应力场的分子动力学模拟方法评价沥青抗老化性能具有如下优点:1)能通过分析沥青老化后含氧官能团的生成以直接反映沥青抗老化性能;2)能对含氧官能团数量进行精确量化;3)这类分子动力学模拟采用的力场参数是根据第一性原理计算的结果优化得到,因此模拟结构反映沥青分子本身的物理化学属性,不受试验条件干扰;4)这种分子动力学模拟方法的计算结果准确性接近量子计算的计算结果,但计算量远小于后者。
发明内容
针对以上问题,本发明提出基于分子动力学模拟的沥青抗老化性能评价方法,本发明的目的在于克服现有试验方法评价沥青抗老化性能研究效果的不足,通过开展基于反应力场的分子动力学模拟,提出一种理论依据可靠、结果简洁明确的方法评价沥青的抗老化性能。
本发明提供基于分子动力学模拟的沥青抗老化性能评价方法,包括如下步骤:
(1)进行沥青组分分离试验将沥青分离成沥青质、胶质、芳香分及饱和分四种组分,并对四种组分进行称重,计算各自所占重量比;
(2)对沥青各组分进行凝胶色谱试验、傅里叶红外光谱试验和元素分析试验,测定所述沥青各组分的相对分子质量、官能团分布和元素组成,以此为依据构建所述沥青各组分的代表性分子,其特征在于每种组分的代表性分子可以是一个,也可以是多个;
(3)根据步骤(1)和步骤(2)的信息,计算建立一个沥青分子模型所需要的每一种代表性分子的数量,并计算出沥青分子模型的摩尔质量;
(4)采用开源软件LAMMPS,将所述沥青代表性分子按步骤(3)计算出的数量与一定数量的氧气分子混合,并按一定的密度构建成一个具有周期性边界条件的立方体盒子,作为沥青老化反应体系以进行老化反应模拟;
(5)采用LAMMPS软件对所述沥青老化反应体系进行模拟,得到老化沥青分子模型,模拟力场采用ReaxFF反应力场,具体步骤是:先进行能量最小化,然后在NVT系综和较低温度下进行较短时长的平衡性模拟,最后在NVT系综和较高的温度下进行较长时长的反应性模拟,该过程中采用LAMMPS软件的fix reax/c/bonds命令按一定频率读取老化反应体系内分子化学键生成与断裂的信息,将结果存储于bonds.reaxc文件中;
(6)采用数据处理软件从所述bonds.reaxc文件中提取羰基、羟基和亚砜基所对应的“C-O”、“H-O”和“S-O”化学键的数量,根据这三种化学键数量计算出模拟得到的老化沥青分子模型中所有氧元素的摩尔质量,以此除以步骤(3)中的沥青分子模型的摩尔质量,即得到老化沥青分子模型的含氧量;
(7)绘制老化反应模拟过程中所述三种化学键数量和老化沥青分子模型的含氧量随老化反应时间的变化图;
(8)选取不同产地、不同品牌或不同规格的沥青,重复步骤(1)-(7),得到它们各自对应的老化特征数据,包括所述三种化学键数量和老化沥青分子模型的含氧量,以此为依据分析和评价所述不同沥青的抗老化性能。
作为本发明进一步改进,所述步骤(5)中,在NVT系综和较高的温度下进行较长时长的反应性模拟时,如果所采用的温度过高,沥青老化反应模拟结果会失真,如果所采用的温度过低,模拟所需时长会大幅增加,1000~1200K的模拟温度是比较适宜的。
本发明具有以下效果:
本发明采用基于反应力场的分子动力学模拟方法,通过精确量化沥青老化后含氧官能团的数量,直接反映沥青抗老化性能,过程规范,结果可靠,克服了实验研究方法无法准确量化或者不能直接反映沥青抗老化性能的不足之处;并且本发明采用的分子动力学模拟方法所用力场参数根据第一性原理计算的结果优化得到,反映了沥青在抗老化方面本身的物理化学属性,不受试验条件的干扰。
附图说明
为了更清楚地说明本发明的技术方案和实施例,本发明进行了附图。
图1:沥青四组分的代表性分子;
图2:沥青老化反应体系;
图3:“C-O”、“H-O”和“S-O”三种化学键随老化反应时间的变化图;
图4:沥青含氧量随老化反应时间的变化图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
本发明提出基于分子动力学模拟的沥青抗老化性能评价方法,本发明的目的在于克服现有试验方法评价沥青抗老化性能研究效果的不足,通过开展基于反应力场的分子动力学模拟,提出一种理论依据可靠、结果简洁明确的方法评价沥青的抗老化性能。本发明将能帮助道路工程师方便快捷地优选出抗老化性能优异的沥青。
本发明基于分子动力学模拟对沥青抗老化性能进行评价,过程如下:
(1)选取某一批次克拉玛依70#基质沥青,记为沥青A,采用溶剂沉淀及色谱柱法将沥青A分离成沥青质、胶质、芳香分及饱和分四种组分;对四种组分进行称重,计算得沥青质、胶质、芳香分及饱和分所占重量比分别为:17.7%、41.4%、28.9%、12.0%。
(2)通过凝胶色谱试验、傅里叶红外光谱试验和元素分析试验测定所述沥青各组分的相对分子质量、官能团分布和元素组成,以此为依据构建所述沥青各组分的代表性分子,每种组分各有一个代表性分子,如附图1所示。
(3)根据步骤(1)中沥青四组分各自所占重量比及步骤(2)所构建的沥青各组分的代表性分子,计算出建立一个沥青分子模型需要沥青质、胶质、芳香分及饱和分四种组分的代表性分子的数量分别为8、25、30、10。
(4)采用LAMMPS软件建立沥青的老化反应体系,首先建立一个具有周期性边界条件的立方体盒子,将8个沥青质、25个胶质、30个芳香分、10个饱和分的代表性分子,以及4025个氧气分子导入其中,其中4025个氧气分子恰好能将体系内沥青的代表性分子完全氧化,体系处于氧气充足状态。体系密度设置为0.1g/cm3,最终得到一个边长为
Figure BDA0002572861120000031
的立方体老化反应体系,如附图2所示。其中深色分子为沥青质代表性分子,浅色分子为氧气分子。
(5)采用LAMMPS软件对所述沥青老化反应体系进行老化反应模拟,模拟采用ReaxFF反应力场,力场参数根据LAMMPS软件自带的ffield.reax.rdx力场参数确定。具体步骤如下:首先对所述老化反应体系进行能量最小化,然后在NVT系综和298K温度下进行200ps的模拟以使体系达到动力学平衡,最后在NVT系综和1200K温度下进行2000ps的模拟以预测沥青的老化反应。
(6)在步骤(5)中采用LAMMPS软件的fix reax/c/bonds命令,按1次/100ps的频率读取老化反应模拟过程中沥青代表性分子的化学键生成与断裂的信息,读取结果存储于命名为“asphalt.bonds.reaxc”的bonds.reaxc文件中。
(7)采用Matlab软件,对“asphalt.bonds.reaxc”文件进行分析,提取羰基、羟基和亚砜基所对应的“C-O”、“H-O”和“S-O”化学键数量,再根据三种化学键数量计算沥青的含氧量。
(8)选取某一批次SK 70#基质沥青,记为沥青B,重复步骤(1)-(7)。
(9)绘制A、B两种沥青老化反应模拟过程中所述三种化学键数量和沥青含氧量随老化反应时间的变化图,分别如附图3和附图4所示。
(10)根据附图3和附图4,对A、B两种沥青的抗老化性能进行评价:在老化反应模拟过程中,沥青A生成的“C-O”、“H-O”和“S-O”化学键数量和含氧量均高于沥青B,因此表明在相同条件下沥青A的老化程度高,沥青B的抗老化性能优于沥青A,即该批次SK 70#基质沥青抗老化性能优于该批次克拉玛依70#基质沥青。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作任何其他形式的限制,而依据本发明的技术实质所作的任何修改或等同变化,仍属于本发明所要求保护的范围。

Claims (2)

1.基于分子动力学模拟的沥青抗老化性能评价方法,其特征在于,包括如下步骤:
(1)进行沥青组分分离试验将沥青分离成沥青质、胶质、芳香分及饱和分四种组分,并对四种组分进行称重,计算各自所占重量比;
(2)对沥青各组分进行凝胶色谱试验、傅里叶红外光谱试验和元素分析试验,测定所述沥青各组分的相对分子质量、官能团分布和元素组成,以此为依据构建所述沥青各组分的代表性分子,其特征在于每种组分的代表性分子可以是一个,也可以是多个;
(3)根据步骤(1)和步骤(2)的信息,计算建立一个沥青分子模型所需要的每一种代表性分子的数量,并计算出沥青分子模型的摩尔质量;
(4)采用开源软件LAMMPS,将所述沥青代表性分子按步骤(3)计算出的数量与一定数量的氧气分子混合,并按一定的密度构建成一个具有周期性边界条件的立方体盒子,作为沥青老化反应体系以进行老化反应模拟;
(5)采用LAMMPS软件对所述沥青老化反应体系进行模拟,得到老化沥青分子模型,模拟力场采用ReaxFF反应力场,具体步骤是:先进行能量最小化,然后在NVT系综和较低温度下进行较短时长的平衡性模拟,最后在NVT系综和较高的温度下进行较长时长的反应性模拟,该过程中采用LAMMPS软件的fix reax/c/bonds命令按一定频率读取老化反应体系内分子化学键生成与断裂的信息,将结果存储于bonds.reaxc文件中;
(6)采用数据处理软件从所述bonds.reaxc文件中提取羰基、羟基和亚砜基所对应的“C-O”、“H-O”和“S-O”化学键的数量,根据这三种化学键数量计算出模拟得到的老化沥青分子模型中所有氧元素的摩尔质量,以此除以步骤(3)中的沥青分子模型的摩尔质量,即得到老化沥青分子模型的含氧量;
(7)绘制老化反应模拟过程中所述三种化学键数量和老化沥青分子模型的含氧量随老化反应时间的变化图;
(8)选取不同产地、不同品牌或不同规格的沥青,重复步骤(1)-(7),得到它们各自对应的老化特征数据,包括所述三种化学键数量和老化沥青分子模型的含氧量,以此为依据分析和评价所述不同沥青的抗老化性能。
2.根据权利要求1中所述基于分子动力学模拟的沥青抗老化性能评价方法,其特征在于;所述步骤(5)中,在NVT系综和较高的温度下进行较长时长的反应性模拟时,如果所采用的温度过高,沥青老化反应模拟结果会失真,如果所采用的温度过低,模拟所需时长会大幅增加,1000~1200K的模拟温度是比较适宜的。
CN202010645196.5A 2020-07-07 2020-07-07 基于分子动力学模拟的沥青抗老化性能评价方法 Pending CN111721693A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010645196.5A CN111721693A (zh) 2020-07-07 2020-07-07 基于分子动力学模拟的沥青抗老化性能评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010645196.5A CN111721693A (zh) 2020-07-07 2020-07-07 基于分子动力学模拟的沥青抗老化性能评价方法

Publications (1)

Publication Number Publication Date
CN111721693A true CN111721693A (zh) 2020-09-29

Family

ID=72573113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010645196.5A Pending CN111721693A (zh) 2020-07-07 2020-07-07 基于分子动力学模拟的沥青抗老化性能评价方法

Country Status (1)

Country Link
CN (1) CN111721693A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112179936A (zh) * 2020-10-12 2021-01-05 南京林业大学 纳米氧化石墨烯改性沥青制备及抗老化性能评价的方法
CN112365936A (zh) * 2020-10-21 2021-02-12 西安理工大学 原子氧环境下沥青氧化老化的分子动力学研究方法
CN112415180A (zh) * 2020-12-10 2021-02-26 广西交科集团有限公司 一种基于灰关联分析的橡胶沥青用基质沥青的选择方法
CN113834839A (zh) * 2021-08-20 2021-12-24 北京工业大学 一种基质沥青及其组分的平均分子结构表征方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105468840A (zh) * 2015-11-20 2016-04-06 长安大学 基于分子动力学的沥青质与二氧化硅界面能评价方法
CN109872778A (zh) * 2019-01-30 2019-06-11 东南大学 一种再生剂在废旧沥青混合料中扩散程度的评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105468840A (zh) * 2015-11-20 2016-04-06 长安大学 基于分子动力学的沥青质与二氧化硅界面能评价方法
CN109872778A (zh) * 2019-01-30 2019-06-11 东南大学 一种再生剂在废旧沥青混合料中扩散程度的评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
DONGLIANG HU 等: "Modeling the Oxidative Aging Kinetics and Pathways of Asphalt: A ReaxFF Molecular Dynamics Study", 《ENERGY FULES》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112179936A (zh) * 2020-10-12 2021-01-05 南京林业大学 纳米氧化石墨烯改性沥青制备及抗老化性能评价的方法
CN112365936A (zh) * 2020-10-21 2021-02-12 西安理工大学 原子氧环境下沥青氧化老化的分子动力学研究方法
CN112415180A (zh) * 2020-12-10 2021-02-26 广西交科集团有限公司 一种基于灰关联分析的橡胶沥青用基质沥青的选择方法
CN113834839A (zh) * 2021-08-20 2021-12-24 北京工业大学 一种基质沥青及其组分的平均分子结构表征方法

Similar Documents

Publication Publication Date Title
CN111721693A (zh) 基于分子动力学模拟的沥青抗老化性能评价方法
Yan et al. Chemical and rheological evaluation of aging properties of high content SBS polymer modified asphalt
Wang et al. Experimental investigation on aging characteristics of asphalt based on rheological properties
Wahhab et al. Prediction of asphalt rheological properties using HP-GPC
Qian et al. Simulation of the field aging of asphalt binders in different reclaimed asphalt pavement (RAP) materials in Hong Kong through laboratory tests
Lin et al. Understanding the low temperature properties of Terminal Blend hybrid asphalt through chemical and thermal analysis methods
Dinegdae et al. Mechanics-based top-down fatigue cracking initiation prediction framework for asphalt pavements
Singh et al. Influence of particle and interaction effects of different sizes of crumb rubber on rheological performance parameters of binders
CN104297200B (zh) 一种红外光谱结合高温模拟蒸馏技术鉴别沥青品牌的方法
CN112034064A (zh) 一种基于分子动力学模拟的沥青老化机理分析方法
Diab et al. Investigating the mechanisms of rubber, styrene-butadiene-styrene and ethylene-vinyl acetate in asphalt binder based on rheological and distress-related tests
CN107421844B (zh) 基于热分析-红外联用的沥青各组分短期老化模拟方法
Shi et al. Nonlinear modelling of selected micro-and macro-properties of weathered asphalt mixtures
Ding et al. Investigation of aging behavior of asphalt under multiple environmental conditions
CN104777288B (zh) 一种利用dsr测定sbs改性沥青中sbs含量的方法
CN109765128B (zh) 基于耗散能的沥青抗老化性能评价方法
Zhao et al. Short-term aging performance evaluation of asphalt based on principal component and cluster analysis
Tang et al. Molecular dynamics simulation and experimental analysis on fluidity improvement of liquid rubber modified asphalt binder
Asi et al. Performance modeling of Arabian asphalt using HP-GPC
CN110411838A (zh) 一种通过应力响应分析评价沥青温度敏感性的方法
Hossain et al. An FT-IRS Approach for the Evaluation of SBS Degradation in SBS Modified Asphalt Binder
CN112415180B (zh) 一种基于灰关联分析的橡胶沥青用基质沥青的选择方法
CN114927170A (zh) 一种基于分子动力学模拟的环氧沥青相容性评价方法
Li et al. Influence of thermo-oxidative aging on the properties of SBS-modified asphalt binders
Hossain et al. Degradation of SBS Polymer during Laboratory Aging of Asphalt Mixture

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200929

RJ01 Rejection of invention patent application after publication